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ABSTRACT 


This  thesis  explores  the  physical  and  mathematical  limitations  of  two  common  attitude 
control  systems:  one  based  on  reaction  wheels  and  another  based  on  control  moment 
gyroscopes  (CMGs).  The  dynamics  are  derived  from  first  principles  and  control 
algorithms  for  achieving  maximum  reaction  wheel  potential  are  discussed.  The  shaped 
eigenaxis  input  is  utilized  to  establish  baseline  maneuver  performance.  A  time-optimal 
shaped  input  is  introduced  and  implemented  in  a  feedback  setting,  subject  to  the 
limitations  of  the  Moore-Penrose  pseudo-inverse  control  allocation.  Finally,  a  feed¬ 
forward  plus  feedback  controller  is  introduced  to  implement  the  time-optimal  torque 
inputs  directly  to  the  reaction  wheels.  This  obviates  the  need  for  the  pseudo-inverse 
control  allocation,  and  therefore  exploits  the  total  capacities  of  both  the  reaction  wheel 
momentum  envelope  and  torque  envelope.  These  reaction  wheel  control  approaches  are 
compared  with  CMG  performance  to  establish  spacecraft  size  and  slew  parameters  that 
make  the  use  of  reaction  wheels  a  reasonable  choice. 
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I.  INTRODUCTION 


A.  MOTIVATION 

Spacecraft  in  low  earth  orbit  (LEO)  traverse  around  the  earth  at  very  high  rates. 
For  instance,  a  spacecraft  in  a  circular  orbit  at  200  nautical  miles  altitude  travels  at  speeds 
in  excess  of  17,000  mph.  This  equates  to  approximately  4  degrees  of  the  orbit  traveled 
every  minute,  completing  a  full  orbit  in  90  minutes.  From  the  point  of  view  of  an  imaging 
satellite,  the  spacecraft  may  only  be  able  to  see  a  ground  target  for  3-4  minutes,  but  this 
amount  of  time  requires  the  orbit  to  pass  directly  over  the  target.  Because  much 
resolution  is  lost  when  taking  low  angle  photographs,  only  the  20  degrees  off-nadir  is 
conventionally  used  [1,  Table  1-1],  [2].  The  off-nadir  constraint  significantly  reduces  the 
amount  of  time  the  target  is  within  the  line  of  sight.  Depending  on  how  far  left  or  right  of 
the  orbit  track  the  target  is,  the  satellite  may  only  have  seconds  to  obtain  an  image.  The 
problem  is  further  complicated  when  multiple  targets  he  along  that  small  target  zone. 

This  space-based  problem  is  even  passed  forward  to  the  military  ground 
commander.  The  battlefield  owner  requires  timely  battlespace  situational  awareness  for 
his  mission  but  may  not  have  the  state-sponsored  weight  in  order  to  gain  access  to 
national  assets  within  a  timely  manner.  The  Army’s  new  imaging  satellite  concept, 
Kestrel  Eye,  is  a  new  approach  to  fill  this  void.  The  intent  is  to  put  many  small, 
inexpensive  satellites  into  orbit  to  obtain  images  and  transmit  them  to  the  battlespace 
owner  more  quickly.  These  satellites  must  be  agile  to  get  multiple  images  per  target  zone 
and  responsive  to  last-minute  request  changes.  Therefore,  conventional  slew  rate 
limitation  of  about  2° /sec  must  be  overcome  [2].  The  mechanisms  that  induce  rotation 
(momentum  exchange  devices)  must  therefore  be  used  to  their  highest  potential.  In  this 
way,  the  spacecraft  can  move  from  target  to  target  as  quickly  as  possible  to  minimize 
otherwise  unproductive  slew  time.  This  thesis  explores  techniques  to  maximize  actuator 
performance  and  minimize  slew  time  without  any  hardware  changes. 
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B.  BACKGROUND 

Control  engineers  typically  employ  two  different  types  of  electro-mechanical 
momentum  exchange  devices  (MEDs)  to  perform  spacecraft  attitude  control.  The  first  is 
reaction  wheels,  which  cause  a  reactionary  body  torque  and  momentum  by  counter¬ 
rotating  a  small  rotor.  Reaction  wheels  are  simple  to  control  and  have  a  high  momentum 
capacity,  which  translates  to  high  angular  velocity  rates  for  the  spacecraft.  The  second 
MED  is  a  control  moment  gyroscope  (CMG),  which  houses  a  spinning  wheel  within  an 
intermediate  gimbal  frame.  The  gimbal  frame  is  allowed  to  rotate,  changing  the 
orientation  of  the  wheel’s  momentum  vector.  The  gimbal  torque  input  is  in  turn  amplified 
(as  the  product  of  the  gimbal  rate  and  the  magnitude  of  the  momentum  vector).  The  high 
torque  output  of  CMGs  is  complicated  by  mechanical  complexity  of  multiple  rotating 
frames  and  the  presence  of  zero-output  singularities  that  must  be  avoided  [3].  The 
elevated  torque,  and  therefore  higher  acceleration  of  the  spacecraft  body,  has  led  CMGs 
to  become  the  de  facto  solution  to  generate  agile  spacecraft.  Although  reaction  wheels  are 
more  prevalent  [4]  (Appendix  B),  cost  less,  are  mechanically  simpler,  weigh  less 
(Appendix  B),  and  are  easier  to  control,  the  torque  output  of  CMGs  simply  cannot  be 
matched  by  reaction  wheels. 

A  spacecraft  slew  is  the  reorientation  from  one  orientation  in  three  dimensional 
space  to  another.  For  a  basic  slew,  the  initial  and  final  spacecraft  orientations  are  given  in 
tenns  of  their  quaternion  reference  points.  To  reconcile  the  difference,  the  most  basic 
approach  is  to  treat  the  difference  in  quaternion  vectors  as  a  single  arc.  The  arc  of 
minimum  distance  rotates  about  a  single  axis,  known  as  the  Euler  axis  [5,  p.  155]  or 
eigenaxis  [6,  pp.  329-338],  A  basic  approach  of  control  therefore  is  to  treat  the  three- 
dimensional  reorientation  problem  as  a  one-dimensional  single  angle  slew  along  this  arc, 
about  the  eigenaxis.  The  eigenaxis  attitude  control  maneuver  requires  a  single  quaternion 
to  be  called  for  and  error  is  regulated  to  zero  as  the  spacecraft  gets  closer  to  the  desired 
quaternion.  This  maneuver  can  be  done  as  a  step  input,  where  the  error  is  constantly 
corrected  for  by  the  controller  [7]-[9].  The  step  input  approach  is  simple  to  implement 
and  stability  is  dependent  upon  only  a  few  gains.  Thus,  the  approach  is  popular  in 
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spacecraft  attitude  control  [10,  pp.  212-267]  [5,  pp.  152-194]  [6,  pp.  403-452]  [11, 
pp.  351-430], 

One  of  the  challenges  of  the  step  input  approach  is  gain  selection.  The  gains 
designed  for  large  maneuvers  may  act  slowly  for  small  maneuvers,  causing  poor 
performance.  Another  approach  is  to  use  a  shaped  trajectory.  Instead  of  passing  a  single 
final  quaternion  to  the  controller,  a  quaternion  path  for  the  maneuver  is  calculated  and  a 
set  of  discrete  intermediate  quaternions  is  passed  as  a  target  trajectory.  Here,  the  gains  are 
tuned  to  give  acceptable  performance  and  their  values  are  agnostic  to  the  slew  size. 
Creamer  et  al.  [12]  developed  a  shaped  bang-coast-bang  maneuver  about  the  eigenaxis, 
used  on  board  the  Clementine  spacecraft.  Time-optimal  quaternion  shaping  was  used  on 
the  TRACE  spacecraft  [13].  Additionally,  [14]  suggests  shaping  both  trajectory  and  body 
rates  and  using  these  as  inputs  to  the  spacecraft  attitude  control  (ACS)  system. 

A  typical  ACS  uses  four  or  more  wheels  for  redundancy.  However,  the  use  of  four 
or  more  wheels  presents  a  control  allocation  problem.  This  is  caused  when  torque  or 
momentum  required  in  the  body  frame  must  be  produced  by  the  redundant  set  of  MEDs. 
A  common  method  of  torque  and  momentum  allocation  is  to  use  the  Moore -Penrose 
pseudo-inverse,  which  provides  a  least  squares  (L2)  solution  [5,  p.  169]  [6,  p.  440]  [11, 
p.  416]  [15].  According  to  this  approach,  the  full  magnitude  of  the  requested  body  torque 
may  not  be  available  in  every  direction.  The  minimum  norm  solution  is  devised, 
providing  a  unique  allocation  of  wheel  torques  [16,  pp.  49-57].  The  pseudo-inverse  does 
not  take  physical  wheel  limitations  into  account  and  may  request  more  torque  than  is 
available.  Once  any  wheel  is  saturated,  then  the  remaining  allocated  torque  is  simply  lost. 
This  is  because  to  remain  on  eigenaxis,  the  integrity  of  the  least  squares  solution  must  be 
maintained.  Thereby,  proportional  saturation  across  all  wheels  is  necessary.  Because  of 
the  nature  of  the  Moore  Penrose  pseudo-inverse,  maximum  momentum  and  torque  will 
be  unavailable  about  certain  axes.  The  pseudo-inverse  is  described  within  [15]  in  great 
detail,  providing  detailed  illustrations  of  the  L2  solution  in  three-dimensional  space  given 
various  reaction  wheel  array  configurations. 

A  way  to  bypass  the  L2  limitation  is  to  operate  the  ACS  in  the  open-loop 

(feedforward).  One  approach  is  to  calculate  the  optimal  torque  solution  and  command  the 
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wheels  directly.  This  is  far  from  conventional  thinking,  because  the  individual  wheel 
torque  solution  is  not  generally  available.  The  shape  of  the  torque  signal  can  be  derived 
over  the  course  of  the  maneuver.  This  feedforward  concept  has  been  around  for  decades 
[  17]— [19],  but  the  approach  is  susceptible  to  instability  due  to  external  disturbances  and 
uncertainties  of  the  spacecraft  inertia  or  initial  states  [9],  For  the  solution  to  work 
properly,  all  states  and  parameters  must  be  established  exactly  prior  to  reorientation,  with 
no  uncertainties  throughout  the  maneuver.  This  ideal  scenario  is  unlikely,  so  a  merger  of 
open-loop  and  closed-loop  (feedback)  must  be  used  to  improve  stability  and  perfonnance. 
Closed-loop  implementation  of  open-loop  maneuvers  was  flown  on  TRACE  [20]  with 
great  success,  but  the  performance  gains  were  somewhat  limited  by  the  spacecraft 
preexisting  software.  Based  on  the  literature  available,  it  appears  that  industry  may  not  be 
utilizing  reaction  wheels  to  their  full  potential,  fostering  a  greater  disparity  between  the 
capabilities  of  reaction  wheels  and  CMGs. 

C.  LIMITATIONS  ON  AGILITY 

In  practice,  a  spacecraft  slew  rate  is  restricted  for  a  variety  of  reasons.  The  most 
obvious  limitations  are  due  to  the  physical  limitations  of  the  MED  array.  For  a  reaction 
wheel  system,  an  individual  wheel  is  limited  on  how  fast  it  can  spin,  thus  the  momentum 
it  can  produce.  Positive  rotation  at  maximum  angular  velocity  results  in  a  maximum 
momentum  vector,  while  counter-rotating  results  in  a  momentum  vector  of  equal 
magnitude  and  opposite  direction.  Given  a  skewed  array  of  several  reaction  wheels,  the 
vector  sum  of  the  maximum  momentum  vectors  in  any  arbitrary  direction  creates  the 
momentum  envelope  (Chapter  IV).  The  shape  of  the  reaction  wheel  momentum  envelope 
is  a  polyhedron  [15],  involving  sharp  edges  and  vertices.  This  signifies  that  more 
momentum  is  available  in  certain  directions  than  others.  This  momentum  envelope 
produces  the  physical  rate  limit  of  the  system  about  any  given  axis. 

Many  attitude  control  systems  involve  restricting  the  axis  of  rotation  to  the 
eigenaxis,  thereby  minimizing  the  angle  of  rotation.  Since  the  momentum  space  is  not 
spherically  uniform,  some  axes  allow  faster  rotation  rates.  This  causes  additional 
problems  with  slew  planning.  Not  only  are  two  quaternions  in  space  given,  but  there  is 
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also  an  eigenaxis-specific  slew  rate.  Not  to  mention  the  increase  in  eigenaxis-specific 
slew  acceleration  (derived  from  the  torque  envelope,  Chapter  IV)  to  be  concerned  about. 

When  operators  want  an  imaging  satellite  to  obtain  an  image,  the  rate  and 
acceleration  variations  over  47t  steradians  require  them  to  focus  on  the  temporal  aspect  of 
the  maneuver.  Typical  questions  are:  Given  the  slew  angle  and  specific  eigenaxis,  will 
the  spacecraft  arrive  at  the  desired  quaternion  at  the  correct  time?  Is  there  another  axis 
the  spacecraft  can  slew  about  more  quickly  to  image  something  else?  Even  with 
automation,  trajectory  planning  becomes  a  time-consuming  task.  It  may  be  far  easier 
from  an  operational  point  of  view  to  operate  within  a  single  capability,  even  if  it  means  a 
loss  of  perfonnance.  This  procedural  limitation  further  reduces  the  operable  reaction 
wheel  momentum  space  polyhedron  to  the  largest  inscribed  sphere,  which  provides 
guaranteed  capability  achievable  in  every  direction.  A  simple  example  is  a  sphere 
inscribed  inside  of  a  cube  (Figure  1).  The  envelope  (cube)  restricts  the  magnitude  of  the 
vector  to  1.0  in  some  directions,  while  the  maximum  magnitude  is  1.73  at  the  vertices. 
This  amounts  to  a  non-uniform  momentum  space,  offering  values  of  momentum  about 
some  axes  that  are  considerably  less  than  others.  If  operating  at  the  envelope,  the  motion 
is  difficult  to  predict  without  correlating  the  specific  eigenaxis  with  the  available 
momentum.  Therefore,  the  engineer  might  design  to  the  lowest  magnitude  (sphere)  and 
provide  this  momentum  limit  to  the  planners  as  their  planning  constraint  for  travel  about 
any  axis.  The  same  logic  can  be  applied  for  torque  limitations,  but  the  main  point  here  is 
that  some  performance  is  lost  due  to  the  need  to  simplify  other  aspects  of  spacecraft 
operations. 
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Max  Sphere  =  1.0 
Max  Envelope  =  1.73 


Figure  1 .  Sphere  Inscribed  within  a  Cube 

D.  OBJECTIVES  AND  SCOPE  OF  THIS  THESIS 

This  thesis  is  an  attempt  to  close  some  of  the  performance  gaps  between  reaction 
wheels  and  CMGs,  with  the  hope  of  identifying  some  key  regions  where  reaction  wheels 
may  actually  outperform  CMGs.  The  true  capabilities  of  reaction  wheels  are  therefore 
explored  to  detennine  whether  or  not  such  an  improvement  is  even  feasible.  The  shaped 
trajectory  used  on  the  Clementine  spacecraft  [12]  is  developed  (Chapter  III)  and  analyzed 
(Chapter  V).  This  maneuver  becomes  the  basis  for  comparison  of  two  different  reaction 
wheel  control  algorithms  that  were  developed  as  part  of  this  thesis. 

The  main  development  is  the  time-optimal  maneuver,  which  merges  shaped 
torque  feedforward  with  shaped  quaternion  and  spacecraft  body  rate  feedback  (Chapter 
V).  This  approach  avoids  the  mathematical  limitations  of  pseudo-inverse  control 
allocation  by  applying  torque  signals  directly  to  the  reaction  wheels.  The  caveat  of  using 
feedforward  plus  feedback  is  that  a  margin  must  be  built  in  to  account  for  uncertainties. 
This  adds  to  the  body  of  knowledge  because  the  use  of  optimal  quaternion  and  body  rate 
trajectories  has  the  ability  to  correct  for  spacecraft  parameter  uncertainties  and 
disturbances.  An  approach  for  adding  system  margin  based  on  a  level  of  uncertainty  is 
discussed.  Optimal  solutions  are  solved  using  DIDO  [21]  and  then  propagated  through  a 
computer  simulation  of  the  spacecraft  dynamics  to  verily  their  feasibility. 
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The  second  reaction  wheel  enhancement  is  an  extension  of  the  time-optimal 
trajectory,  but  limited  to  the  pseudo-inverse  space.  This  method  again  breaks  away  from 
the  eigenaxis,  but  operates  entirety  within  the  confines  of  the  L2  volume.  This  allows  the 
control  engineer  to  upload  the  time-optimal  quaternion  trajectory  instead  of  simply  the 
eigenaxis  quaternion  trajectory  [12].  Since  the  only  change  is  the  shape  of  the  trajectory, 
it  is  not  very  intrusive  on  the  current  employment  of  a  standard  shaped  quaternion 
trajectory  implementation.  This  method  will  be  introduced  at  the  end  of  Chapter  V. 

To  aid  in  the  setup  of  the  problem,  the  spacecraft,  reaction  wheel,  and  CMG 
dynamics  will  be  derived  from  first  principles.  This  is  done  in  Chapters  II,  IV,  and  VI, 
respectively.  The  equations  of  motion  are  then  implemented  within  Simulink,  which  will 
act  as  a  computer  simulation  of  the  actual  nonlinear  spacecraft  dynamics.  Various  control 
algorithms  will  be  applied  to  the  same  simulation  to  retain  the  structure  of  the  model. 
Either  the  reaction  wheel  or  CMG  simulations  can  be  selected  as  the  current  MED  device 
within  the  model.  Simulink  also  offers  a  graphical  depiction  of  mathematical  changes  to 
the  full  dynamic  setup,  whereas  changes  can  be  hidden  within  loops  if  utilizing 
MATLAB  alone. 

The  sphere  concept  is  applied  to  the  time-optimal  case.  For  a  series  of  time- 
optimal  maneuvers  between  any  two  quaternions,  it  is  shown  how  the  minimum 
guaranteed  momentum  and  torque  can  be  derived  for  any  direction  of  travel  about  any 
effective  eigenaxis.  This  adds  the  ability  for  the  planner  to  devise  a  single  slew  rate  over 
any  rotation  angle,  because  the  sphere  remains  independent  of  rotation  axis  and  angle 
size.  The  effective  sphere  increases  in  size  twice  from  the  smallest  inscribed  sphere,  once 
for  each  of  the  two  optimal  control  techniques  mentioned  above.  Therefore,  two  levels  of 
ACS  improvement  are  given  as  part  of  this  work.  Chapter  V  presents  the  concept  of 
effective  eigenaxis  for  an  optimal  maneuver  and  the  details  of  the  derivation. 

A  new  technique  for  visualizing  slews  is  introduced  in  Chapter  VII.  This  method 
utilizes  the  nonnalized  momentum  and  torque  envelopes  as  the  surface  to  display 
momentum  and  torque  vectors  throughout  the  maneuver.  In  the  literature  review, 
eigenaxis  and  quaternion  traces  were  found,  which  only  accounts  for  the  path  traveled. 

This  new  visualization  tool  allows  the  behavior  of  the  rate  and  acceleration  to  be 
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observed  with  respect  to  the  available  momentum  and  torque  capacity.  The  tool 
demonstrates  fundamental  features  of  the  optimal  control  maneuvers. 

Finally,  the  perfonnance  of  all  three  reaction  wheel  control  approaches  are 
compared  against  that  of  a  like-sized  CMG  in  Chapter  VII.  Data  from  the  market  research 
from  over  75  reaction  wheels  and  CMGs  is  used  to  generate  trends  and  differences 
between  these  different  systems.  A  simple  method  of  comparing  reaction  wheels  and 
CMGs  is  developed.  Momentum  and  torque  capacities  are  compared  with  slew  angle 
sizes  to  detennine  regions  where  reaction  wheels  outperform  CMGs  under  the  various 
control  algorithms  tested.  Concluding  remarks  and  suggestions  for  future  work  are  given 
in  Chapter  VIII. 
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II.  SPACECRAFT  KINEMATICS  AND  DYNAMICS 


The  purpose  of  this  chapter  is  to  develop  the  spacecraft  kinematics  and  dynamics. 
Spacecraft  attitude  kinematics  will  be  modeled  using  quaternions.  Spacecraft  dynamics 
will  be  derived  from  first  principles,  showing  the  dynamic  relationship  between  the 
momentum  exchange  device  (MED)  action  and  the  spacecraft  body  reaction.  The  MED 
assembly  (subscript  medA  in  this  chapter)  is  any  generic  arrangement  of  momentum 
exchange  devices.  It  can  be  3,  4,  or  n  reaction  wheels,  n  CMGs,  or  any  combination  of 
thereof. 

A.  MODELING  A  SPACECRAFT 

The  spacecraft,  MED  assembly,  and  feedback  can  be  modeled  within  Simulink, 
shown  in  Figure  2. 
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Figure  2.  Spacecraft  Simulation  Model 


Given  a  desired  (commanded)  quaternion  and  a  current  quaternion  state,  a 
quaternion  error  vector  is  produced.  The  controller  attempts  to  drive  that  error  to  zero  by 
producing  an  appropriate  torque  on  the  spacecraft.  The  MED  Assembly  attempts  to 
produce  that  torque.  The  Spacecraft  Dynamics  model  the  spacecraft  response,  with  torque 
disturbances  (Td)  possible,  but  these  are  not  considered  in  this  thesis.  The  spacecraft 
kinematics  detennines  the  current  quaternion  state.  In  this  thesis,  orbital  motion  is  not 
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considered,  so  N(0°  =  0 ,  but  it  is  shown  to  clarify  the  relationship  between  ncqb  and 
°(0B .  The  rate  of  the  body  with  respect  to  the  orbit  is  utilized  within  the  spacecraft 
kinematics.  This  is  detennined  by  subtracting  the  orbital  motion  NC0°  from  v (0 11 . 


N(0B  =  N(0°  +  °G)B 


°(0B  =  N(0B-N(0° 


(1) 


Since  orbital  motion  is  not  considered,  NC0B  =  °(0B .  The  angular  rate  of  the  body 
with  respect  to  inertial  space  is  thus  defined  in  Equation  (2). 


w  = 


co_ 


(2) 


B.  QUATERNION  KINEMATICS 

Quaternion  kinematics  can  be  written  as  a  set  of  four  differential  equations.  The 
first  three  quaternions  can  be  written  as  qm ,  while  q4  is  kept  as  a  separate  dynamic 
equation  for  simplicity,  shown  in  Equation  (3)  [6,  Algorithm  5.76]. 


9i23  =  |(«4  °(Ob-°CObx 

1  (  n  n\T 

^4=-2(  (O  )  qm 


(3) 


where  qm  - 


7 1  0_2  ^3 


.  The  entire  quaternion  vector  can  be  written  as 


q  = 


*123 


^4 


.  Another  thing  to  note  about  quaternions  is  that  the  two-norm  is 


always  equal  to  one.  This  allows  one  to  solve  for  q4 ,  given  qm  .  The  origin  for  this  thesis 
will  be  defined  as#  =  [o  0  0  1  ]  . 


\\q\l  =  Jql2+q22+q32+q/  =1 

ft  2  2  2 

74  =  Vl”7i  ~q2  ~q3 


(4) 
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c. 


MODELING  SPACECRAFT  KINEMATICS 


Equation  (3)  is  modeled  into  the  Spacecraft  Kinematics  Simulink  block,  shown  in 
Figure  3.  The  direction  cosine  matrix  (DCM)  mentioned  is  only  used  to  determine  the 
disturbance  torques  and  is  not  used  here  as  the  disturbances  are  taken  to  be  zero. 


Spacecraft  Kinematics 


dq  =  1/2  (q4  °u)B  -  °wB  x  q) 
dq4  =  1/2  (°uB)T  q 


Create 

DCM 


D.  SPACECRAFT  DYNAMICS 


For  a  rigid  body  spacecraft  with  no  external  torques  applied,  the  total  angular 
momentum  must  be  conserved.  Since  the  total  angular  momentum  h^t  must  be  constant, 
then  a  positive  variation  in  the  MED  assembly  angular  momentum  hhmedA  will  have  an 
opposite  impact  on  spacecraft  angular  momentum  hNlc  and  vice  versa.  All  momenta  are 
written  with  respect  to  the  inertial  frame. 


hl  =  Klc  +  h 


where  \h^t  =  constant 


(5) 
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Angular  momentum  is  equivalent  to  the  product  of  inertia  and  angular  rate 
(. h  =  J(0 ).  Applied  specifically  to  the  spacecraft,  Jslc  represents  the  total  structural 

inertia  of  the  spacecraft.  The  inertia  of  the  MED  assembly  is  included  within  the 
spacecraft  inertia.  All  configuration  changes  of  the  MED  assembly  with  respect  to  the 
relatively  small  spacecraft  body  angular  velocity  will  be  considered  negligible 
(  jmedA  N COB  ~  0  ) .  The  angular  rate  of  the  spacecraft  is  expressed  in  the  body  frame  relative 

to  the  inertial  frame  co  =  NooB . 

Kc  J s/c^slc  s/cbody  'J  >!:/C\  )  ^ 


In  this  thesis,  equations  often  are  written  in  dyadic  fonn  to  ensure  that  notation  of 
the  bases  is  clearly  understood.  Some  derivations  involving  multiple  frames  can  be  easier 
manipulated  (and  understood)  when  utilizing  dyadic  form.  Following  this  notation  and 
assuming  a  diagonal  inertia  tensor,  Equation  (6)  can  be  rewritten  as 


Kc  =  (jsujA  +  istcjlk  +  JsIcJA  )  •  (®*4  +  ®  A  +  0):K  ) 

=  A/cn®  A  (4  *4)  +  jslcu<°yb  1  (4  *  4  )  +  j,lcy<°A  (4  *  4  ) 

+ Jsu:jjA  (4  *  4 ) + jsicj°y  4  (4  •  4 ) + j,/cj°z  4  (4  •  4 ) 
+ Js/cjiA  (4  •  4 ) + A/C,«v4  (4  •  4 ) + A/c33®-4  (4  •  4 ) 


(7) 


Since  bx  ,b2  ,b3  are  orthogonal,  b,  •  bf  =  1  and  bi  •bj  =  0  .  This  leads  to 


Kc  =  jsIc^  A+js/c^A+jsIc^cK  ■ 


(8) 


Sometimes  it  will  be  more  convenient  to  show  an  equation  in  matrix  form, 
Equation  (9),  but  either  form  can  be  used  universally. 


Js/cn 

0 

0 

B 

(»x 

B 

A/q®, 

II 

0 

J  s!  c22 

0 

COy 

- 

J  s/c22  ^ y 

0 

0 

J s/c33 

°>z 

jslc330)z 

(9) 
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Since  angular  momentum  must  be  constant,  the  derivative  with  respect  to  time 
must  be  zero. 


K  =  Ke  +  hN=  o 


(10) 


External  torques  can  be  added  as  disturbance  torques  per  Equation  (1 1). 


(ii) 


Disturbance  examples  in  LEO  are  gravity  gradient,  aerodynamic  loads,  and  solar 
radiation.  These  and  other  external  torques  will  be  considered  negligible  in  this  thesis  and 
will  be  set  to  zero  ( rf,  =  0 ). 


Equation  (10)  reveals  that  the  spacecraft  and  MED  assembly  angular  momentum 
rates  must  be  equivalent  in  magnitude,  opposite  in  direction.  Also,  the  torque  required  to 
manipulate  the  spacecraft  corresponds  to  the  spacecraft  angular  momentum  derivative  in 
the  inertial  frame  h^c  =  T^eq  .  Thus,  the  two  terms  can  be  used  interchangeably  to  express 

the  torque  applied  to  the  spacecraft  by  the  MED  assembly.  Accordingly,  the  torque 
required  of  the  MED  assembly  is  equivalent  to  the  negative  torque  in  the  body  frame, 
shown  in  Equation  (12). 


hN 


=  -h 


sic 


(12) 


To  obtain  the  angular  momentum  rate  with  respect  to  the  reference  frame,  the 
angular  rate  between  the  frames  is  necessary.  Therefore,  the  angular  momentum  rate 
(torque)  with  respect  to  the  inertial  frame  is  equal  to  the  sum  of  the  torque  in  the  body 
frame  and  the  cross  product  between  angular  rate  and  angular  momentum. 


(13) 


The  time  derivative  of  Equation  (6)  is  shown  in  Equation  (14).  The  jslc 
component  is  only  applicable  if  fuel  is  being  expended  during  the  maneuver  or  if  flexible 
elements  are  causing  change.  It  is  assumed  that  any  changes  would  be  very  small; 
therefore  the  spacecraft  can  be  modeled  as  a  rigid  body.  Moreover,  it  is  assumed  that  the 
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inertia  of  the  spacecraft  is  not  changing  over  the  course  of  a  maneuver  and  the  tenn  is 
zero,  shown  in  Equation  (15). 


j'B  —  T  N  ,  T  N-B 

*^dr  —  J  dr  ®  J dr  CO 


j’B  j  N  •  B 

"• s/c  —  J  s/c  ® 


Substituting  Equation  (6)  and  Equation  (15)  into  Equation  (13)  yields 


VN  j  N  -B  .  N  B  .  j  N  B 

h  <dr  ~  J  dr  co  co  x  J  ,  co 


Equation  (17)  shows  the  derivation  in  matrix  fonn. 


Kc  = 


js,jh  +  [jslc33  -  j,c12  )C0z(0y 
is/cJVy  +{jS/cu  ~  hie,,)03.03:. 
JsIcA  +  ( Jslc22  ~  Jslcn  )(0yC0x 


(14) 

(15) 


(16) 


(17) 


The  derivation  may  also  be  carried  out  in  dyadic  fonn: 

Kc  =  jsuJlM  +  j +  7,/c*®  A 

(18) 

+(®. A  +  ®A  +  ®A)  X  (js,cuC0.A  +  jslCi2CDyb2  +  7,/C9©A) 

Kc  =  Acn®A  +  jsujKy  +  Ac*,®  A 

+  {js,c33  -  j,lcn  A  +  (jslcu  ~  AcJ®c®  A  +  (jslc22  ~  llcn  )®y®  A 

=  [jSlcna*+(jslc33-  jslcn)®  Z(°y)k  (19) 

+  (  Ac22  A  +  (j,cn  ~  Js/c33  )(°X(Oz)b2 

+  (jsie33(bz+(jslC22-jslCu)coYcox)b3 


Conespondingly,  by  refening  back  to  Equation  (12),  the  torque  of  the  MED 
assembly  is  related  to  the  spacecraft  torque  per  Equation  (20).  The  MED  portion  is 
written  in  the  inertial  frame  and  includes  the  stored  momentum  in  the  body  frame,  to  be 
later  derived  in  Chapters  IV  and  VI  for  reaction  wheels  and  CMGs,  respectively. 

l-N  r  N  •  B  (V_.fi..  j  N  B 

Kd=-Js,c  (0  -  co  X Jsk  CO  (20) 
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Rewriting  the  angular  momentum  rate  in  Equation  (20)  as  its  corresponding 
torque  yields 


medA 


^.B  N  .  B  v .  t  N  B 

^  CO  XV.  CO 


(21) 


E.  MODELING  SPACECRAFT  DYNAMICS 

As  a  slight  variation  of  Equation  (21),  Equation  (22)  is  modeled  within  the 
Spacecraft  Dynamics  Simulink  block,  shown  in  Figure  4. 


-B  _ N  N  -B  ..  T  N--B 

T  =  — T  -  CO  xj.  CO 

sic  medA 


Note  that  the  disturbance  torques  Td  =  0  .  Thus,  this  term  is  not  included  in  (22). 


Spacecraft  Dynamics 


Figure  4.  Spacecraft  Dynamics  Model 
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III.  EIGENAXIS  REORIENTATION 


While  traveling  at  several  thousand  kilometers  per  second,  a  spacecraft  is 
commanded  to  move  from  one  quaternion  state  to  another.  This  reorientation  time  may 
take  a  second  or  two  for  small  maneuvers  or  tens  of  seconds  to  minutes  for  large 
maneuvers.  For  LEO  Earth-imaging  spacecraft,  this  reorientation  must  occur  swiftly  and 
precisely.  Reorientation  (slew)  time  is  time  lost  as  it  is  time  not  conducting  mission- 
essential  targeting.  Since  spacecraft  can  traverse  key  areas  in  minutes  (or  possibly 
seconds),  slew  time  ultimately  limits  the  number  of  targets  per  pass.  Since  these  multi¬ 
million  or  billion  dollar  spacecraft  are  only  productive  when  they  are  not  slewing,  slew 
time  is  productivity  lost. 

A.  SHORTEST  ANGLE  MANEUVER 

A  simple,  effective  slew  maneuver  is  the  eigenaxis  reorientation  or  eigenaxis 
slew.  The  eigenaxis  slew  is  a  rotation  along  the  shortest  arc  between  two  attitudes.  This 
maneuver  is  used  due  to  its  simplicity,  predictability,  and  reasonably  rapid  speed,  as 
opposed  to  performing  successive  rotations  about  the  individual  body  axes.  Given  three- 
axis  control,  there  are  infinite  paths  that  may  be  taken  to  transition  between  two  attitudes. 
The  eigenaxis  slew  utilizes  the  shortest  of  these  paths.  The  axis  of  rotation  is  known  as 
the  eigenaxis,  whose  direction  remains  constant  and  orthogonal  to  the  direction  of 
angular  travel  throughout  the  maneuver.  Assuming  that  the  spacecraft  can  be  reoriented 
while  maintaining  rotation  about  the  eigenaxis,  the  problem  can  be  simplified  into  a 
single  angle  slew,  shown  in  Figure  5. 
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eigenaxis 
(out  of  page) 


0  rotation  angle 


Figure  5.  Eigenaxis  Slew 


For  a  reorientation  from  one  angle  to  another,  a  direction  cosine  matrix  (DCM) 
B  Ra  (discussed  in  greater  detail  in  Chapter  IV)  can  be  used  to  show  the  necessary 
relative  reorientation  from  reference  A  to  reference  B.  There  is  at  least  one  eigenvector 
associated  with  the  DCM  that  produces  a  unity  eigenvalue  shown  by  Sidi  [5,  Algorithm 


4.1]. 


RAe  =  le 


B 


(23) 


Equation  (23)  represents  the  fact  that  instead  of  breaking  a  DCM  into  a  series  of 
three  consecutive  rotations,  the  reorientation  can  be  created  by  using  a  single  rotation 
about  a  specific  axis  or  eigenvector.  Unity  eigenvalue  occurs  when  rotating  about  the 
principle  eigenvector  or  eigenaxis.  Because  of  this,  the  eigenaxis  slew  is  the  shortest 
angle  maneuver. 

Assume  the  spacecraft’s  orientation  is  originally  described  by  quaternions  such 
that  q0  =  [  r/0  q(]_  q0  q(^  ]'  with  a  requirement  to  slew  to  some  final  quaternion 

state  q ,  =  [  qf  qfi  qh  qfi  ]7  .  Note,  q4  is  the  scalar.  The  quaternion  difference  qd 
can  be  described  by  Equation  (24)  [6,  Algorithm  7.69]. 
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(24) 


9/4 

9/3 

-9/3 

-9/, 

9o, 

<u2 

-9/3 

9/4 

9/, 

-9/3 

9o2 

4d3 

9/2 

-9/ 

9/4 

-9/3 

9o3 

C1f 

9/3 

9/3 

9/4 

9o4 

Equation  (25)  shows  the  relationship  between  the  quaternion  difference,  the 
eigenaxis,  and  the  slew  angle  [5,  Algorithm  4.1 1]. 


qd  =  e  sin 


cLl 


■  cos 


6  2 A 
V  2  j 


(25) 


V  J 


Rearranging  Equation  (25),  the  desired  slew  angle  can  now  be  solved  for. 

6d  =  2cos~'  [q(h ) 


(26) 


The  eigenaxis  can  now  be  solved  as 

1 

C~  sm{edl2)qdm 


(27) 


B.  QUATERNION  FEEDBACK 

Using  an  eigenaxis  slew,  the  reorientation  problem  can  be  simplified  into  a  single 
angle  slew,  like  the  tracing  of  an  arc  of  a  circle  on  a  2D  plane.  This  section  will  explain 
how  the  complicated  spacecraft  dynamics  derived  in  Chapter  II  can  be  modeled  as  a 
linear  double  integrator,  providing  the  spacecraft  travels  precisely  along  this  eigenaxis 
arc.  The  double  integrator  model  will  be  obtained  from  the  eigenaxis  quaternion  feedback 
(EQF)  control  logic  often  used  in  spacecraft  ACS  [6,  Algorithm  7.74]. 

«"  =  -kjs/c  qem  ~°J sic N<0B  +  V  X  Js,c NcoB ,  (28) 

Recall  the  nonlinear  spacecraft  torque  from  Equation  (16).  This  inertial  torque  is 
what  the  EQF  controller  is  attempting  to  control  («"  =  . 
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ii  —N  j  N  -B  .  N  B  .  .  T  N  ,^B 
11  =^slc  =  Jslc  °>  +  °J  XJsIc  ® 


(29) 


Equations  (28)  and  (29)  both  contain  the  same  gyroscopic  coupling  term,  a  term 
not  conducive  to  solving  a  linear  trajectory  along  an  arc.  The  nonlinearity  can  be 
eliminated  and  a  new  control  u’  is  detennined  in  Equation  (30).  This  shows  that  the 
gyroscopic  coupling  term  was  simply  added  to  the  EQF  to  eliminate  this  nonlinearity  in 
the  control. 


«'  =  «"- WxJs/cW 

=  -kJslcqein-cJslcN<oB  (30) 


Since  tbc  =  Js/c  N(bB  from  Equation  (15),  the  control  u  ’  is  effectively  a  torque 
control  in  the  body  frame.  Equation  (30)  can  be  rearranged  as  follows. 


«'  +kJslcqem  +  cJslcNC0B=  0 

JslcN(bB  +  kJslcqein+cJslcNCOB  =  0 


(31) 


The  constant  inertia  can  then  be  divided  out  to  assume  acceleration  as  the  control 
variable. 


u  =  u'/Jslc 


(32) 


u  +  k  Q ,  c  ^  0)B  —  0 

•*^123 

N  B  -  /  .  ~  N  ,-^B 

6)  +  k q,  +c  co  =0 

^123 


(33) 


The  problem  thus  far  is  measured  in  three-dimensional  space.  But  the  purpose  of 
an  EQF  is  to  maintain  the  trajectory  along  the  eigenaxis.  Since  the  arc  is  measured  along 
a  plane,  the  quaternion  trajectory  can  be  measured  approximately  as  a  single  angle,  while 
the  body  rate  can  be  measured  as  a  single  angular  rate,  both  measured  as  an  error.  The 
simplification  is  explained  in  [6,  Example  7.15]. 

eE  =  e-ed  (34) 

u  +k{e-e(l)+c{e-ed)  =  o 
eE  +  k{e-ed)+c(e-ed)= o 
20 


(35) 


In  Equation  (33),  N coB  was  not  considered  an  error  term  because  the  control 
assumes  no  shaped  angular  velocity.  Instead,  it  regulates  it  to  zero.  Therefore,  the  rate 
error  is  equal  to  body  rate  0E=0  =  N coB .  The  fundamental  fonnula  for  a  proportional 
plus  derivative  controller  (PD)  is  shown  as 

0E  +  k0E  +  c0E  =  0  (36) 


Equation  (36)  is  now  written  in  terms  of  a  double  integrator,  where  the  control  is 
equal  to  the  second  derivative  of  the  signal  u  =  dE  =  a . 


u  =  a 

=  —kOE  -cOE 


(37) 


Assume  a  desired  angle  and  rate  are  given.  The  controller  can  now  regulate  the 
acceleration  signal  to  drive  the  error  to  zero.  As  stated  previously,  the  desired  rate  is 
zero,  therefore  it  is  eliminated  from  the  equation.  This  means  that  Equation  (28)  was 
treated  as  a  rest-to-rest  maneuver.  The  double  integrator  in  Equation  (38)  is  implemented 
in  Figure  6. 

u  =  a  =  -k(0-dd)-cco  (38) 


Desired  Input  Controller  Spacecraft 

Input  Error  a  1 


Figure  6.  PD  Control  of  a  Double  Integrator  Model 

Therefore,  the  EQF  controller  is  essentially  just  a  PD  control  system  in  disguise. 
The  gyroscopic  coupling  term  can  be  added  back  in  to  the  PD  controller,  just  to  be 
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subtracted  back  out  by  the  spacecraft  physics.  Therefore,  the  spacecraft  can  be  slewed  as 
a  double  integrator,  as  shown  in  the  simplistic  spacecraft  simulation  in  Figure  6. 


C.  GAIN  DESIGN  FOR  QUATERNION  ERROR  FEEDBACK 

If  the  goal  is  to  perfonn  a  rest-to-rest  reorientation,  0d  may  be  fixed  and  act  as  a 
step  input,  which  moves  from  zero  to  the  desired  angle.  As  with  any  step  input,  the  error 
is  very  large  at  first  and  decreases  significantly  over  time.  Proportional  and  Differential 
Gains  (/c,c)  must  therefore  be  chosen  carefully. 


Underdamped  systems  are  the  quickest  to  respond,  but  must  be  carefully  designed 
so  the  torque  (acceleration)  and  overshoot  do  not  become  too  large.  The  percent 
overshoot  ( Mp )  is  a  measure  of  the  overshoot  compared  to  the  starting/ending  condition 
difference.  For  the  following  example,  10%  overshoot  is  used.  Settling  time  ts  is  the  time 
it  takes  the  system  to  fall  and  remain  within  a  certain  band,  commonly  detennined  by 
imaging  and  targeting  requirements  (e.g.,  2%  of  the  final  or  steady-state  value).  An 
example  of  gain  selection  is  shown  in  (39)  such  that  the  percent  overshoot  and  settling 
time  is  the  selection  criteria  for  damping  ratio  and  natural  frequency  [22,  Algorithms 
4.39,  4.42], 


c 


-ln(Mp  / 100) 
Jn2+ln(Mp/100) 


a ) 


n 


(39) 


The  previously  mentioned  k  and  c  are  the  feedback  gains  for  the  second  order, 
homogeneous  system  shown  in  Equation  (40).  The  same  canonical  expression  is  shown 
in  Equation  (41)  in  terms  of  the  damping  ratio  and  natural  frequency. 


e+ce+ke= 0 

(40) 

e  +  2£cQne  +  co2e  =  0 

(41) 
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Now,  it  can  be  seen  that 


c  =  2£w„ 
k  =  o)n2 


(42) 


Using  the  feedback  control  system  mentioned  in  Section  C,  the  ideal  gains  for  a 
quick  10°  slew  will  be  much  different  than  those  for  a  60°  slew  due  to  the  fact  that  the 
control  system  operates  proportional  to  the  angle  error,  given  the  same  hardware.  As  an 
example,  Figure  7  shows  the  time  history  of  the  system  variables  for  the  underdamped 
case.  Both  10°  (solid  lines)  and  20°  (dashed  lines)  reorientations  occurred  using  the  same 
set  of  gains.  Notice  that  they  both  settle  at  the  same  time.  The  20°  slew,  however,  shows 
significantly  higher  acceleration  and  rate  commands.  The  phase  space  indicates  that  both 
slews  follow  essentially  the  same  trajectory,  and  the  magnitudes  are  scaled  since  the 
system  is  linear. 


Phase  Space 


10  15 

time  (sec) 


Figure  7.  Underdamped  Reorientation  (solid  10°,  dashed  20°  slews) 


It  would  be  expected  that  a  10°  slew  could  occur  about  twice  as  fast  as  a  20°  slew 
for  a  given  value  of  peak  acceleration,  but  naturally,  the  linear  control  system  prevents 
this.  As  a  further  example,  the  same  10°  slew  is  compared  with  a  60°  slew  to  show 
further  exaggeration  of  this  restriction.  Figure  8  demonstrates  that  the  required 
acceleration  and  rates  increase  significantly  with  the  slew  size. 
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Control  and  State  Variable  Trajectories 


Phase  Space 


o 

Figure  8. 


5  10  15  20  0  20  40  60  80 

time  (sec)  9(t) 

Underdamped  Reorientation  (solid  10°,  dashed  60°  slews) 


At  some  point,  the  actuators  will  not  be  able  to  produce  the  required  torque 
(acceleration)  and/or  angular  momentum  (rate).  Suppose  the  10/20°  slews  from  Figure  7 
are  now  subject  to  torque  and  momentum  saturation  limits  that  are  equivalent  to  5°/s2  and 
4°/s,  respectively  (values  representing  a  very  agile  spacecraft).  This  means  the  MED 
array  cannot  produce  more  than  these  limitations  around  the  desired  eigenaxis.  Figure  9 
demonstrates  the  saturations  by  flat  regions,  where  maximum  torque  and  acceleration  are 
held  for  extended  periods  of  time.  The  phase  space  indicates  both  trajectories  were 
saturated,  but  the  20°  slew  was  saturated  twice.  Within  the  same  20  seconds,  the  20°  slew 
was  not  complete  due  to  the  torque  and  momentum  limits. 
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Phase  Space 


Control  and  State  Variable  Trajectories 


0  5  10  15  20  0  10  20  30 

time  (sec)  9(t) 

Figure  9.  Underdamped  Reorientation  (solid  10°,  dashed  20°  slews)  a,  co  Limited 


Notice  that  the  spacecraft  began  at  maximum  torque  and  held  it  until  the 
momentum  limit  was  reached.  At  that  point,  maximum  momentum  was  held  until 
crossing  the  target  angle,  after  which  the  spacecraft  decelerated  at  maximum  torque. 

Engineers  typically  want  to  optimize  trajectories,  minimizing  overshoot  and 
settling  time  to  reduce  wasted  transient  distance  and  time.  To  perform  a  slew  quickly 
with  no  overshoot,  the  critically  damped  case  may  be  tried  (Figure  10). 


Control  and  State  Variable  Trajectories 


Phase  Space 
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Figure  10.  Critically  Damped  Reorientation  (solid  10°,  dashed  20°  slews) 
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Notice  that  in  Figure  10,  both  slews  do  not  require  anywhere  near  maximum  torque 
or  momentum.  In  this  simulation,  using  the  critically  damped  gains  seems  like  a  good 
approach.  In  a  real  spacecraft,  the  mass  (inertia)  will  decrease  over  time  as  fuel  is 
expended.  Thus,  space  system  engineers  must  plan  in  a  buffer  factor  in  order  to  account 
for  the  change  in  mass.  As  a  “rule  of  thumb,”  an  engineer  may  select  an  initial  damping 
ratio  around  £  =  0.7.  As  mass  decreases,  the  damping  ratio  increases  since  the  two 

values  are  inversely  proportional  (£  °=  1  /  m) .  The  damping  ratio  will  eventually  approach 

critical  damping  and  move  into  overdamped  towards  the  end  of  the  life  of  the  spacecraft. 
Thus,  one  could  argue  that  the  critical  damping  or  near-critical  damping  scenarios  would 
not  be  useful  in  spacecraft  design. 

A  simple  approach  to  control  system  design  is  to  devise  a  set  of  gains  that  work 
across  a  large  range  of  slew  angles.  To  satisfy  all  angles,  the  control  engineer  must  design 
the  gains  for  the  largest  slew  angle  the  spacecraft  is  regularly  expected  to  perform.  If  this 
is  done  using  a  linear  model,  the  system  may  perform  very  poorly  for  small  angle 
reorientations.  Therefore,  the  engineer  can  be  forced  to  overdesign  the  hardware  to  satisfy 
the  torque  and  momentum  requirements  for  smaller  angle  maneuvers. 

D.  TRAJECTORY  SHAPING 

Gain  design  is  very  delicate  because  the  engineer  must  tune  the  gains  special  for 
the  spacecraft,  regardless  of  slew  size.  Instead,  input  shaping  allows  a  fixed  set  of  gains 
that  must  correct  for  only  a  discrete  minimal  error  to  be  utilized  for  multiple  slew  sizes. 
Optimal  trajectories  can  be  pursued  without  major  hardware  upgrades. 

Input  shaping  is  an  alternative  to  using  a  step  command  as  the  input  to  the  attitude 
control  system.  Input  shaping  can  be  introduced  in  the  form  of  a  set  of  shaped  torque 
commands,  represented  as  an  acceleration  input  in  Figure  13  [17]— [18].  This  is  ideal, 
because  it  completely  avoids  the  need  for  a  controller.  Torque/acceleration  input  shaping 
is,  however,  unrealistic  because  the  system  would  be  operated  in  the  open  loop.  This 
means  that  it  does  not  have  a  feedback  mechanism  to  account  for  uncertainties  in  the 
spacecraft  inertia  and  external  disturbances  on  the  spacecraft,  which  is  undesirable. 
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Desired  Spacecraft 

Torque  1 

Trajectory 


F  rv  T 

1 

CO 
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6 

s 
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Figure  1 1 .  Torque  Input  Shaping  Model 


A  better  approach  is  to  shape  the  trajectory  of  the  input  to  a  closed-loop  system: 
the  input  angle  itself.  Rather  than  simply  inserting  a  desired  angle  0d  and  letting  the 

onboard  controller  calculate  how  to  respond.  The  desired  spacecraft  response  can  be 
decided  ahead  of  time  and  then  implemented  as  a  series  of  small  steps.  For  example,  a 
10°  slew  can  be  divided  into  100  or  1,000  incremental  steps.  Now,  the  gains  can  be 
optimized  for  the  step  size  (fixed  or  extremely  small)  rather  than  the  entire  maneuver  size 
(variable). 


Desired 

Input 

Trajectory 


Very 

Small  Controller  Spacecraft 

Input 

Error  I  1  I  l  1  l 


Figure  12.  Angle  Input  Shaping  Model 


As  long  as  the  spacecraft  model  is  reasonably  accurate,  the  path  can  be 
predetennined  accurately.  The  feedback  loop  then  adjusts  the  path  to  accommodate 
spacecraft  uncertainties  and  external  disturbances.  Similar  methods  have  been  used  in  at 
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least  a  few  cases  in  the  industry  [  12]— [  14].  The  input  shaped  trajectory  0d  is  fed  to  the 
system  much  like  that  in  Figure  12. 


The  ideal  input  is  the  trapezoidal  input  (referring  to  shape  of  CO  )  shown  in  Figure 
13.  This  input  accounts  for  both  the  torque  and  momentum  limits  a  priori  to  control  the 
error  0d  so  that  the  control  system  does  not  saturate. 


f 

time  angle  6d 

Figure  13.  Rate  Limited  Eigenaxis  Slew 

Due  to  the  shape  of  the  angular  rate  trajectory,  the  torque/momentum  limited 
input  is  also  referred  to  as  a  bang-coast-bang  maneuver.  This  is  exactly  the  command¬ 
generated  trajectory  used  for  the  Clementine  spacecraft  that  was  used  for  lunar  mapping 
in  1994  [12,  Figure  2].  The  eigenaxis  and  trajectory  angle  were  first  calculated  by  the 
onboard  computer  to  transit  from  the  current  quaternion  to  the  desired  quaternion.  In 
Chapter  V,  it  will  be  confirmed  that  this  trapezoidal  input  is  in  fact  the  time-optimal 
solution  for  eigenaxis  slews. 


28 


IV.  REACTION  WHEELS 


Reaction  wheels  have  their  physical  limitations.  They  have  maximum  angular 
rates  and  accelerations  that  must  be  modeled.  This  will  ensure  that  even  if  a  large  torque 
is  demanded,  only  the  maximum  MED  assembly  torque  will  be  applied.  In  order  to  apply 
the  saturations,  a  model  of  the  reaction  wheels  must  be  created. 

A  basic  model  of  a  reaction  wheel  is  shown  in  Figure  14.  The  reaction  wheel  is 
shown  in  its  own  local  frame,  referenced  as  the  wheel  frame  with  its  spin  axis  aligned 
with  the  m>3  axis.  A  motor,  attached  to  the  body  frame,  drives  rotation  of  the  wheel.  The 
principle  of  operation  is  for  the  motor  to  apply  a  torque  to  the  wheel.  This  causes  the 
wheel  to  rotate  quickly.  The  resultant  equal  and  opposite  torque  is  applied  back  through 
the  motor  to  the  body  frame.  This  causes  a  rotation  of  the  much  larger  spacecraft  body. 
The  momentum  and  torque  vectors  align  with  the  w3  axis. 


/V 


Figure  14.  Reaction  Wheel  Schematic 
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A. 


ORIENTATION 


The  angular  momentum  of  a  single  reaction  wheel  (RW)  within  its  (ith)  wheel 

frame  W  with  respect  to  body  can  be  shown  in  the  spacecraft  body  frame  B  through  a 
simple  rotation.  The  fundamental  formula  for  this  change  of  basis  is  given  as 

h‘=["Rw‘]h,w'  (43) 


In  dyadic  fonn,  (43)  can  be  looked  at  as  the  projection  of  the  wheel  frame  onto 
the  body  frame  (component- wise  dot  product). 


4*4i 

4*4-2 

4*4-3 

4*  4. 

4*4-2 

4*4-3 

h'v- 

(44) 

4  •  4i 

4*4-2 

4*4-3 

This  transfonnation  can  be  completed  as  a  matter  of  three  fundamental  rotations, 
and  represented  as  a  direction  cosine  matrix  (DCM).  A  DCM  is  an  invertible  square 
matrix  where  the  transpose  is  the  inverse  such  that 

bRw*  =  [w‘RbJ  or  V[Vj  =7  (45) 

This  property  allows  a  much  more  straightforward  approach  of  first  placing  the 
RW  in  the  B  frame,  converting  to  the  Wj  frame  via  W‘RB ,  and  then  using  the  transpose 
bRw‘  to  convert  back.  In  the  example  shown  in  Figure  15a,  the  MED  is  placed  directly  on 

/V  /V 

the  B  xyz  origin,  rotating  about  the  b3  axis  per  Figure  14.  This  produces  an  angular 
momentum  hBw  directly  in  the  b3  direction  for  a  reaction  wheel. 
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Figure  15.  Reaction  Wheel  Transformation 


First,  a  fundamental  3  rotation  about  the  shared  b3  /  w3  axis  must  be  completed  to 
orient  the  wx  and  w2  axes  for  subsequent  rotations  (Figure  15b).  This  will  allow  for 
multiple  reaction  wheels  placed  along  the  same  skew  angle  (/3),  forming  a  pyramid.  This 
first  rotation  can  be  represented  in  (46). 


*3  («,-)  = 


cos  a;  sina;  0 
-sina(.  cosa,.  0 
0  0  1 


(46) 


Assuming  a  =  constant  for  the  initial  3  rotation  of  the  representative  wheel,  a 
2  rotation  is  now  completed  about  the  w2  axis  (Figure  15c).  This  rotation  is  used  for  both 
reaction  wheels  and  CMGs,  and  is  referred  to  as  the  skew  angle.  The  skew  angle  is 
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normally  assumed  the  same  across  all  devices,  so  the  subscript  will  not  be  used  to 
distinguish  which  MED  it  refers  to.  This  rotation  can  be  represented  in  Equation  (47). 


KjO 3)  = 


f  cos P  0  -sin p  ^ 
0  1  0 

sin  / 3  0  cos  1 3 


(47) 


The  previous  series  of  rotations  can  now  be  combined  into  a  single  transformation 
converting  a  vector  from  the  body  frame  to  the  wheel  frame. 


w<RB  =  R2(P)R2(ai) 


1 

c> 
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1 

Co 

_ 1 

ca; 

sat 

0 

0 

1 

0 

-sat 

cai 

0 

1 

"CD 

0 

cp 

2 

0 

0 

1 

(48) 


catcP 

-SCCi 

COl'S  (3 


SCCtCfi  -sp 
cat  0 
sa, sfi  cP 


It  follows  from  Equation  (45)  that  BRW<  is  now  the  transpose,  which  is  equivalent 
to  the  reverse  sequence  via  the  negative  angle  of  each  rotation.  Equation  (49)  is  the 
physical  representation  of  the  dot  product  matrix  shown  in  Equation  (44). 


BRW‘  =R3(-ai )R2(-P)  =  [w'RbJ 


B 


I W: 


ca^P  -saj  casP 
sa^P  caj  sasP 


-sP  0  cP 


(49) 


Note  also  that 


b^wn 

bx*wi2 

V^3 

caxP 

-sai 

casp 

b2*wn 

b2*wi2 

4*^3 

= 

satcp 

caj 

sasp 

1 

u> 

• 

b2*wi2 

4  •  W',3 

-Sp 

0 

cP 

(50) 
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As  an  example,  place  the  initial  MED  at  angle  a  =  0° .  The  resultant 

transformation  from  the  W  frame  to  the  B  frame  is  shown  in  Equation  (51)  for  a 
reaction  wheel. 


B  DW 

A, 


B 


r3(-o°)r2(-P)= 


c/3  0  s/3 

0  1  0 
-s/3  0  c/3 


(51) 


Recall  that  the  RWs  in  Equation  (51)  were  stationed  at  a  =  0° .  A  classic 
configuration  is  using  four  MEDs  in  a  square  pyramid  configuration  in  which  each  of  the 
MEDs  are  placed  90°  apart  [15].  Arranged  at  the  proper  skew  angle,  this  tricetrix 
arrangement  maximizes  torques  in  all  directions,  minimizing  torque  cancellations.  This 
also  provides  redundancy  in  case  of  a  single  wheel  failure.  The  arrangement  used  will  be: 
a,  =  0 °,a2  =  90°, a3  =  180°, a4  =  270° .  Figure  16  demonstrates  this  common  placement. 
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Based  on  the  wheel  axis  description  mentioned  above,  the  angular  momentum  in 
the  wheel  frames  can  be  written  as 


h 


wi 

rWj 


0 

0 

1 


I  w. 


(52) 


Transforming  angular  momentum  from  the  W  frame  to  the  B  frame  using 
Equation  (43)  and  a  =  0°  yields  (53). 


or  in  general  form 


I B 


= 


bi'\ 
b2'\ 
b,  •  w. 


(53) 


(54) 


The  total  angular  momentum  of  the  system  is  simply  the  sum  of  each  of  the 
MEDs  in  the  body  frame. 


h 


B 

medA 


(55) 


For  the  reaction  wheel  case  in  Equation  (54),  the  individual  rotation  matrix  and 
angular  momentum  components  can  be  combined  into  the  product  of  two  matrices, 
equivalent  to  Equation  (55).  The  first  matrix  represents  the  orientation  of  each  wheel  with 
respect  to  the  body  frame,  separated  as  individual  columns.  For  n  MEDs,  this  3 xn 
orientation  matrix  is  multiplied  by  an  nx  1  vector  containing  the  momentum  in  the  wheel 
frame  as  shown  in  Equation  (56). 
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(56) 


The  3 xn  matrix  that  relates  n  MEDs  to  the  principle  body  axes  will  be  labeled  this 
point  forward  as  the  orientation  matrix  BZW  .  Since  it  contains  only  elements  from  the  3rd 
column  of  the  original  rotation  matrix,  it  will  be  further  described  with  a  subscript  3: 

BrjW 

z,3  . 


BZ 


w 

3 


B 


brwh 

brwh 

brw, 

b2'Wh 

b2  •  w2} 

••  b2-wj 

b3'K 

63 

Now,  Equation  (56)  can  be  rewritten  as 


C = = 'z,"'*"' 


where  hw  = 


K 

h2 


(57) 


(58) 


(59) 


For  the  square  pyramid  mentioned,  the  orientation  matrix  is  shown  in  Equation 
(60)  for  reaction  wheels. 
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(60) 
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Any  RW  orientation  configuration  Z  can  be  modeled  by  specifying 
the  appropriate  a  and  /?  rotations.  Note  that  Equation  (60)  assumes 

a  =  [  0  90  180  270 

Other  than  the  case  where  n=3,  Z  will  always  be  a  3xn  matrix  where  n>3.  This  is 
caused  because  the  system  is  over-actuated,  i.e.,  more  wheels  exist  than  are  necessary  to 
control  the  spacecraft  attitude.  A  common  approach  for  performing  control  allocation  in 
the  over-actuated  case  is  to  use  the  pseudo-inverse  to  map  three  orthogonal  commands  in 
the  body  frame  to  n  commands  for  the  individual  wheels. 

Z+=ZT(zZTf'  (61) 


B.  DIFFERENTIAL  EQUATIONS  OF  MOTION 

Like  the  spacecraft,  the  momentum  of  the  reaction  wheel  can  be  defined  as  the 
product  of  wheel  inertia  and  angular  velocity  of  the  wheel,  all  with  respect  to  the  inertial 
frame. 


ZJ  —  T 

H  h  —  Jwi  (Q  j 


(62) 


Previously,  the  wheel  inertia  was  combined  into  the  spacecraft  inertia  in  Equation 

(6),  where  h^,c  -  {js/chody  +JmedA ) N (On  .  Therefore,  the  wheel  inertia  rotating  with  the 

spacecraft  body  must  be  pulled  out,  leaving  only  the  configuration  of  the  wheel  with 
respect  to  the  body. 


Hw=Jj<o,B  +  JwiB(Dw‘ 


(63) 
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(64) 


Shown  in  dyadic  fonn: 
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(65) 


K,W'  =  (jWil  1 W  +  jyrinW  +  jyAilW  )  *  ) 

=  jy*3fli\  (w/3*  Wfj) 

=  j,nPAh 


Change  of  basis  to  the  body  frame: 


hJ  =  •  (kk +kk+kk) 

=  Li,p,  (w<3  ’k)k  +  jwhPi  (wh  •  b2  )b2  +  ( w,.3  •  b.)k 


(66) 


The  dot  product  of  two  vector  components  can  be  flipped  interchangeably  (i.e., 
hj  •  wa  =  w(l  *6, ).  This  is  important  to  note  because  the  transfonnation  matrix  can  now  be 
written  in  the  more  conventional  notation  shown  in  Equation  (44).  The  resultant  rotation 
matrix  is  shown  in  Equation  (67). 
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(67) 


Equation  (67)  can  be  simplified  into  a  version  containing  only  the  rotation  matrix 
portion  multiplied  by  its  momentum,  expressed  in  the  body  frame. 


h  r  = 


brA3 

b2>wh 

km\ 


i  .  Q. 

J  wi33  i 


(68) 


From  Equation  (56),  the  total  angular  momentum  for  n  wheels  in  the  body  frame 
can  be  rewritten  as 
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hB  = 


b2>w. 


b^w 


brw. 


b2  •  w2i 

k  •  w. 


jw  Q, 
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(69) 


Equation  (69)  is  equivalent  to  Equation  (58),  where  hB  =  bZ"hwv  .  Within  the 
model,  the  inertia  and  rates  are  expanded  out  in  the  following  manner  to  keep  track  of 
them  separately. 
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(70) 


The  torque  with  respect  to  the  wheel  frame  is  the  time  derivative  of  Equation  (70). 


T1*  =  ■ 


dt 


(K)=J,.nw 


(71) 


The  torque  with  respect  to  the  body  will  be  defined  as  the  change  in  angular 
momentum  with  between  the  wheel  and  body  frames. 


<  =  ~£(h*)  =  kl  where  hB  =  BZ"h 


W 


(72) 


Expanding  out  these  terms  following  the  chain  rule  yields  Equation  (73).  Since 
this  is  a  reaction  wheel,  the  orientation  is  fixed  and  therefore  RZf  =  0 .  The  inertia  is 
assumed  to  not  change  with  time,  Jrw  =  0  . 


b^w  r  ntt' 


JM 


(73) 

(74) 


In  dyadic  form,  a  single  reaction  wheel  torque  can  be  found  using  the  derivative 
of  Equation  (66). 
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-^(A/ )  =  UPi  (wh  *4)4  + ihPt  (wh  •  4  )4 + ]m?p,  (wh  •  4)4 

+  PA (wh  *4)4  +  PA  (w/3  •  4 )b2  +  jwA  (^3  •  4)4  (7S) 

+ PA  (™h  •  4  )4 + PA  (A  •  4  )4 + iw,,,0/  (^3  •  4)4 

The  inertia  of  the  wheel  is  not  changing  with  respect  to  the  wheel  frame.  Also,  the 
orientation  of  the  3  axis  of  wheel  frame  is  not  moving  with  respect  to  the  body  frame. 


.4*„  =  0  and  w,  =  0 


Equation  (75)  now  becomes 


) = p3 Pi  (^3  *4)4+ pA  (v4  *4)4+ pA  (43  •  4  )4 
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J  xvi^ 3  (  1  3  / 


V*),  PA 


V*), 


By  Euler’s  transport  theorem  adopted  from  [23,  Algorithm  8.18a],  the  inertial 
derivative  of  Equation  (74)  is 


•d{h:)_Bd{K) 


+  No)Bxh“ 


/,)'  =  ht  +  N(0B  x  ht 


Equation  (79)  can  be  rewritten  by  substituting  the  body  torque  from  Equation  (74) 
and  body  angular  momentum  terms  as 


-N  BryW  J  rtW  .  N-.B  ..  BryW  T 

th  =  Z3  J  Q  +  (0  x  Z3  J  J2 


By  substituting  Equation  (80)  into  Equation  (12),  knowledge  that  torque  is  in  the 
wheel  frame  as  defined  in  p  =  J  O  from  Equation  (71),  and  a  rearrangement  of  terms, 
the  result  is  an  equation  in  terms  of  the  torque  of  the  reaction  wheel.  This  equation  now 
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involves  the  pseudo-inverse  to  translate  torque  commands  given  in  the  body  frame  to  n 
wheels,  shown  in  Equation  (82). 

-N  BryW  -W  .  N  ,-B  ..  BryW  J  ^  W  /Q  U 

V,+  ®  X  Z3  J  J2  (81) 

V  BrjW~\+ 1  _ N  N,^B.,BrjWj  f\W\ 

T»'  =  L  Z  J  [~Treq  ~  (0  X  Z,  J  J2  j  (82) 

Further,  Equation  (82)  can  be  reduced  to  Equation  (83),  solving  for  the  reaction 
wheel  dynamics  in  terms  of  the  spacecraft. 

£2  =JW  l  Z3  J  [-Treq-  0)  X  Z  J  Q  )  (83) 

C.  MODELING  REACTION  WHEELS 

In  order  to  apply  the  physical  limits  of  the  reaction  wheel  actuator,  the  torque 
required  of  the  assembly  must  be  broken  down,  torque  and  wheel  speed  limit  saturations 
applied,  and  then  the  torque  of  the  reaction  wheel  assembly  rebuilt  in  Figure  17.  Equation 
(82)  is  used  to  break  down  the  torque  required  and  (80)  is  used  to  rebuild  it. 
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Figure  17.  Reaction  Wheel  Model 


This  model  is  contingent  that  momentum  in  the  wheel  frame  is  derived  directly 
from  wheel  torque  from  Equation  (71)  such  that 
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hw  =  f Twdt  . 

w  J  w 


(84) 


The  wheel  torque  limitations  occur  in  the  individual  wheel  frames.  The  rate 
limitations  are  modeled  as  wheel  momentum  saturations.  Note,  this  model  does  not  shut 
off  the  wheel  torque  when  |/zi  1 1  =  /imax .  It  is  not  necessary  for  the  comparisons  in  this 
thesis  because  only  torque  and  momentum  limited  signals  are  used. 


D.  REACTION  WHEEL  MOMENTUM  SPACE 

By  measuring  the  maximum  momentum  in  any  given  direction  over  a  2 K 
steradian  for  a  particular  reaction  wheel  configuration,  a  boundary  can  be  defined  that 
describes  the  maximum  momentum  capacity  of  the  array.  This  is  the  physical  limitation 
of  the  system  defined  by  the  vector  sum  of  the  maximum  capacities  of  the  individual 
momentum  wheels.  If  this  boundary  is  created  for  every  possible  direction  around  a  unit 
circle  in  the  body  frame,  the  surface  of  boundaries  is  typically  known  as  the  momentum 
envelope  [15,  p.  1607].  To  determine  the  best  possible  mechanical  configuration  of 
systems,  the  volume  must  be  maximized,  ideally  forming  the  most  spherical  shape 
possible.  Since  each  of  four  wheels  is  maximized  individually,  the  resultant  shape  has 
radii  larger  than  hmax .  If  all  four  wheels  were  aligned,  a  maximum  radius  of  4  could  be 
achieved,  leaving  zero  momentum  in  other  directions.  The  maximum  momentum  space 
for  an  evenly-spaced  (Aa  =  90°)  four  wheel  system  is  achieved  when 


[ 3  =  cos  1 


f  _n 

vV3. 


rad  =  54.74°  . 


(85) 


The  maximum  angle  in  Equation  (85)  is  shown  in  [15,  p.  1610]  as  the  conjugate 
35.26°.  The  resulting  envelope  is  shown  in  Figure  18.  The  figure  was  created  using  a 
mesh  of  triangles  between  all  surface  points. 

Assuming  all  four  wheels  have  equal  momentum  storage,  Figure  17  is  calculated 
as  a  function  of  the  maximum  momentum  available  to  any  given  wheel  hmia  .  The 
envelope  scales  symmetrically  with  the  size  of  the  wheel  ( /?max  ).  Therefore,  larger  wheels 
would  push  the  surface  outward,  allowing  more  momentum  capacity,  but  the  shape  would 
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be  retained.  This  is  why  the  dodecahedron  is  generic  for  any  four-wheel  system,  but 
unique  to  only  four-wheel  systems.  A  five-wheel  system  would  have  a  much  different 
shape  with  a  maximum  at  a  different  angle  [5  and  a  larger  volume  (see  [15]). 
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Figure  18.  Momentum  Space  for  4  Wheel  Reaction  Wheel  (a=  0,90,180,270°) 


The  momentum  envelope  shown  represents  the  angular  momentum  available 
about  any  given  rotation  axis.  During  times  of  maximum  momentum  use,  the  momentum 
vector  would  touch  the  outside  surface  of  the  envelope  until  the  system  begins  to 
decelerate  and  momentum  decreases  back  to  the  center  of  the  shape  where  momentum  is 
zero.  The  vertices  of  the  momentum  envelope  indicate  optimal  axes  that  allow  locally- 
rapid  rotation  rates.  Although  fourteen  vertices  exist,  only  six  are  global  maxima  (marked 
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with  blue).  The  centers  of  each  facet  indicate  the  minima  for  each  respective  plane. 
These  minima  are  also  the  global  minima,  each  marked  in  magenta. 

The  four-wheel  polyhedron  is  again  shown  in  three  orthogonal  views  in  Figure  19 
to  clarify  its  shape.  The  blue  maximum  vertices  (global  maxima)  occur  at  magnitudes 

equal  to  /tfax  =  2.309/zmax .  The  black  minor  vertices  (local  maxima)  occur  at 

^max  =  2-0/zmax .  The  minimum  hB  is  actually  zero,  but  !\]m  represents  the  minimum  across 
the  momentum  envelope.  The  minimum  radius  occurs  at  the  middle  of  the  flat  facets 
where  h^in  =  1 .633 /zmax ,  interestingly  equal  for  every  axis. 


Max  Vertices  1*1  =  2.309/?^ 
Minor  Vertices  |/i|  =  2.000/?max 
Minimum  1/(1  =  1 .633 


-  "i 

Figure  19.  Momentum  Space  of  a  4  Reaction  Wheel  System 
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Table  1  lists  all  of  the  possible  configurations  of  momentum  vertices  and  facet 
minimums.  Once  normalized,  these  represent  the  respective  axes  with  the  momentum 
magnitudes  listed  above. 


Max  Vertices  Minor  Vertices  Minimum  Facets 

[  1  1  0  ]  [  sin(P)  0  cos(P)  ]  [  1  0  0  ] 

[  1-1  0  ]  [  sin(P)  0  -cos(P)  ]  [  -1  0  0  ] 

[  -1  1  0  ]  [  -sin(P)  0  cos(P)  ]  [  0  1  0  ] 

[  -1  -1  0  ]  [  -sin(P)  0  -cos(P)  ]  [  0  -1  0  ] 

[  0  0  1  ]  [  0  sin(P)  cos(P)  ]  [  sin(P)  sin(P)  2cos(P)  ] 

[  0  0  -1  ]  [  0  sin(P)  -cos(P)  ]  [  sin(P)  sin(P)  -2cos(P)  ] 

[  0  -sin(P)  cos(P)  ]  [  sin(P)  -sin(P)  2cos(P)  ] 

[  0  -sin(P)  -cos(P)  ]  [  sin(P)  -sin(P)  -2cos(P)  ] 

[  -sin(P)  sin(P)  2cos(P)  ] 

[  -sin(P)  sin(P)  -2cos(P)  ] 

[  -sin(P)  -sin(P)  2cos(P)  ] 

[  -sin(P)  -sin(P)  -2cos(P)  ] 

Table  1.  Reaction  Wheel  Significant  Axes 

From  Figure  17,  the  wheel  torque  vectors  directly  align  with  the  wheel 
momentum  vectors.  This  causes  a  torque  envelope  equal  in  shape  to  that  in  Figures  18 
and  19.  Instead  of  using  hmm  as  the  scaling  factor,  the  dodecahedron  would  be  instead 
scaled  by  the  maximum  individual  wheel  torque  Trnax .  All  maximum  and  minimum 
vertices  are  co-aligned  with  the  momentum  vertices,  just  with  different  units. 

E.  LARGEST  INSCRIBED  SPHERE 

As  mentioned  in  Chapter  I,  sometimes  it  is  necessary  to  create  artificial 
limitations  on  a  control  system.  Recall  from  Figure  18  that  the  twelve  facet  minima  occur 
at  /? f  n  =  1.633hmax .  Although  they  are  considered  minima,  they  still  lie  on  the  momentum 

surface.  For  example,  on  axis[  1  0  0  ],  momentum  up  to  1.633 hmax  (torque  up  to 
1  -633  Tmax )  is  available.  This  is  the  physical  system  limit  that  cannot  be  exceeded  for  an 
eigenaxis  rotation  about  that  specific  axis.  Since  this  is  the  global  envelope  minimum,  a 
sphere  can  be  created  with  radius  1.633 hmax  that  includes  all  twelve  points.  Such  a  sphere 
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is  the  largest  possible  sphere  within  the  momentum  space,  shown  in  Figure  20. 
Restricting  commanded  momentum  to  lie  within  the  inscribed  spherical  volume  ensures 
that  1 ,633/zmax  is  available  about  every  possible  axis. 


Figure  20.  Reaction  Wheel  Momentum  Space  with  Largest  Inscribed  Sphere 


Of  course,  40%  more  momentum  is  available  about  the  maximum  axes  shown 
in  Table  1,  but  in  a  conventional  attitude  control  system  design,  this  additional 
torque/momentum  capacity  is  sacrificed.  This  sacrifice  is  made  to  ensure  a  designed  slew 
rate  is  available  for  any  given  slew.  The  inertia  of  a  spacecraft  is  rarely  symmetric  and 
equal  about  every  major  axis.  The  resultant  is  an  inertia  ellipsoid,  further  reducing  the 
momentum  space  envelope.  This  occurs  due  to  the  need  to  ensure  the  momentum  does 
not  change  with  a  non-symmetric  inertia  matrix.  The  maximum  spacecraft  inertia  value  is 
also  used  about  every  axis,  further  limiting  the  system  to  enforce  consistent  slew  rates. 
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F. 


REACTION  WHEEL  PSEUDO-INVERSE  STEERING  LOGIC 


The  physical  limitation  of  a  reaction  wheel  system  was  previously  discussed. 
Another  limitation  associated  with  reaction  wheels  is  the  mathematical  limitation  of  the 
Moore-Penrose  pseudo-inverse,  Equation  (61).  When  reorienting  coordinates  from  the  n 
wheels  in  the  wheel  frame  W  to  three  orthogonal  coordinates  of  the  body  frame  B  ,  there 
is  a  unique  solution. 


hB  =  BZ™hw 
tb  =  bZ™tw 


(86) 


The  pseudo-inverse  offers  a  single  solution  for  translating  the  actuation  from  the 
body  coordinates  ( 3x1  vector)  unifonnly  across  n  wheels  (nxl  vector).  The  pseudo¬ 
inverse  is  desirable  because  it  is  the  minimum  least  squares  (L2)  solution,  minimizing  the 
sum  of  the  squares  of  all  of  the  wheels’  individual  momenta/torque  [15,  p.  1609]. 


hw=[BZ?JhB 
t  w=[BZf]+rB 


(87) 


Since  the  system  is  underdetermined,  there  are  infinite  possible  nxl  wheel 
configurations  that  will  satisfy  the  body  coordinate  solution,  but  the  pseudo-inverse  will 
give  only  one  such  configuration. 


Pseudo-Inverse  Control  Allocation 


Using  the  reaction  wheel  orientation  matrix  for  a  4  wheel  system 


(60)  and  /3  =  cos  1 


f  1  ^ 


vV3y 

inverse  used  in  this  example  is 


rad  =  54.74° ,  the  resultant  orientation  matrix 


mentioned  in 
and  pseudo- 


BrjW 


zw  = 


o 

l 


o 

V2 

1 


-V2  0 

0  -y/2 

1  1 


(88) 
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V2  0  1 

V3  0  V2  1 
4  -V2  0  1 

0  -V2  1 


(89) 


To  illustrate  the  implications  of  [ "Zl1  ]+,  consider  the  following  possible  wheel 
momentum  configuration  where  1 .0  hmax  represents  a  scalar  of  the  maximum  possible 
momentum  for  a  given  wheel. 

1 

h possible  ~  ^rnax  (90) 

1 

Using  Equation  (90)  leads  to  the  following  possible  momentum  in  the  body 

frame. 

1.633 
0 

1.155 

Assume  this  possible  hB  was  commanded  rather  than  derived.  Now,  the  pseudo¬ 
inverse  would  be  required  to  derive  the  individual  wheel  momentum  commands,  h "  . 

Applying  the  least-squares  solution  [  /fZl'  ]+  leads  to  a  different  solution  than 

1  i  -i  i  ]r).  The  pseudo-inverse  solution  is  actually  physically  unachievable 
because  the  largest  momentum  command,  1 .5  hrma  ,  is  higher  than  the  maximum  hmax  for  a 
given  wheel: 
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h 


W 

pseudo-unacheivable 


B 

possible 


1.633 

0 

1.155 


1.5 

0.5 

-0.5 

0.5 


|  W 


h 


max 


(92) 


In  order  to  implement  the  least  squares  control  allocation,  the  proportionality  of 
the  solution  must  remain  intact  [24,  pp.  880-883].  The  wheel  momentum  vector  in 
Equation  (92)  is  therefore  nonnalized  by  the  highest  element  (greatest  in  magnitude). 
This  allows  appropriate  control  allocation  to  rotate  the  spacecraft  predictably  about  the 
desired  axis. 


h 


w 

pseudo-  lim 


f 


max 


r  ByW  ~|+  iB 

-3  J  ndesired 


h 


W 

pseudo-unacheivable 


max 


d" 


pseudo-unacheivable 


) 


1.000 

0.333 

-0.333 

0.333 


I  w 


h 


max 


(93) 


The  normalized  solution  Equation  (93)  is  far  different  than  the  wheel  momentum 
configuration  at  the  beginning  of  this  example,  Equation  (90).  Translating  back  to  the 
body  frame  reveals  that  the  pseudo-inverse  limited  solution  in  Equation  (94)  can  only 
achieve  67%  of  the  possible  momentum  shown  in  Equation  (91). 


h 


B 

pseudo- lim 


ByW 

^  n  pseudo-  lim 


BrjW 


1.000 

0.333 

-0.333 

0.333 


I  w 


h 


max 


1.089 

0 

0.770 


I B 


h 


(94) 
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The  previous  example  can  be  redone  using 


(95) 


with  a  as  anything  from  0.333  to  1.0.  The  system  returns  the  same  limit  at 

1.089  T 

Kseudo-  Jim  =  0  /rmax  (96) 

0.770 

which  is  only  67%  of  the  possible  body  momentum  in  this  example. 

r  ~\B 

1.633 

k possible  ~  0  ^max  (97) 

1.155 

The  same  performance  limitation  results  in  a  loss  of  torque  per  Equation  (87). 


2.  Performance  Loss 

The  momentum  vector  in  the  previous  example  corresponds  to  a  minor  vertex  (see 
Figure  19). 

B  f  -i* 

1.633  c/3 

Kossibie=  0  hm  ax  =  o  2 .0  /zmax  (98) 

1.155  J  sp 

The  pseudo-inverse  limited  vector  is  67%  of  hBpossjble  giving  a  magnitude  of 
1*1  =  1-33  A_. 

B  f  ~\B 

1.089  cp 

Keudo-Km=  0  /jmax  =  0  1 -33/imax  =  0.67  (99) 

0.770  J  sp 
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Graphically,  h^seudo  Vim  would  be  well  within  the  reaction  wheel  momentum 

envelope.  Shown  in  Figure  21,  the  red  shape  represents  the  full  reaction  wheel 
momentum  space.  It  is  the  physical  limitation  of  the  reaction  wheel  system  utilizing  the 
maximum  momentum  possible  to  rotate  about  any  given  axis.  The  black  dot  represents 
the  original  point  (minor  vertex)  shown  in  Equation  (98).  The  red  dot  represents  the 
pseudo-inverse  limited  system  shown  in  Equation  (99).  The  blue/black  lines  represent  all 
possible  magnitudes  that  share  the  same  eigenaxis  as  the  minor  vertex.  The  blue  line 
specifically  represents  the  set  of  all  achievable  configurations  by  the  pseudo-inverse.  The 
black  line  represents  the  set  of  all  unachievable,  yet  still  possible,  configurations. 


Minor  Vertex 
Possible  Magnitudes 
Pseudo-Inverse  Limited  Max 

Possible  Pseudo-Inverse 
Limited  Magnitudes 


b\ 


Figure  2 1 .  Reaction  Wheel  Momentum/Torque  Space  Showing  Pseudo-Inverse  Limit 


The  same  is  true  for  all  other  minor  vertices,  as  they  are  unachievable  using  the 
pseudo-inverse.  On  the  other  hand,  the  minimum  points  and  the  maximum  vertices  are  all 
achievable,  shown  in  Figures  22  to  23. 
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Figure  22.  Reaction  Wheel  Pseudo-Inverse  Limited  Space 


The  representation  shown  in  Figure  22  makes  it  clear  that  the  pseudo-inverse  is  a 
mathematical  limitation  on  the  system  performance  and  not  a  physical  system  limitation. 
Indeed,  this  is  true  because  the  pseudo-inverse  minimizes  the  sum  of  the  squares  of  the 
moment/torque  in  a  given  direction.  The  least  squares  solution  discards  some  of  the 
available  moment  capacity  in  favor  of  an  L2  control  allocation.  In  reality,  the  spacecraft 
has  the  ability  to  rotate  faster  (higher  momentum)  about  many  axes,  but  the  pseudo¬ 
inverse  does  not  necessarily  allow  this.  Going  from  wheel  space  to  body  space,  there  is 
only  one  mathematical  answer.  Going  from  body  space  to  wheel  space,  there  are  an 
infinite  number  of  potential  arrangements.  The  pseudo-inverse  only  provides  the  least- 
squares  solution,  a  mathematically  limited  solution. 

Figure  23  shows  the  pseudo-inverse  limited  momentum  space  from  Figure  22 
overlaid  into  the  full  reaction  wheel  space  from  Figure  19.  This  makes  it  clear  how  much 
volume  is  missing  across  the  different  orthogonal  views.  An  interesting  note  is  that  the 
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pseudo-inverse  solution  intersects  all  six  of  the  maximum  vertices  and  all  twelve  of  the 
minimum  points,  but  none  of  the  eight  minor  vertices. 


Figure  23.  Reaction  Wheel  Pseudo-Inverse  Limited  Space  -  Orthogonal 


G.  PSEUDO-INVERSE  LARGEST  INSCRIBED  SPHERE 

As  previously  shown  with  the  full  reaction  wheel  space,  the  full  momentum 
envelope  can  be  further  limited  to  an  enclosed  sphere  (Figure  19).  Similarly,  the  pseudo¬ 
inverse  limited  space  can  be  confined  to  a  sphere  to  ensure  a  constant  value  of  achievable 
momentum  in  all  directions.  The  minimum  pseudo-inverse  envelope  surfaces  occur  at 
1 .333hmax .  Therefore,  the  pseudo-inverse  space  can  be  confined  to  a  sphere  with  a  radius 
of  1.333 hmm  (Figure  24).  This  will  ensure  that  this  momentum  is  available  in  all 
directions. 
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Figure  24.  Pseudo-Inverse  Momentum  Space  with  Largest  Sphere 

Figure  25  shows  a  comparison  of  the  reaction  wheel  momentum  space  under  the 
various  control  limitations.  Here  it  is  obvious  the  huge  loss  of  momentum  space  that 
occurs  under  the  typical  spherical  simplification  utilized  in  the  design  of  a  spacecraft 
attitude  control  system.  Later  chapters  will  be  used  to  discuss  methods,  particularly  the 
use  of  optimal  control,  to  improve  the  usage  of  the  available  system  capacity. 
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Figure  25.  Momentum  Space  Comparison 
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V.  OPTIMAL  REACTION  WHEEL  ATTITUDE  CONTROL 


This  chapter  will  be  used  to  derive  the  time-optimal  slew  for  a  reaction  wheel 
system.  First,  the  time-optimal  eigenaxis  slew  will  be  generated.  This  standard  heritage 
maneuver  is  best  only  if  the  engineer  wants  the  spacecraft  to  stay  on  the  eigenaxis  to 
minimize  the  overall  rotation  angle.  The  spacecraft  often  has  more  momentum  and  torque 
available  to  exploit,  but  it  is  necessary  to  deviate  from  an  eigenaxis  path  in  order  to 
accomplish  this.  The  remaining  sections  of  the  chapter  will  expand  the  linear  derivation 
from  Section  A  into  the  nonlinear  reaction  wheel  and  spacecraft  dynamics  developed  in 
Chapters  II  and  IV.  DIDO,  an  optimal  control  solver  [21],  will  be  used  to  obtain  the 
optimal  control  solution.  An  attitude  control  architecture  will  also  be  presented  in  order  to 
implement  the  optimal  solution. 

A.  TIME-OPTIMAL  EIGENAXIS  REORIENTATION 

The  time-optimal,  rest-to-rest  eigenaxis  slew  solution  is  fairly  straightforward. 
Not  all  of  the  following  steps  are  necessary  to  solve  this  problem,  but  all  steps  will  be 
covered  for  completeness.  The  optimal  control  methodology  is  necessary  to  understand 
prior  to  Section  B,  where  a  much  more  complicated  derivation  is  presented. 


1.  Problem  Formulation 


In  general,  figures  of  merit  are  used  to  determine  how  well  the  slew  meets  the 
requirements  important  to  the  engineers  and/or  customer.  These  are  broken  up  into  costs 
(stuff  to  minimize)  and  constraints  (boundary  conditions  that  must  not  be  exceeded).  The 
maneuver  is  defined  as  rest-to-rest  with  the  initial  and  final  angles  predetermined. 

Initial  time  is  zero  and  the  final  time  is  free  (tenn  to  be  optimized),  if  it  is  desired 
to  design  a  rapid  slew  maneuver. 


t0=t=  0 


tf  =  r 


(100) 
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As  shown  in  Chapter  III,  the  spacecraft  model  can  be  simplified  as  a  double 
integrator.  In  this  case,  the  state  variables  and  control  are 
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u  -  [m(0]  =  [ct(0] 

The  state  space  representation  of  Equation  (101)  is  shown  as 
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x  =  f(x(t),u(t)) 
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(103) 


First  and  foremost,  the  solution  must  meet  the  boundary  conditions.  The  purpose 
of  slewing  is  to  accurately  point  at  a  new  target  (end  condition).  If  the  spacecraft  slews 
but  does  not  get  the  desired  end  condition,  then  it  may  not  be  possible  to  perfonn  the 
mission  requirements.  The  final  angle  will  be  defined  as  the  difference  from  start  to 
finish.  Therefore,  the  initial  angle  will  be  set  as  zero.  Initial  and  final  rates  are  zero 
because  the  maneuver  is  considered  rest-to-rest.  The  boundary  conditions  are  defined  as 
follows. 


0(C)  =  #0  =  ° 
co(t0)  =  co0  =  0 

e(tf)=df  =  ed 

co(tf)  =  cof  =0 


These  boundary  conditions  can  be  rewritten  as 


*(f0)  =  *°  = 


x(tf)  =  xf 
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A  convention  to  be  used  to  define  the  endpoint  conditions  e  will  be  to  set 
Equation  (106)  equal  to  zero. 
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(*('/)) 


Of  od 

0 

(Of 

0 

(107) 


The  Mayer  cost  E  is  a  fonn  of  endpoint  cost.  Assuming  the  boundary  conditions 
“must”  be  met,  this  is  the  measure  of  optimality  that  quantifies  the  solution  in  terms  of 
the  final  state.  Since  this  is  a  time-optimal  slew,  the  endpoint  cost  is  the  final  time. 

E{x(tf)}  =  tf  (108) 


The  Lagrange  cost  F  is  a  fonn  of  running  cost.  This  measure  quantifies  the  cost  of 
the  maneuver  itself.  A  common  form  is  the  quadratic  cost  F  =  f fu2(t)dt ,  which 

J'o 

represents  the  “energy”  or  efficiency  of  the  system.  Since  the  control  u  is  directly  related 
to  torque,  this  cost  can  be  used  to  minimize  torque-related  functions  like  power  during  a 
maneuver.  For  the  simplified  case  of  individual  time-optimal  slews,  the  Lagrange  cost 
will  not  be  used.  The  superscript  (D.N.E.)  will  be  used  to  show  that  the  term  is  not  only 
zero,  but  it  does  not  exist  for  this  problem. 

I" 1  F(x(t),u(t))dt  =  0(i"Vi')  (109) 

^  ro 

The  cost  functional  3  can  be  composed  of  both  Mayer  and  Lagrange  costs. 


5[x(-),M(-),t/]  =  £’(x(t/))  +  J  F(x(t),u(t))dt 
3\x(-),u(-),tf]  =  tf 


(110) 


Rewriting  this  problem  in  the  format  used  by  Ross  [25,  p.  43],  we  get 
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Given 

Minimize 


M 


S[x(-),K(-)>f/]  =  f/ 


Subject  to 


{e0,co0,t0)=(o,o,o) 
(df-ed,cof)={  o,o) 


(in) 


2.  The  Hamiltonian 

For  every  state,  a  costate  will  be  defined  with  subscripts  to  match  their  respective 

states. 


K 

A, 


(112) 


The  control  Hamiltonian  is  given  by 

H[X,x,u )  =  F(jc,«)+  A,1  f(x,u) 


q(d.n.e.)  + 

(0 

L  -1 

U 

(113) 


(114) 


=  \co  +  Aojm 


3.  The  Hamiltonian  Minimization  Condition 

The  Hamiltonian  Minimization  Condition  states  that  the  Hamiltonian  H  must  be 
minimized  over  the  control  u  [26,  p.  108].  This  leads  to  the  Euler  Lagrange  equation. 


d_ 

du 


(?ieC0  +  \0u)  =  K  =  0 


(115) 


Since  u  does  not  appear,  Pontryagin’s  Principle  must  be  used.  From  [26, 
Algorithm  3.8.13],  Pontryagin’s  Principle  states  the  following. 
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<)H 

du 


>0, 
=  0, 
<0, 


— u  <u  <u 

max  max 


U  =  Un 


(116) 


Applied  to  this  problem,  the  switching  function  in  Equation  (1 16)  becomes 


K  =  < 


>0, 
=  0, 
<0, 


-u  <  u  <  u 

max  max 


u-un 


(117) 


By  rearranging  Equation  (117)  in  terms  of  the  costate  to  solve  for  the  control,  the 
rule  in  Equation  (118)  results.  Since  the  value  of  the  control  does  not  matter  when  the 
costate  -  0 ,  the  control  can  be  set  to  zero.  This  is  not  necessary,  but  it  is  acceptable 

since  the  use  of  any  control  during  this  state  will  not  affect  the  optimal  control  problem. 
This  is  stated  simply  as  an  energy  conservation  technique.  Thus,  the  control  is  given  as: 


"“max  K  >  0 

0  for  K  =  0 


“max  K  <  0 


(118) 


Ultimately,  Equation  (118)  states  that  the  control  must  be  at  its  positive  or 
negative  maxima.  The  function  sgn  will  be  defined  as  the  sign  of  its  independent  variable. 


sgn(.): 


1  (•)>« 

0  for  (•)  =  0 

-1  (*)<o 


(119) 


Equation  (118)  can  now  be  written  compactly  as 

u  =  -sgn(A(0)«max 


(120) 


Equation  (120)  matches  the  formulation  by  [25,  Algorithm  3.36]. 


4.  The  Costate  Dynamics  and  the  Adjoint  Equation 

From  [25,  Algorithm  2.2],  the  adjoint  equation  states 
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(121) 


,■ ,  ,  dH 

-A(0  =  3- 

dx 


(122) 


The  differential  costate  equations  shown  in  Equation  (122)  can  be  integrated  to 
create  equations  for  the  costates  in  terms  of  unknown  constants. 


— Cjf  +  c2 


(123) 


Transversality 


The  endpoint  Lagrangian  is  given  by  [25,  Algorithm  2.9]. 


E(xf,tf)  =  E(x(tfj)  +  vTe(x(tfj) 


(124) 


E(xf,tf)  =  tf  + 


Va 


6f  6d 


~ff  + 


V0(Of-Od)  +  V(o(cQf) 


(125) 


The  tenninal  transversality  conditions  are  given  by  [25,  Algorithm  2.8]. 


(126) 


Utf)  = 


(127) 


or  more  simply  Xf  =  V 


(128) 
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Since  the  final  conditions  of  the  costates  simply  equal  new  unknown  variables, 
nothing  new  is  learned  by  analyzing  the  transversality  condition.  Since  the  final  time  is 
free,  the  Hamiltonian  Value  Condition  is  necessary.  This  states  terminal  conditions  for 
the  Hamiltonian  [25,  Algorithm  2.10]. 

H[tf)  =  - 1  (130) 

The  Hamiltonian  is  minimized  when  the  dH/du  =  0  and  d2H/du2>  0. 
Equation  (130)  provides  a  fixed  point  for  the  terminal  point  of  the  Hamiltonian.  Since 
dH(t)/du  =  0  ,  then  H(t)  =  -1  for  all  time. 

6.  Bang-Bang  Control 

Ultimately,  the  primary  lesson  learned  above  is 

w  =  -sgn(A(>max  (131) 

K  =  ~C\t  +  C2  (132) 

Equations  (131)  and  (132)  can  be  combined  as 

u  =  -sgn(-c,t  +  c2)nmax .  (133) 

Since  Xm  is  linear,  it  can  only  intersect  =  0  once.  This  may  or  may  not  occur 
within  the  range  [o,U],  therefore  it  can  be  stated  that  at  most,  only  one  sign  change 

(switch)  occurs  within  (133).  [19,  pp.  193-200]  additionally  confirms  that  only  one  sign 
change  occurs  for  initial  rest-to-rest  maneuvers.  For  arbitrary  initial  or  final  conditions,  it 
is  possible  for  no  switches  to  occur  (e.g.,  high  initial  velocity  with  constant  deceleration 
to  slow  down).  Therefore,  for  the  rest-to-rest  maneuver  in  this  example,  the  quickest 
method  of  rotation  about  this  single  eigenaxis  would  be  to  use  bang-bang  control  where 
one  switch  occurs  between  [o,U].  This  means  to  use  maximum  torque  available  along 
that  axis  until  half  way  there,  then  reverse  torque  in  order  to  slow  down  to  zero  right  at 
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the  point  in  which  the  desired  slew  angle  0d  is  achieved.  The  resulting  trajectory  is 
shown  in  Figure  26. 


Figure  26.  Acceleration  Limited  Eigenaxis  Slew 


Since  spacecraft  attitude  reorientations  may  also  be  momentum-limited  in 
addition  to  torque-limited,  large  angle  maneuvers  may  look  closer  to  that  in  Figure  29 
(heritage  maneuver  from  Chapter  III).  Momentum  is  proportional  to  angular  rate, 
therefore  the  trajectory  is  shown  as  rate-limited.  The  maneuver  is  initially  torque-limited 
Tmax  (°Cax )  ur>td  the  maximum  momentum  /imax  (o)nm )  is  reached.  During  this 
momentum-limited  section,  the  spacecraft  coasts  at  C0max ,  until  it  can  decelerate  at  its 
maximum  torque  -Tmax(-amax).  Figure  27  is  the  solution  to  the  optimal  control  problem 
with  an  additional  rate  limit. 
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Torque  /  Momentum  Limited  Maneuver 


-—o' 


Figure  27.  Rate  Limited  Eigenaxis  Slew 


B.  OPTIMAL  CONTROL  EQUATIONS  OF  MOTION 

Since  the  input  is  already  going  to  be  determined  offline  per  the  shaped  heritage 
maneuver,  the  time-optimal  solution  should  instead  be  utilized.  With  theoretically 
unlimited  computing  power  at  the  ground  station,  the  optimal  solution  can  be  computed 
for  every  individual  slew.  The  equations  of  motion  for  the  spacecraft  and  reaction  wheel 
system  from  Chapters  II  and  IV  will  now  be  consolidated.  Although  the  equations  are 
listed  in  matrix  form  for  simplicity,  the  number  of  dynamic  equations  reflects  the  number 
of  independent  equations  that  can  be  extracted  from  the  equation  set. 


1.  Spacecraft  Kinematics 

The  spacecraft  kinematics  are  repeated  from  Chapter  II.  Since  the  orbital  motion 
is  set  to  zero  (A  co  0  =  0 ,  see  Chapter  III),  then  ncjb  =  °(0B .  This  simplifies  the  dynamics 
slightly  allowing  the  removal  of  three  states.  This  change  is  reflected  immediately  in 
Equation  (134). 


$123  =  ^/®*  -Vxfe) 

CU  =-\(NCOB)T  qn 3 


(134) 


where  q 


123  ' 


q\  t/ 2  q?> 
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2. 


Reaction  Wheel  Dynamics 


The  torques  of  the  individual  reaction  wheels  follow  from  Equation  (71). 

<  =  J,nw  (135) 

Equation  (135)  can  be  solved  in  terms  of  the  reaction  wheel  acceleration. 

aw=j?  <  (136) 


3.  Spacecraft  Dynamics 

The  spacecraft  dynamics  were  developed  in  Chapters  II  and  IV. 


_ /V  j  N  -B  .  N  ^^B  .  .  j  (V_.fi 

*req=Js,c  (O+Q)  X  Js/c  CO 


-tn  =  hZ'l tw  +  Ncou  X  °Z'U  Q 

req  3  w  3  uxv 


(137) 

(138) 


Combining  Equations  (137)  and  (138)  yields  the  following. 

T  N-B  ,  N  B  .  j  N  B  brjH-W  N  B  ..  b  rjw  j  g<W 

Jjr  CO  +  CO  x7  ,  CO  =  -  Z,  T  -  CO  X  Z,  7  12 


T  N  'B  brjwW  N  ,_fi  w  /  T  (V„fi  ,  brjw  J  r^w\ 

J,c  ©  =  -  Z3  T„  -  ffl  X  ( Js/C  ©  +  Z3  ) 


(139) 

(140) 


Solving  Equation  (140)  in  terms  of  the  spacecraft  angular  velocity  rate  yields 


X'-Wxp,,. 


(V  __fi  .  b  rjw  j  r~\  W  \\ 

co  +  Z3J  J2  )) 


(141) 


4.  State  and  Control  Variables 

The  state  vector  now  has  7+n  states  and  the  control  vector  has  n  states,  where  n 
refers  to  the  number  of  reaction  wheels  in  the  array. 
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<7i 

q2 


Qa 

N..B 

6), 


N..B 

My 

N 

<°z 

w 


Q 


and  u 


(142) 


Q 


W 

n 


5.  Actuator  Constraints 


Torque  and  momentum  limitations  will  be  included  in  the  optimal  control 
problem,  as  these  constraints  represent  the  physical  limitations  of  the  attitude  control 
system. 


^  -W  ^  -W 

T  S  T  ST 

max  i  max 


i  =  1,2 


(143) 


-Q. 


<QT  <£Y  , 

max  i  max  ‘ 


i  = 


(144) 


6.  Initial  and  Final  Conditions 


The  spacecraft  will  be  conducting  rest-to-rest  slews.  This  means  initial  and  final 
quaternions  are  known.  Initial  and  final  spacecraft  rotation  rates  are  zero.  For  this 
analysis,  initial  and  final  reaction  wheel  rates  will  be  considered  zero.  In  reality, 
engineers  operate  the  reaction  wheel  at  some  nominally  non-zero  value.  This  is  because 
Q("  =0  can  cause  jitter  in  the  presence  of  “stiction”  (friction  from  static  positions)  [13, 
p.  5].  The  initial  and  final  states  can  be  represented  as  follows. 


(<?V3V<?30,<?40, 


N..B  N..B  N 

C0X  ,  CO 

*o  To 


<,q^--,q^)=(?10,?20,^,?40,o,o,o,o,--,o) 


{% >q2f  >% ’N(0xf’ N(0yf’ N(0l&Wyf  ’• ' ) =  [q{  ,q{,q{,q{ 


(145) 

(146) 
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7. 


Cost  Function 


The  objective  is  to  minimize  the  overall  maneuver  time.  Therefore,  the  cost 
function  is 

3[x(-),u(-),tf]  =  tf  (147) 


8.  Optimal  Control  Problem 

By  gathering  Equations  (134)  to  (147),  the  optimal  control  problem  can  now  be 
defined: 


Given 


x  = 


N 

<?1  <?2  <?3  <?4  (Ox 


N  B 

tu. 


N,^B  nlf 

cox  Q, 


Q 


W 


Minimize 


Subject  to  qm  =  ^[qA  N(0B  -  N (0°  x  qm) 

1  /  J V  B \T 

<74=--(  ®  )  ?123 


N(bE 

<  _  j  —  1  /  b  ryw 

~  J  s/c\  ^3 

w  N  R 

w  _  a 

w 

(*As/c 

NcaB 

+  hzu,raw)) 

nw=j;X 

’ch0 

N N  B 

,  0)x  ,  CQy  , 

xo  To 

r\W  \  (  o  o  o 

A,0 )  =  Wi  ,q2  ><h 

0,0, 0,0, •••,()) 

(% 

Wy 

.  0J,f  ’  ®v, 

'U")  =  (9if 

f  f  f 
Al 

,o,o,o,o,-,o) 

-TW 

max 

<ZW  <TW 
i  max 

-Qw 

max 

<  aw  <  o.w 

—  i  —  max 

C.  SCALING  THE  PROBLEM 


(148) 


The  purpose  of  scaling  is  to  generalize  all  states  to  within  similar  orders  of 
magnitude.  Without  scaling,  one  state  may  take  on  large  values  ( 1()6  j  while  others  are 

comparatively  small  ( 1 0  6  j .  The  relative  difference  would  be  is  (l()  12  j ,  which  causes 
numerical  problems  due  to  operating  near  machine  epsilon.  By  simply  multiplying  each 
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state  by  the  appropriate  “scaling  factor”  their  values  can  be  made  the  same  order  of 
magnitude,  and  numerical  issues  successfully  sidestepped. 

The  common  formulation  for  scaling  units  is  done  through  “canonical  units.” 
Canonical  refers  to  a  specific  definition  where,  for  example,  a  velocity  unit  (VU)  is  equal 
to  the  quotient  of  a  distance  unit  (DU)  and  a  time  unit  (TU). 

VU  =  —  (149) 

TU 

Within  [21],  “designer  units”  are  chosen  instead  to  keep  the  states  and  costates  in 
relative  orders  of  magnitude.  Instead  of  looking  at  the  costates,  nor  at  the  coupling 
between  position  and  velocity,  this  thesis  outlines  a  different  kind  of  designer  unit.  They 
will  be  denoted  by  a  superscript  *  (or  state  followed  by  “z”  within  any  listed  code). 


1.  Designer  Units 

As  scaling  factors,  the  following  designer  units  are  used  to  bring  all  variables 
close  to  the  first  order  of  magnitude. 

Wheel  torque:  r  =  Tmax  (maximum  wheel  torque  for  chosen  RW) 

Time:  T  =  Is  (use  real  time) 

Quaternion:  No  further  scaling.  Assumes  -1  <  q  <  1  is  sufficiently  scaled 

S/C  Spin  Rate:  co  =  3°  /  s  =  0.0524  rad  /  s  (slightly  faster  than  conventional  rate) 

Wheel  Rate:  Q  =  Qmm  (maximum  wheel  spin  rate  for  chosen  RW) 

For  example,  the  relationship  between  the  torque  T  and  the  scaled  torque  r  will 
be  the  proportionality  constant  T  . 

T  =T  T  (150) 

Therefore,  the  scaled  torque  would  always  operate  in  the  range  [— l,l] . 


(151) 
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Scaling  of  constants  within  a  canonical  unit  representation  can  be  useful  because 
many  terms  will  cancel  out  leaving  simpler  equations  to  manipulate  and  code.  If  all 
scaling  terms  are  directly  coded  into  the  software,  scaling  of  constants  does  nothing.  This 
is  because  the  constant  must  be  scaled  down  first,  then  scaled  back  up  to  its  original  value 
within  the  dynamic  equation.  In  the  case  of  designer  units  where  many  mismatched  units 
are  being  manipulated,  it  is  unlikely  that  convenient  units  will  cancel  leaving  an  elegant 
equation.  Therefore,  scaling  is  not  required  and  not  completed  for  constants.  An  example 
of  designer  unit  scaling  is  as  follows  for  the  spacecraft  dynamics.  The  full  set  of  scaled 
equations  is  shown  in  Appendix  A. 


N(bB 


t. 


br*w-w  N(oBx(j 


N(0B  +  b ' 


NcoB 


f  CO 

j 


=  j  1 

u  sic 


z;j„s2wj) 


(152) 


fe oB  =  J 
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f 

b  ry\\>—W 

—  Zo  T 

3  "  w 

rz*T*' 

V 

{  CO*  ) 

'a,(T’)x(j„/m,(o>')+bz;j,.aw(Sl')) 


2.  Scaled  Optimal  Control  Problem 

The  scaled  optimal  control  problem  is 
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Given 


x  = 


u  = 


—  —  —  —  n  —d  n  —  b  /V  —  o  rrw 

<h  <h  <?3  94  ®,  Q!  ■■■  Q„ 

nr 


Tf 


t5 


Minimize  S[x  '(•), «(•),?/]  =  tf 

Subject  to  qn3  =  ^(q4  N  cdB  -  N (QB  x  ^123)( 7*© *) 


^4  =  -i(W)r?I23(rV) 


N  tzB  t-1 

®  =4 


6  f  r*r*  3 

bryW—W  1  1 

-  Z  Th,  — 

V  ®  y 


v  (r* )  x  (ys/c  "©» hz;jKQw  (Q*)) 


=y;'T 
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Vr*^ 


o 
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A r—B  N  rrrB  F\W 

CO  .  CO  .£2,  .• 

<h_ 

To  ' 

to'  l0  ' 

’  "o  / 

\q' 

( _/ 

n  ttB  n  F\w 

CO  .  CO  .£2,  .* 

••,QVV)  = 
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’  nf) 

,q 

q  q  q  q 


-  1  <  fW  <  1 
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D.  DIDO  IMPLEMENTATION 

DIDO  is  a  MATLAB-based  software  package  for  solving  optimal  control 
problems  [21].  A  huge  advantage  of  DIDO  is  that  the  optimal  control  problem  from 
Equation  (153)  is  the  only  necessary  derivation.  Since  Pontryagin’s  Principle,  discussed 
in  Section  A,  is  automatically  integrated  by  the  DIDO  software,  further  analysis  is  not 
required.  But,  Pontryagin’s  Principle  can  be  checked  to  verify  the  numerical  solution. 
The  code  simply  requires  the  optimal  control  problem  to  be  inputted  into  a  series  of  four 
functions  run  from  a  single  m-filc. 

1.  Cost  Function 

The  cost  function  comes  from  Appendix  A,  Equation  (275),  implemented  in 
Figure  28. 

3  [*(•),  «(•),*/]=  tf  (154) 
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%  Cost  function  for  the  Min  Time  Slew  using  RWs 

function  [ endpointCost,  runningCost]  =  MinTimeSlewRWCost ( primal) 

tf  =  primal . nodes ( end) ; 

endpointCost  =  tf; 
runningCost  =  0; 


Figure  28.  Cost  Function  Implemented  as  a  DIDO  Function 


2.  Endpoint  Function 

The  endpoint  function  contains  the  boundary  conditions  from  Equations  (280)  and 
(281),  implemented  in  Figure  29. 


( 77  77  77  77  N  ttB  N  ttB  N  ttB  nlf  r\W  \  _ 

(—  —  —  —  NttB  n —b  n—b  \ 

y-i\f  ’(h)  .’fof  ’cl\f  ’  ®xf’  ® yf ’  — 


f  0  0  0  0  '\ 

yq  q  q  q  j 

«L  «L  0  0  0  0..V 

*9  *  9  *  9  *  )U9V,9U)  w9  9 W 

yq  q  q  q  ) 
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%  Events  function  for  the  Minimum  Time  Slew  using  RWs 
function  endpointFunction  =  MinTimeSlewRWEndpoint( primal) 

global  CONSTANTS 

Zn  =  CONSTANTS. Zn  ;%  number  of  RWs 

qO  =  primal . states ( 1 : 3 , 1 )  ;  qf  =  primal . states ( 1 s 3 , end) 

q40  =  primal . states ( 4 , 1 )  ;  q4f  =  primal . states ( 4 , end) 

wO  =  primal . states ( 5 s 7 , 1 )  ;  wf  =  primal.states(5:7,end) 

WO  =  primal . states ( 8 : 7+Zn, 1 ) ;  Wf  =  primal. states ( 8 : 7+Zn, end) 

endpointFunction  =  [ qO ;  q40;  wO;  WO;  qf;  q4f;  wf ;  Wf]; 


Figure  29.  Endpoint  Function  Implemented  as  a  DIDO  Function 


3.  Path  Function 

The  path  is  not  necessary  for  this  problem,  so  no  path  function  script  is  needed. 

4.  Dynamics  Function 

The  dynamic  equations  are  from  Equations  (276)  through  (279). 
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Within  the  dynamics  function  shown  in  Figure  30,  state  scaling  occurs  at  the  point 
of  state  assignment  and  is  then  accounted  for  in  the  equations  of  motion. 
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%  Dynamics  function  for  the  Minimum  Time  Slew  using  RWs 
function  xdot  =  MinTimeSlewRWDynamics (primal) 


global  CONSTANTS 
Js  =  CONSTANTS. Js 

Jw  =  CONSTANTS . Jw 

Z  =  CONSTANTS . Z 

Zn  =  CONSTANTS. Zn 


;%  spacecraft  inertia  matrix 
;%  wheel  inertia  matrix 
;%  RW  orientation 
; %  number  of  RWs 


Tauz 

=  CONSTANTS. Tauz 

Wz 

=  CONSTANTS. Wz 

wz 

=  CONSTANTS. WZ 

Tz 

=  CONSTANTS. Tz 

qz 

=  CONSTANTS. qz 

; %  torque  unit 
;%  wheel  rate  unit 
;%  wheel  rate  unit 
; %  time  unit 
; %  quaternion  unit 


%  States  and  controls 


q 

=  primal . states ( 1 : 3 ,:) *qz 

;% 

quaternions  1-3 

q4 

=  primal . states ( 4 ,:) *qz 

;  % 

quaternion  4 

w 

=  primal . states ( 5 : 7 ,:) *wz 

;% 

spacecraft  rate 

W 

=  primal . states ( 8 : 7+Zn, :) *Wz 

;% 

wheel  rates 

tau 

=  primal. controls*Tauz 

;% 

wheel  torques 

N  =  length(q4);  %determine  size  and  preallocate  matrix  sizes 

qdot  =  zeros(4,N);  q4dot  =  zeros(l,N)  ; 
wdot  =  zeros ( 3, N);  Wdot  =  zeros (Zn,N); 


%  Equations  of  Motion: 
for  i  =  1:N 


qdot(:,i)  =  l/2*(  q4(i)*w(:,i)  -  cross(w( : ,i) ,q(  :  ,i)  )  )*Tz/qz  ; 

q4dot(:,i)=  -l/2*w( : , i) ' *q( : , i) *Tz/qz  ; 

wdot(:,i)  =  Js\eye( 3 ) * ( -Z*tau( : , i)  ... 

-cross (w(:,i),Js*w(: , i) +Z* Jw*W( : , i) )  ) *Tz/wz ; 

Wdot(:,i)  =  Jw\eye( Zn) *tau( : , i)  *Tz/Wz; 

end 

xdot  =  [qdot  ;  q4dot;  wdot  ;  Wdot]; 


Figure  30.  Dynamics  Function  Implemented  as  a  DIDO  Function 


5.  Main  Function 

With  only  a  few  dozen  lines  of  code,  the  main  m-lilc  runs  the  four  functions  and 
DIDO.  This  is  where  the  CONSTANTS  are  defined.  Boundaries  are  set  up  for  the  states, 
controls,  and  endpoint  conditions  (see  Figure  31). 
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tfMax  =  40 


;%[sec]  upper  bound  for  scaled  problem 


bounds . lower • time  =  [0;  0]; 

bounds. upper. time  =  [0;  tfMax]; 

%  Scaled  initial  conditions 

q0  =  q_0/qz  ; 

wO  =  [  0 ; 0 ; 0 ] /wz  ; 

WO  =  zeros ( Zn, 1 ) /Wz  ; 

hO  =  Jw*W0  ; 


%  final  conditions 

qf  =  q_f/qz 
wf  =  [0;0;0]/qz 
Wf  =  zeros ( Zn, 1 ) /Wz 

%  boundary  conditions  [  ql-q4  ;  wx;wy;wz 

bounds . upper . states  =  [ ones ( 4 , 1 ) /qz ;  ones (3,1) 
bounds . lower . states  =  -bounds. upper. states; 

bounds . upper . controls  =  ones ( Zn, 1 ) *tau_max/Tauz ; 
bounds. lower. controls  =  -bounds. upper. controls; 

bounds. upper. events  =  [q0;  wO ;  WO;  qf;  wf;  Wf];  %(20xl) 
bounds . lower . events  =  bounds. upper. events; 

%  Tell  DIDO  the  bounds  on  the  problem 

MinTimeSlewRW. bounds  =  bounds; 


Wl-Wn  ] 

ones ( Zn, 1 ) *W_max/Wz  j ; 


Figure  3 1 .  DIDO  Bounds  Implemented  within  a  DIDO  script 


DIDO  is  initially  run  with  no  guess  and  as  little  as  20  nodes  to  detennine  if  a 
feasible  solution  will  arise  (see  Figure  32). 


%%  RUN  1  -  NO  GUESS  MADE 

%  Number  of  nodes  for  the  spectral  algorithm 

algorithm. nodes  =  20; 

[cost,  primal,  dual]  =  dido(MinTimeSlewRW,  algorithm); 

Figure  32.  DIDO  No  Guess  (20  nodes)  Implemented  within  a  DIDO  script 

The  results  are  then  seeded  back  into  DIDO  as  initial  guesses  along  the  entire 
dynamic  range  using  the  following  code.  This  time  a  much  more  refined  set  of  120  nodes 
is  used  (see  Figure  35) 
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%%  RUN  2  -  SEED  RESULTS  OF  RUN  1  AS  GUESS 
algorithm. nodes  =  120; 
algorithm. guess  =  primal; 

[cost,  primal,  dual]  =  dido(MinTimeSlewRW,  algorithm); 

Figure  33.  DIDO  Seeded  Guess  (120  nodes)  Implemented  within  a  DIDO  script 

This  completes  all  steps  required  to  run  DIDO  under  both  the  “no-guess”  and 
“seeded  guess”  scenarios.  The  results  from  run  1  can  be  looked  at  as  an  initial  estimate  to 
put  DIDO  on  track. 

E.  OPTIMAL  MANEUVER  IMPLEMENTATION 

The  optimal  control  solution  is  fundamentally  open  loop  and  the  control  solution 
is  applied  directly  to  the  individual  wheels  (see  Figure  34).  Additionally,  the  optimal 
quaternion  path  is  presented  as  a  shaped  angular  trajectory  at  the  forefront  of  the  loop. 
Since  body  rates  are  known,  the  EQF  can  follow  the  rate  trajectory  instead  of  being 
regulated  to  zero  angular  velocity.  When  plant  uncertainties  or  external  disturbances 
exist,  the  instantaneous  quaternion  error  generates  work  for  the  EQF  Controller,  applying 
any  additional  wheel  torque  as  required.  With  these  changes,  Figures  2  and  17  become 
Figures  34  and  35,  respectively. 


Orbital  Motion 


Figure  34.  Architecture  for  Optimal  Maneuver  Implementation:  Spacecraft  Model 
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Z (3xn) 


i2-' 


r"  ,  Z‘  (-tN  -  NuB  X  Zhw) 

rw  rw  req  rw 


r"  =  Z  rw  +  Nu.BXZhw 


Figure  35.  Architecture  for  Optimal  Maneuver  Implementation:  Reaction  Wheel 

Subsystem 


Figure  36  shows  the  integration  of  the  rate  error  into  the  eigenaxis  quaternion 
feedback  controller,  a  slight  modification  from  Equation  (28). 

^req  —  J s/c  Qei23  ~  C  J s/c  "F  ®  ^  J  sic  ®  ( 


Figure  36.  Architecture  for  Optimal  Maneuver  Implementation:  Eigenaxis  Quaternion 
Feedback  Controller 
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By  modifying  the  control  architectures  as  per  Figures  34  to  36,  the  full  reaction 
wheel  array  torque  and  momentum  is  available  for  use.  If  the  state  quaternions  q  feed 
back  disturbances  and/or  uncertainty  error,  quaternion  error  is  corrected  by  creating  the 
quaternion  error  from  the  optimal  quaternion  trajectory.  The  EQF  Controller  then  creates 
a  torque  requirement.  Body  rate  error  is  propagated  through  both  the  EQF  controller  and 
wheel  dynamics  in  the  fonn  of  additional  required  torque.  The  total  required  torque  x'lq  is 

then  passed  through  the  pseudo-inverse  to  sum  with  the  optimal  torque  . 


T 


w 

W 


(158) 


The  system  will  naturally  require  some  correction  torque  t*  ,  because  the 

spacecraft  inertia  and  external  forces  (if  present)  cannot  be  modeled  perfectly.  This  leads 
to  the  requirement  for  the  pseudo-inverse  to  perform  the  allocation  of  the  feedback 
torques.  Since  the  pseudo-inverse  uses  the  minimum  least-squares  solution  for  converting 
body  torque  into  n  wheel  torques,  the  feedback  provided  follows  a  minimum  energy 
control  allocation. 


To  accommodate  the  need  for  feedback,  it  is  necessary  to  reduce  the  control 
authority  of  the  optimal  control  solution.  This  is  because  if  the  optimal  torque  solution 

T™  is  already  utilizing  the  maximum  system  torque  r^ax ,  the  system  will  not  have  the 

authority  needed  to  make  the  feedback  corrections.  The  available  control  authority  for 
optimal  control  should  be  reduced  inversely  proportional  to  the  knowledge  of  the  system. 
If  the  system  knowledge  is  high,  then  the  control  authority  should  be  reduced  only 
slightly.  For  instance,  if  the  system  knowledge  is  high,  an  optimal  control  authority  of 
95%  may  be  appropriate.  If  the  system  knowledge  is  low,  an  optimal  control  authority  of 
75%  may  be  appropriate.  Within  the  optimal  control  problem  formulation  presented  here, 
the  adjustments  would  resemble  that  in  Equation  (159). 


-0.95t1<tw  <0.95  r 


w 

max 


(159) 
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F.  EXAMPLE  TIME-OPTIMAL  MANEUVER 

In  this  section,  the  feasibility  and  optimality  of  a  representative  time-optimal 
example  will  be  documented.  As  a  representative  generic  spacecraft,  the  following 
figures  will  be  used. 

-w 

0  -V2  0 

V2  0  ->/2  (160) 

1  1  1 


2.54  0  0 


0  2.54  0  kg-m 2 

(161) 

0  0  2.54 

C„.max  =  ±8.57m/V-m 

(162) 

=  ±  3000  rpm  =  +  3 14.2  rad  /  s 

(163) 

jw  0  0  0 

Jwl  33 

o  L  o  o 

J  w2  33 

0  0  jw  0 

0  0  0  j  w 

J  w 4  33 

The  designer  units  selected  are: 

Wheel  torque:  x *  =  ±  8 .57  mN  ■  m 

Time:  T*=lsec 


(164) 


Quaternion:  q  =  1  (no  further  scaling) 

S/C  Spin  Rate:  oj  =  3deg/  5  =  0.0524 rad  /  s 

Wheel  Rate:  Q  =  3 14.2  rad  /  s 
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A  120°  rest-to-rest  slew  will  be  performed  about  the  eigenaxis 


e  = 


V2V3  0  VT 73 


from  q°  = 


-Vl 72  0  -1/2  1/2 


-\T 


to  the  origin 


qf  =  [  0  0  0  1  ]  .  The  boundary  conditions  are  shown  in  Equation  (165). 

K>$20.?V?4o>  03 03 '.v0’  ’^20  ’^30  ’^40  j 

=  (-Vl72,  0  ,  - 1/2,  1/2,0,  0,  0,  0,  o,  o,  oj 


(165) 


:(  0  ,  0  ,  0  ,  1  ,0,0, 0,0, 0,0,0) 


1.  Optimal  Solution 

The  time  history  of  the  control  allocation  determined  by  DIDO  is  shown  in  Figure 
37.  For  each  of  the  120  time  nodes,  DIDO  determined  the  optimal  control  solution  for 
each  wheel.  Linear  interpolation  was  used  to  determine  the  feedforward  control  allocation 
that  would  be  used  for  the  simulation. 


RW  Torque  (mN-m-s)  for  [-0.7071 1 ,0, -0.5, 0.5]  to  [0,0,0,1]  slew  (120  nodes) 
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Figure  38  shows  the  time-optimal  quaternion  trajectory.  The  final  node  exists  at 
qf  =  [  0  0  0  l]  as  expected  in  the  boundary  conditions. 


0  .4  - 
0.2  - 
Oi 
-0.2 
-0.4 
0 

-0.2 

CO 

^  -0.4 
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-0.6 
1.2  F 
1  - 


0.6( 

0.4* 


Quaternions  for  |-0.7071 1,0, -0.5, 0.5]  to  [0,0, 0,1]  slew  (120  nodes) 


Mil  in 


ammwwwBiin—  i 
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Figure  38.  Optimal  Quaternion  Trajectory  (120  Nodes) 

Additionally,  the  final  nodes  for  the  angular  rates  all  end  at  zero  per  the  boundary 
conditions. 

2.  Feasibility 

A  feasible  control  solution  is  one  that  propagates  through  the  system  dynamics  as 
expected  and  satisfies  the  boundary  conditions.  The  feasibility  of  both  solutions  is 
checked  by  propagating  DIDO’s  control  solution  through  the  system  dynamics.  First,  the 
state  values  at  each  of  the  individual  time  nodes  must  align  with  the  propagation  of  the 
dynamics  at  those  specific  times.  Second,  the  dynamics  must  propagate  to  the  appropriate 
final  conditions.  If  both  of  these  have  been  satisfied,  then  it  is  detennined  that  enough 
nodes  have  been  used  and  the  solution  is  indeed  feasible. 
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The  interpolated  control  signal  (Figure  37)  is  propagated  through  the  spacecraft 
dynamics,  modeled  within  Simulink  per  Figures  34  to  36.  For  DIDO  propagation,  EQF  is 
turned  off.  The  DIDO  torque  signal  is  fed  into  the  reaction  wheel  dynamics  for  open-loop 
propagation.  Since  the  knowledge  of  the  spacecraft  inertia  is  considered  perfect,  the 
feedback  portion  will  only  affect  the  conventional  slew,  later  used  for  comparison. 

States:  Once  propagation  is  complete  for  the  full  120  nodes,  the  control  and  states 
for  both  runs  and  their  respective  control  propagations  can  be  plotted  together.  If  the 
curves  for  each  propagated  state  pass  through  its  respective  DIDO  nodes,  the  DIDO- 
produced  control  solution  can  be  deemed  feasible. 

Figures  39-41  show  is  a  graphical  comparison  of  the  initial  DIDO  solution  vs.  the 
dynamically  propagated  solution  from  Simulink  for  all  1 1  states.  By  inspection,  it  appears 
that  every  single  numerically  calculated  point  is  directly  on  top  of  its  propagated 
trajectory.  This  shows  that  the  DIDO  solution  is  not  only  feasible,  but  also  very  accurate. 


Quaternions  for  [-0.7071 1,0, -0.5, 0.5]  to  [0,0, 0,1]  slew  (120  nodes) 
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t#BN  (deg/s)  for  [-0.7071 1 ,0, -0.5, 0.5]  to  [0,0, 0,1  ]  slew  (120  nodes) 


RW  Rate  (rpm)  for  [-0.7071 1 ,0, -0.5, 0.5]  to  [0,0, 0,1  ]  slew  (120  nodes) 
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Scaled  States  and  Costates:  Another  important  check  to  detennine  the  efficiency 
of  code  is  the  relative  magnitude  of  the  scaled  states  and  costates.  In  Figures  42  and  43, 
the  scaled  states  and  costates  (3c, A)  are  represented  in  the  selected  designer  units  prior  to 
their  rescaling  back  into  SI  units.  All  states  and  controls  are  on  the  same  order  of 
magnitude,  which  demonstrates  good  scaling.  Only  a  1  order  of  magnitude  difference 
shows  up  in  the  costates,  which  are  very  sensitive  to  scaling.  The  relative  comparison 
shows  that  reasonable  scaling  was  achieved  during  problem  formulation.  If  the  costates 
were  to  differ  from  each  other  by  1-2  orders  of  magnitude,  scaling  of  the  time  constant 
T*  or  turning  the  knobs  on  the  state  parameters  could  adjust  this.  Additionally,  a  key 
point  to  note  is  that  the  scaled  states  only  differ  from  their  costates  by  one  to  two  orders 
of  magnitude.  A  well-scaled  problem  will  have  them  roughly  the  same  order  of 
magnitude  [21,  p.  33],  therefore  tweaking  of  the  designer  units  could  remedy  this. 

1 

cr  0 

-1 

1 

3  0 

-1 

1 

CS  0 

-1 

0  5  10  15  20  25  30  35  40 

time  (sec) 

Figure  42.  DIDO  Scaled  States 
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Scaled  Costates  for  [-0.7071 1,0, -0.5, 0.5]  to  [0,0, 0,1]  slew 


Figure  43.  Scaled  Costates 


3.  Optimality 

The  first  place  to  check  for  optimality  is  the  Hamiltonian.  To  ensure  the  optimal 
solution  is  found,  the  Hamiltonian  must  be  held  constant  with  respect  to  time. 
Additionally,  the  final  time  is  free,  so  the  Hamiltonian  Value  Condition  is  necessary, 
refer  back  to  Equations  (124)  through  (130)  from  Section  A.  Similarly,  it  can  be  shown 
that 

(166) 

Equation  (166)  reveals  a  key  point;  for  all  time-optimal  cases,  the  Hamiltonian  is 
held  constant  at  -1.  Figure  44  is  a  zoomed-in  plot  that  shows  that  the  time-history  of  the 
Hamiltonian  on  approximately  -1  ±  .007 ,  which  is  reasonable. 
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Hamiltonian 


Figure  44.  Hamiltonian 


Although  DIDO  does  not  require  inputs  from  the  optimality  analysis  in  Section  A 
to  run,  this  information  can  help  further  analyze  the  optimality  of  the  solution. 
Pontryagin’s  switching  functions  are  easily  derived.  For  the  four  control  torques,  the 
switching  functions  are  as  follows. 

s,  =  4,  / J. -  VT73 X„,  /  j,  - sl2/3 XM  /  j, 
s,  =  K 1 1„  -  VT73  X„ !  j,  -  -Jlii  XM  /  j, 

I —  , —  (16y) 

S ,  =  Xa>  /  jw  i,  -  x„  / 1 

s.  =  4.  /  iw  - Vu3  xml  /  j, + -Jin  x.2  /  j, 

where  S  =  A-  (168) 

Per  Pontryagin’s  Principle  (refer  to  Section  A),  these  switching  functions 
represent  the  slope  of  the  Hamiltonian  with  respect  to  the  control  used.  This  defines  the 
control  space  by  further  determining  what  control  effort  should  be  used  at  any  instant  in 
time.  If  the  slope  is  negative,  i.e.,  5,  <  0  ,  the  maximum  torque  should  be  used  (minimum 

torque  for  positive  slope).  When  the  slope  =  0,  the  value  of  the  control  does  not  really 
affect  the  Hamiltonian.  The  following  switching  logic  represents  the  optimal  use  of  the 
controls. 


If  SA 


<0  Tl  —  ^  max 

=  0  ,  then  use  tmm  <  r,  <  T„ 
>0  T;  =  T,„;„ 


(169) 
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Figure  45  shows  the  dynamics  of  the  reaction  wheel  assembly  alongside  the 
solution  for  the  switching  function  at  each  of  the  120  nodes.  The  torque  responds  very 
closely  to  the  switching  function  logic  in  Equation  (169),  thereby  further  validating  the 
optimality  of  the  solution. 


RW  Torque  (mN-m-s)and  Switching  Function  (120  nodes) 


G.  EIGENAXIS  SLEW  COMPARISON 

For  comparison  purposes,  the  Eigenaxis  Quaternion  Feedback  (EQF)  Controller  is 
used  to  perfonn  a  conventional  eigenaxis  slew  maneuver  via  a  shaped  bang-coast-bang 
trajectory.  EQF  is  governed  by  Equation  (170),  repeated  from  Chapter  III,  Equation  (28). 

=  ~kJqem -cJffl+Vx/j/cV  (170) 

For  this  example,  a  0.1  sec  settling  time  was  chosen.  Recall,  the  step  size  is  a 
function  of  the  distance  between  discrete  points  used  within  the  shaped  quaternion 
trajectory.  Therefore,  a  small  settling  time  can  now  be  used.  The  damping  ratio  was 
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chosen  as  0.9  so  that  it  does  not  allow  for  much  overshoot.  Using  these  criteria,  the  gains 
k  and  c  can  be  found. 

°-lsec 
£  =  0.9 

®„=4/(wO  =  44.4  (171) 

k  =  co2n  =  \  .91 E3 
c  =  2C(y„  =  80 

A  sanity  check  can  be  done  by  comparing  the  time-optimal  slew  against  the 
classic  EQF  mentioned  above.  Figure  46  demonstrates  the  different  paths  taken  by  each 
slew.  Notice  that  q2  changed  for  the  optimal  slew  only.  Acceleration  also  occurred  faster 
for  the  optimal  slew.  A  detailed  analysis  of  how  acceleration  and  momentum  are  utilized 
differently  occurs  within  Chapter  VII. 


time  (sec) 

Figure  46.  Optimal  and  Eigenaxis  Comparison:  Quaternions 

Figure  47  represents  this  comparison  with  respect  to  the  Euler  angles.  The  DIDO 
run  took  only  40.5  seconds,  21.8%  faster  than  the  51.8  second  EQF  slew.  It  does  not 
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matter  which  axis  is  used,  the  time-optimal  slew  will  always  be  just  as  fast  or  faster  than 
the  eigenaxis  slew. 


EQF  Slew  Comparison  (Optimal  Slew  is  21 .8%  faster) 


Figure  47.  Optimal  and  Eigenaxis  Comparison:  Euler  Angles  and  Rates 


The  optimal  maneuver  achieved  the  final  boundary  conditions  much  faster  for  two 
reasons.  It  was  able  to  take  advantage  of  off-eigenaxis  motion,  achieving  much  greater 
magnitudes  in  both  torque  and  momentum.  Figure  48  illustrates  the  additional  torque 
(steeper  slope)  and  momentum  (maximum  rate).  The  body  rate  of  the  eigenaxis  slew 
confined  by  the  pseudo-inverse  was  restricted  to  3°/s.  On  the  other  hand,  the  optimal 
maneuver  could  reach  momentum  and  torque  outside  the  pseudo-inverse  space,  obtaining 
body  rates  of  up  to  4.24°/s,  a  41.5%  increase  in  maximum  rate. 
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Figure  48.  Optimal  and  Eigenaxis  Comparison:  Body  Rates 


Table  2  presents  a  maneuver-time  comparison  of  a  few  more  120°  slews.  Each  is 
moved  from  the  initial  quaternion  back  to  the  origin  [0,0,0, 1].  As  the  initial  quaternion 
changed,  eigenaxis  slew  rate  increased.  This  is  because  the  eigenaxis  was  limited  only  by 
the  pseudo-inverse-limited  octahedron,  not  the  spherical  surface  (see  Figure  24).  As  the 
eigenaxis  approached  a  maximum  vertex  in  run  2,  the  optimal  slew  only  improved  the 
time  by  4.2%.  Interestingly,  the  time-optimal  maximum  body  rate  was  consistently 
between  4.0-4.3°/s,  but  this  accounts  for  much  off-eigenaxis  movement. 


q0 

Eigenaxis 

Eigenaxis  Slew 
time  wmax 
(sec)  (0/s) 

Optimal  Slew 
time  wmax 
(sec)  (7s) 

time 

improvement 

[-0.707,0,-0.5, -0.5] 

[  0.8165,  0.0000,  0.5774] 

51.8 

3.0 

40.5 

4.2 

21.8% 

[-0.5, -0.5, -0.5, 0.5] 

[  0.5774,  0.5774,  0.5774] 

44.0 

3.7 

42.1 

4.3 

4.2% 

[-0.612,-0.433,-0.433,0.5] 

[0.7071,0.5000,0.5000] 

46.7 

3.5 

41.5 

4.0 

11.0% 

[0.612,0,0.612,0.5] 

[-0.7071,  0.0000,-0.7071] 

51.2 

3.1 

39.8 

4.3 

22.3% 

Table  2.  Optimal  vs.  Eigenaxis  Comparison  (120°  slews) 


H.  EFFECTIVE  EIGENAXIS  OF  TIME-OPTIMAL  MANEUVER  - 
ANALYTICAL  APPROACH 

It  is  important  to  introduce  the  concept  of  effective  eigenaxis.  Given  any  two 
quaternions,  the  optimal  slew  deviates  off-eigenaxis,  varying  its  rotation  axis  and  rate 
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throughout  the  maneuver.  The  effective  eigenaxis  is  defined  as  the  axis  of  rotation  for  a 
minimum  angle  slew  (original  eigenaxis  from  qo  to  qf),  regardless  of  path  actually 
traveled  during  the  maneuver.  For  example:  for  Table  2,  run  1,  the  effective  eigenaxis  of 
the  optimal  slew  was  [0.8165,0,0.5774].  The  actual  distance  traveled  by  the  time-optimal 
slew  was  123°  versus  the  minimum  angle  maneuver  (120°)  about  the  eigenaxis. 

Using  this  concept,  it  is  possible  to  determine  the  relationship  between  the 
effective  eigenaxis  and  the  effective  eigenaxis  torque  and  slew  rate  about  that  effective 
eigenaxis,  the  motivation  for  the  next  two  sections.  This  is  necessary  to  form  a  general 
comparison  between  the  two  very  different  maneuvers.  The  crucial  advantage  of  the 
optimal  maneuver  is  that  it  can  access  all  of  the  dodecahedron  momentum  space, 
including  that  outside  of  the  pseudo-inverse  octahedron  (Figure  49).  This  is  because  the 
optimal  torque  trajectories  for  each  wheel  are  obtained  as  part  of  the  solution  to  the 
optimal  control  problem.  Flence  the  pseudo-inverse  need  not  be  employed  (recall  control 
law  from  Figure  34).  This  momentum  envelope  minimum  begins  at  1.633  hmm  at  the 
center  of  the  twelve  flat  faces,  expanding  to  2.309  hmax  at  the  maximum  vertices  (1.633 
Unax  to  2-309  Tmax  in  torque  space). 


-i  o  1 

A 


Time  Optimal 
can  access 
these  previously- 
unavailable 
regions 


b2 


Figure  49.  Time-Optimal  Accessible  Regions 
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It  is  difficult  to  compare  a  generic  optimal  maneuver  against  an  eigenaxis  slew, 
primarily  because  the  precise  shape  of  the  optimal  solution  is  not  intuitive.  The  optimal 
maneuver  is  tailored  to  maximize  performance  of  an  individual  maneuver.  Thus,  a  120° 
slew  about  one  effective  eigenaxis  might  be  far  different  than  a  120°  slew  about  a 
different  effective  eigenaxis  (see  Table  2  for  slew  time  comparison).  Looking  at  Figure 

50,  the  [  0  1  0  ]  (or  b2 )  axis  yields  a  limit  at  1.633  /?max  ( Trnax  ),  while  the  eigenaxis 
about  [  0  0  1  J  (or  b, )  yields  2.309  hrmx  ( Trnax  ).  This  means  that  the  rotation  about 

the  b  ,  axis  can  potentially  accelerate  and  traverse  41%  faster  than  a  rotation  about  b2  . 

Here,  the  effective  eigenaxis  momentum  will  be  derived  from  the  average 
eigenaxis  momentum  over  47t  steradians.  The  first  step  will  be  to  calculate  the  average 
amplitude  of  momentum  or  torque  available  about  every  possible  eigenaxis.  The  average 
amplitude  of  a  series  of  any  set  can  be  determined  by  the  sum  of  the  set  divided  by  the 
number  of  entries  in  the  set.  Therefore,  the  average  of  magnitude  of  n  (100,  1000,  1M) 
eigenaxis  momentum  vectors  is  the  quotient  of  the  sum  and  n. 


If  h  was  placed  back  into  3D  space  along  the  same  n  eigenvectors,  the  shape 

would  now  be  a  sphere  of  radius  havg .  Consider  each  vector  to  have  some  volume 

associated  with  it.  Since  the  total  length  of  the  sum  is  conserved,  the  sphere  would  have 
the  same  volume  as  the  original  dodecahedron. 


This  approach  will  be  followed  in  reverse  in  order  to  attain  havg ,  the  radius  of  the 

sphere.  Using  MATLAB’s  convhull  command,  the  volume  of  a  shape  is  calculated  as  a 
subset  of  its  three-dimensional  points.  Using  this  command,  the  volume  of  the 
dodecahedron  is  calculated  as  7i4.632  hmax3 .  A  simple  spherical  relationship  yields  the 
radius  of  the  sphere. 


(173) 
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Now,  the  equivalent  spherical  radius  of  the  reaction  wheel  momentum  (or  torque) 
envelope  is  1.805.  The  shape  is  shown  in  Figure  50. 


h  ,  h  =1.805  h 

opt- sphere  ir 

^  opt -sphere  =  ^  ‘^05  Tn 


(174) 


Since  the  sphere  is  not  inscribed  within  the  dodecahedron,  it  would  appear  that 
some  torque/momentum  described  by  the  sphere  is  not  available.  This  may  not  be  the 
case  though  as  the  sphere  simply  represents  the  effective  torque/momentum  about  the 
effective  eigenaxis.  Since  the  trajectory  is  allowed  to  deviate  from  the  effective  eigenaxis, 
the  torque/momentum  space  at  the  vertices  is  exploited  to  raise  the  effective 
acceleration/rate  limits  outside  the  boundaries  of  the  torque/momentum  envelope. 
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Consider  the  optimal  trajectory  using  only  eigenaxis  slews,  but  has  access  to  the 
entire  momentum  space.  The  average  slew  rate  (or  acceleration)  derives  directly  from  the 
newly-created  optimal  momentum  sphere  with  values  per  (174).  Since  both  torque  and 
momentum  limitations  are  assumed,  the  average  time  for  any  given  slew  angle  can  be 
predicted  using  the  simple  trapezoidal  trajectory:  bang-bang  (from  Figure  26)  or  bang- 
coast-bang  (from  Figure  27),  depending  on  whether  the  momentum  limit  is  reached.  The 
same  spacecraft/wheel  configuration  from  Section  G  will  be  used  for  comparison.  Over 
the  slew  regime  extending  to  180°,  the  curve  shown  in  Figure  5 1  is  developed. 


Figure  5 1 .  Reaction  Wheel  Effective  Eigenaxis  for  Time-Optimal  Maneuver  -  Analytical 

(js=25.4  kg-m2) 


The  curve  at  the  beginning  indicates  a  torque  limited  region.  For  slew  angles  of 
30°  or  less,  the  reaction  wheel  array  does  not  have  the  torque  to  reach  its  full  momentum. 
Therefore,  the  maneuver  is  bang-bang.  It  will  be  in  maximum  torque  mode  (edge  of  the 
torque  sphere,  0.1 805 Am),  then  reverse  maximum  torque  to  slow  down  without  ever 
reaching  the  edge  of  the  momentum  envelope.  For  slews  greater  than  30°,  the  maximum 
momentum  is  reached  and  the  system  theoretically  coasts  at  a  constant  rate  (proportional 
to  1 .805  Nms  )  for  a  portion  of  the  maneuver. 
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In  this  section,  the  concept  of  effective  eigenaxis  was  used  as  a  baseline  to 
compare  the  full  spectrum  of  possible  optimal  maneuvers  against  the  predictability  of 
shaped  eigenaxis  slews.  The  relationship  between  effective  eigenaxis  torque/momentum 
for  optimal  maneuvers  was  developed.  For  the  optimal  maneuver,  an  effective 
momentum  of  1 .805  hmax  and  torque  of  1 .805  Tmax  can  be  used  to  represent  the  potential  of 
the  time-optimal  maneuvers  over  47t  steradians.  Additional  comparisons  will  be  carried 
out  in  Chapter  VII. 


I.  EFFECTIVE  EIGENAXIS  OF  TIME-OPTIMAL  MANEUVER  - 
NUMERICAL  VERIFICATION  OF  ANALYTICAL  APPROACH 


The  analytical  approach  from  Section  H  must  be  validated.  To  gather  sample 
data,  100  optimal  simulations  were  completed  using  the  spacecraft/wheel  configuration 
from  Section  G.  The  only  dependent  variable  used  was  the  beginning  attitude.  A  random 
quaternion  was  generated  for  each  trial,  and  the  system  slewed  to  the  origin. 


q0  =  2  (rand  (4 ,1)  -  0 .5) 

qf  =  [  0  0  0  1  J 


(175) 


This  allows  a  random  sampling  of  slew  angles  and  eigenaxes.  The  sample  data  is 
plotted  alongside  the  analytical  curve,  shown  in  Figure  52. 
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Figure  52.  Reaction  Wheel  Effective  Eigenaxis  for  Time-Optimal  Maneuver  -  Numerical 

Approach  with  1.807  xmax  (js=25.4  kg-nr) 


The  majority  of  the  sample  data  lies  above  the  curve,  demonstrating  that  the 
approach  taken  in  Section  H  might  be  too  aggressive.  The  slope  of  the  momentum- 
limited  region  seems  correct,  but  the  intercept  with  the  torque-limited  region  is  too  early. 
Therefore,  the  effective  torque  bound  must  be  reduced.  Several  values  for  the  effective 
torque  are  shown  in  Figure  53. 
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Figure  53.  Reaction  Wheel  Effective  Eigenaxis  for  Time-Optimal  Maneuver  -  Numerical 

Approach  with  [1.333,1.569,1.807]  Tmax  0’s=25.4  kg-m2) 

The  sample  data  points  in  Figure  53  are  seen  to  lie  below  the  1.333rmax  curve, 
demonstrating  that  the  spherical  system  devised  during  the  analytical  approach  made  an 
achievable  configuration.  Since  the  actual  optimal  runs  were  faster  than  the  curve  fit,  it 
shows  that  the  spherical  approach  using  the  reduced  value  of  xBeff  tends  to  give  a  slightly 

conservative  estimate  of  the  true  performance  of  time-optimal  solutions.  This 
demonstrates  that  the  maneuver  times  obtained  by  assuming  a  trapezoidal  input  via  an 
effective  eigenaxis  slew  is  achievable  via  optimal  control  all  the  time.  This  allows  the 
spherical  shape  with  momentum/torque  radii  equal  to  1.805/1.333  of  the  individual  wheel 
momentum/torque  maximum  to  be  used  in  later  in  this  thesis  to  facilitate  comparison. 
This  35%  increase  in  momentum  over  47t  steradians  is  significant. 
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°P’m 


1.805  hn 
1.333  h 


=  1.354 


(176) 


In  Chapter  VII,  it  will  additionally  be  shown  how  the  optimal  maneuver 
maximizes  the  use  of  the  momentum/torque  envelope  to  achieve  better-than-eigenaxis 
performance. 


J.  EFFECTIVE  EIGENAXIS  OF  PSEUDO-INVERSE  LIMITED  TIME- 
OPTIMAL  MANEUVER 

If  adding  optimal  torque  signals  directly  to  the  individual  wheels  is  out  of  the 
question,  then  the  system  can  still  be  optimized  for  trajectory  shaping  by  using  a  slight 
variation  from  the  time-optimal  logic  used  in  Figures  34-36.  The  system  must  still 
employ  the  pseudo-inverse  within  the  reaction  wheel  dynamics,  but  performance  can  be 
improved  within  this  space.  Figure  54  illustrates  the  pseudo-inverse  limited  optimal 
shaped  trajectory  being  used  in  the  spacecraft  model. 


Ortxtal  Motion 
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Quaternion  Quaternion 
Trajectory  Error 
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Body  Rate 
Trajectory 


Figure  54.  Pseudo-Inverse  Limited  Time-Optimal  Spacecraft  Model 


The  accessible  gains  are  highlighted  in  Figure  55.  Granted,  the  blue  regions  are 
accessible  without  the  use  of  optimal  control,  but  not  about  every  eigenaxis.  This  section 
demonstrates  how  optimal  control  can  be  exploited  for  equal  gain  over  47t  steradians. 
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Figure  55.  Pseudo-Inverse  Limited  Time-Optimal  Accessible  Regions 


In  the  same  fashion  as  analytically  extracted  for  the  time-optimal  in  Sections  H 
and  I,  the  analytical  and  numerical  approach  can  be  completed  for  the  pseudo-inverse 
limited  time-optimal  maneuvers.  It  can  be  shown  that  the  pseudo-inverse  limited  time- 
optimal  space  can  be  achieved  by  using  maximum  magnitudes  of  1.577/1.333  of  the 
individual  wheel  momentum/torque  maximums. 


h  =1  577  h 

pseopt— sphere  *  max 

^ pseopt- sphere  ^  »333  ^max 


(177) 
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Figure  56.  Pseudo-Inverse  Limited  Time-Optimal  Momentum  Sphere 


By  shaping  the  trajectory  using  pseudo-optimal  inputs  rather  than  the  standard 
trapezoidal  input,  18%  increases  in  speed  and  acceleration  can  be  achieved. 


hB 

pseoptm 


B 


1 .577  hn 
1.333/7 


1.183 


(178) 


K.  SUMMARY 

This  chapter  was  used  to  develop  and  analyze  the  time-optimal  slews.  First,  the 
time-optimal  slew  for  the  eigenaxis-restricted  case  was  demonstrated  as  a  bang-bang,  or 
as  a  bang-coast-bang  maneuver  when  momentum-limited.  Next,  the  eigenaxis  was 
abandoned  and  the  time-optimal  reorientation  was  developed  and  implemented  in  DIDO. 
Slew  time  improved  over  the  trapezoidal  input  by  up  to  22%  for  120°  slews.  This  optimal 
control  approach  allowed  the  limitations  of  the  pseudo-inverse  to  be  avoided,  increasing 
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the  available  momentum  and  torque  envelope.  Next,  the  effective  eigenaxis  concept  was 
introduced,  which  allowed  a  comparison  between  the  eigenaxis  reorientation  and  the 
time-optimal  slew.  The  effective  eigenaxis  was  further  developed  to  show  how,  on 
average,  the  available  momentum  and  torque  for  the  time-optimal  case  are  1 .805  hmm  and 
1.333  Tmax  about  the  effective  eigenaxis.  For  the  pseudo-inverse-limited  time-optimal 
case,  they  are  1.577  hmax  and  1.333  Tmax ,  respectively. 
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VI.  CONTROL  MOMENT  GYROSCOPES 


Like  reaction  wheels,  CMGs  also  have  their  physical  limitations.  They  have 
maximum  angular  rates  and  accelerations  that  must  be  modeled.  This  ensures  that  only 
the  maximum  individual  CMG  torques  can  be  applied,  regardless  of  demanded  torque.  In 
order  to  appropriately  apply  limitations,  a  model  of  the  CMGs  must  be  created.  This 
chapter  is  a  focus  on  single  gimbal  control  moment  gyroscopes  (SGCMGs).  Dual  gimbal 
and  variable  speed  CMGs  will  not  be  discussed.  The  following  introduction  is  simplistic 
regarding  SGCMGs,  but  representative  of  the  key  interactions.  The  rest  of  this  chapter 
will  be  used  to  derive  these  terms  and  show  the  interaction  with  the  spacecraft  body. 

At  the  center  of  a  CMG  is  a  large  momentum  wheel,  driven  by  its  own  wheel 
motor  at  a  nominally  constant  speed.  The  basic  concept  of  a  CMG  is  for  a  wheel  to  spin 
freely  within  a  secondary  gimbal  frame.  The  gimbal  frame  is  rotated  within  the  body 
frame  by  a  gimbal  motor.  The  change  in  the  direction  of  the  wheel  angular  momentum  as 
a  result  of  gimbal  motion  causes  a  large  torque  orthogonal  to  the  gimbal  and  angular 

momentum  axes,  i.e.  f  =  <5  x  h  .  Since  a  relatively  small  gimbal  torque  input  can  become 
a  much  larger  output  torque,  this  concept  is  known  as  “torque  multiplication”  [27,  p.  1]  or 
“torque  amplification”  [28,  p.  160]. 

A  basic  model  of  a  CMG  is  shown  in  Figure  57.  At  the  center  of  the  CMG  is  a 
momentum  wheel  shown  in  the  wheel  frame  W .  The  wx  and  w2  axes  are  fixed  to  the 
wheel  and  rotate  in  the  gimbal  frame  as  a  function  of  the  constant  wheel  rate  £2  and  time. 
It  is  important  to  note  that  the  spin  axis  of  the  wheel  frame  is  aligned  with  the  gimbal 
frame,  such  that 

G(Ow  =Qw3=Qg,  .  (179) 

The  wheel  motor  is  attached  to  the  gimbal  frame  and  drives  rotation  of  the  wheel. 
The  wheel  is  spun  at  a  constant  rate;  therefore,  the  simple  rotation  of  the  wheel  alone 
does  not  cause  a  torque  transfer  between  the  wheel  and  gimbal  frames. 

<  =  /w£&  =  0  (180) 
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Figure  57.  CMG  Schematic 


Figure  57  shows  the  small  gimbal  momentum  and  torques  as  small  red  and  green 
arrows,  respectively.  The  large  wheel  momentum  and  large  resultant  orthogonal  torque 
are  shown  in  with  larger  arrows.  All  of  these  momentums  and  torques  can  be  labeled 
relative  to  the  gimbal  frame.  For  simplicity,  gyroscopic  coupling  terms  are  not  shown  on 
the  diagram. 


K  =  Ksi 

—  ^wSl 

K  = 

Ts  =  zg8  3 


(181) 

(182) 


Consider  the  spacecraft  body  frame  aligned  with  the  gimbal  frame  shown  in 
Figure  57.  Since  the  gimbal  frame  is  allowed  to  rotate  within  the  body  frame,  the 
following  set  of  time  derivatives  with  respect  to  the  body  frame  are  important  to  note. 
They  will  be  used  in  the  chapter’s  derivations. 
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A.  ORIENTATION 

The  description  of  the  orientation  of  a  CMG  will  be  carried  out  much  like  that  of 
the  orientation  of  reaction  wheels  in  Chapter  IV.  To  avoid  redundancy,  some  portions  of 
the  explanation  will  be  condensed. 

The  fundamental  formula  for  the  angular  momentum  gimbal  to  body  change  of 
basis  is  given  as 

hB  =  [BR°‘]h,G'  (184) 


The  dyadic  form  of  this  projection  onto  the  body  frame  is 


h 


k*gi2 

• 

_ 1 

Gi 

4*<?, 

4  •  £«■ 

hG 

(185) 

4  •  k,. 

The  orientation  will  be  placed  in  the  B  frame,  converting  to  the  G,  frame  via 
g'Rb  ,  and  then  using  the  transpose  BRG‘  to  convert  back.  In  the  example  shown  in  Figure 

/V  A 

58a,  the  MED  is  placed  directly  on  the  B  xyz  origin,  rotating  about  the  bx  axis  and 
gimbaling  around  the  4  axis  per  Figure  58b.  This  produces  an  angular  momentum  hB 
directly  in  the  bx  direction  for  the  wheel. 
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a)  Initial  aligned 
with  body  axis 


b)  After  3  (a)  biAg3 
rotation 


b,  8 , 


b  ->  8  2 


c)  After  2  (P) 
rotation 


d)  After  3  (6) 
rotation 


Figure  58.  CMG  Transformation 

The  first  two  orientation  transformations  are  exactly  like  those  for  the  reaction 
wheel.  A  fundamental  3  rotation  about  the  share  b3  /  g3  axis  is  completed. 


/?3  («,.)  = 


(  r\  \ 

cos  a i  sina,.  0 
-sin  a,,  cos  a,.  0 
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0  1 


(186) 


A  2  rotation  is  now  completed  about  the  g2  axis  (Figure  58c)  to  create  the 
orientation  skew  angle  /l . 


R;(W  = 


cos/3 

0 

-sin/3 
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sin/3 
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(187) 


The  distinct  difference  between  reaction  wheels  and  CMGs  occurs  in  the  third, 
gimbal  rotation.  This  is  where  CMGs  derive  most  of  their  advantage.  The  final  rotation  is 
a  3  rotation  by  angle  8  about  the  new  g3  axis  (Figure  58d).  For  reaction  wheels,  it  is 
assumed  that  8  =  0 . 


cos  (5, 
-sin  (5; 
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sin  <5,  0 
cos<5;  0 


(188) 


The  previous  series  of  rotations  can  now  be  combined  into  a  single  transformation 
converting  a  vector  from  the  body  frame  to  the  gimbal  frame. 


G‘RB  =  R3(8i)R2(P)R3(ai ) 
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It  follows  from  Equation  (45)  that  BRa‘  is  now  the  transpose  of  G‘RB ,  which  is 
equivalent  to  the  reverse  sequence  via  the  negative  angle  of  each  rotation. 


-a,)R2(-p)R3(-8i)  =  [G‘RBJ 

cajc[3cSi  -  sars8i  - sajc8i  -  ca  cj3s8j 

-i  G 

cais/3 

cajs8j  +  salcpc8i 

caic8i  -  sajcfis8l 

SCXjSfi 

-spc8i 

sps8i 

cP 

(190) 


As  an  example,  place  the  initial  CMG  at  angle  a  =  0° .  The  resultant 
transformation  from  the  G  frame  to  the  B  frame  is  shown  in  Equation  (191). 
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An  array  of  four  CMGs  is  typically  arranged  in  a  square  pyramid  (similar  to 
reaction  wheels)  in  which  each  of  the  CMGs  are  placed  90°  apart.  The  arrangement  used 
here  will  be:  a,  =  0°,a2  =  90°, a3  =  180°, a4  =  270°.  Figure  59  demonstrates  this 
placement.  This  orientation  matches  that  shown  in  [29,  p.  3]  given  that  the  author’s  initial 
gimbal  angles  were  each  set  at  90°  rather  than  <5,  =  0°  as  assumed  in  this  thesis. 


Based  on  the  wheel  axis  description  mentioned  above,  the  gimbal  and  wheel 
angular  momenta  in  the  gimbal  frames  can  be  written  as 
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Transforming  angular  momentum  from  the  G  frame  to  the  B  frame  using  (184) 
and  a  =  0°  yields  Equation  (193).  Note  that  the  gimbal  momentum  in  the  body  frame, 
Equation  (193),  is  equal  to  the  reaction  wheel  momentum  in  the  body  frame,  Equation 
(53).  This  is  because  the  orientation  of  the  gimbal  rotation  axis  in  is  equal  to  that  of  the 
reaction  wheel  spin  axis,  and  both  were  subject  to  the  same  two  initial  rotations. 
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or  in  general  form: 


B 

- 

and  h'l  = 

K 

ti-k 

k-h 

ti'k 

(193) 


(194) 


Recall  the  total  angular  momentum  of  the  system  is  the  sum  of  each  of  the  MEDs 
in  the  body  frame.  Unlike  the  single  momentum  set  for  reaction  wheels  in  Equation  (69), 
CMGs  derive  their  momentum  from  two  different  transfonnation  sets.  This  allows  the 
gimbal  momentum  and  gimbal  torques  to  be  translated  separately  to  through  each  of  the  n 
CMGs  to  the  body  frame. 
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B 

and  YjK  = 
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(197) 


(198) 


Now,  Equations  (195)  and  (196)  can  be  rewritten  as 


K = 1" 
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K  = 

4, 

and  4  = 

K 
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(200) 


(201) 


For  the  square  pyramid  mentioned,  the  orientation  matrices  is  shown  in  Equation 
(60)  for  CMGs. 
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i  G 


(202) 


B.  MODELING  CMG  ORIENTATION  MATRICES 

These  CMG  orientation  matrices  can  be  created  within  Simulink.  Figure  60 
demonstrates  how  Equation  (190)  is  implemented.  Since  a,  and  /3  are  fixed  for  the 
CMG,  their  respective  rotation  matrices  will  be  held  constant.  The  CMG  transformation 
matrix  BRG‘  is  only  state  dependent  upon  <5, .  The  column  selectors  simply  extract  the 
first,  second,  and  third  columns  individually. 


ith  CMG  Orientation  Matrix  Column  BzGi 


Figure  60.  CMG  ith  Orientation  Matrix  Column  Model 

Additionally,  Figure  61  illustrates  how  the  R  ,  and  R2  matrices  can  be  created  for 
a  generic  angle. 
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Figure  6 1 .  Rotation  Matrix  Model 


Figure  62  fully  creates  all  three  orientation  matrices.  BZG‘  from  Figure  60  is 
represented  four  times,  once  for  each  CMG.  The  individual  columns  are  independently 
combined  to  create  the  appropriate  3x4  orientation  matrices  BZG  . 
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Orientation  Matrices 


BZG1 


Figure  62.  CMG  Orientation  Matrix  Model 


C.  DIFFERENTIAL  EQUATIONS  OF  MOTION 

The  momentum  of  a  CMG  can  be  defined  as  the  sum  of  the  wheel  angular 
momentum  and  the  gimbal  angular  acceleration. 


HCMG=HNw+HNgi 


Each  the  wheel  and  gimbal  can  be  independently  modeled. 


T J  T  A r^W: 

Hw  =  JW;  CO, 


T  ( .  B^G,  .  G,  W:  \ 

JWi[  ®  ®  j 


/N  _  _  B  .  j  tS ,  _  U:  .  j  Lr-  _  _ 

w,  ®  +JWi  +JWi  °>i 


G: 


(204) 


(205) 
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(206) 


1.  Gimbal  Momentum  and  Torque 

a.  CMG  Gimbal  Momentum  -  Single  Gimbal 
The  gimbal  momentum  is  shown  first  in  matrix  form. 

hin  00  0 

K;=Jg>?‘=  0  ^  o  0 

0  0  4,  J  [ 

Shown  in  dyadic  fonn: 

hg,‘  —  l<?!  l<?il  +  j gi22 8 i 28 i 2  jgi338i3§i3  )  *  ) 

_  Jgi^iSid  (Si3  *  8i3  ) 

_  jgi '33  3 

Change  of  basis  to  the  body  frame: 

K  =  jgijiSn  •  (kk  +  kk  +  kk ) 

=  jgiji  (t,3  •  k  )k + jgji  (g/3  • k)k  +  jgiji  (g,3  •  k)k 

The  resultant  rotation  matrix  is  shown  as 


(208) 


(209) 


(210) 
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b.  CMG  Gimbal  Momentum  -  Gimbal  Array 

Similar  to  reaction  wheels,  the  resultant  orientation  matrix  can  be  created 
using  the  components  from  the  3ld  column  of  the  gimbal  frame  rotation  matrix. 


where 


(213) 


(214) 
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c.  CMG  Gimbal  Torque 

The  torque  of  the  gimbal  can  now  be  determined  with  respect  to  the  body 

frame.  Since  this  is  simply  the  gimbal  axis,  the  orientation  is  fixed  and  Z  =  0  .  The  inertia 
is  assumed  not  to  change  relative  to  the  gimbal  axis. 


A) 


dt 


dt 


(215) 


*  ^ 


WB  BrjG 


(216) 


In  dyadic  form,  a  single  CMG’s  gimbal  torque  can  be  found  using  the 
derivative  of  Equation  (210). 


BA 

dt 


(217) 


The  inertia  of  the  gimbal  is  not  changing  with  respect  to  the  gimbal  frame. 
Also,  the  orientation  of  the  3  axis  of  gimbal  frame  is  not  moving  with  respect  to  the  body 
frame  as  determined  in  Equation  (183). 

Jgm  =  0  and  gi3  =  0.  (218) 


Equation  (2 1 7)  now  becomes 


B_d_ 

dt 


(219) 


(220) 
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Wheel  Momentum  and  Torque 

a.  CMG  Wheel  Momentum  -  Single  Wheel 

The  CMG  wheel  momentum  is  shown  first  in  matrix  form. 
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An  equivalent  derivation  is  shown  in  dyadic  fonn. 

K‘  =(jWlJil8ii  +jWiJJi2  +  44, ) 

=  j^P-ih  (&,  •  8h )  +  jWiJi8h  ( 8 /,  *  8h ) 

_  Jwnpi8il  +  Jwi33^i8i3 

Change  of  basis  to  the  body  frame: 

K,  =  (jWnPi8h + jWi3A8i)'{kk +44+44) 

=  jWnPi{h  •bl)bl  +  jWiipi(gh  •b2)b2+jWnPi(gil  *4)4 
+  jWiJi(8i3'bl)bl+jWlJi(8i3  *4)4  +  uji(gi3  *4)4 


(221) 


(222) 


(223) 


The  resultant  rotation  matrix  is  shown  as 
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(225) 


b.  CMG  Wheel  Momentum  -  Wheel  Array 

The  resultant  orientation  matrix  can  be  created  using  the  components  from 
the  1st  and  3ld  columns  of  the  gimbal  frame  rotation  matrix. 
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iB  B  rjG  j  oG  ,  B  rjG  j  G 

K  =  Z3  J  d  +  Z,  JJ2 


(227) 


c.  CMG  Wheel  Torque  -  Wheel  Array 

The  torque  of  the  gimbal  can  now  be  determined  with  respect  to  the  body 
frame.  Since  this  is  simply  the  gimbal  axis,  the  orientation  is  fixed  and  BZ^  =  0 .  The 
inertia  is  assumed  not  to  change  relative  to  the  gimbal  axis.  Also,  the  wheel  is  kept  at  a 
nominally  constant  speed. 


T 


B 

IV 


(228) 


(229) 


Using  partial  differentiation  of  the  second  term  within  Equation  (229),  the 
following  is  obtained. 


Z?JWQ?  =  ~(‘Z?)JA 


BrjG 


dS  dt 
dd  d 
dt  dd 


Therefore,  BZ°  =  —(BZ°)S 


where 


can  be  interpreted  as  a  CMG  Jacobian  matrix. 


cajcj5cdi  -  sajsdj 

-sajcdj  -  ca  c[isdj 

ca^P 

cajsdj  +  saic/5cdi 

caic8!  -  saicf3s8i 

satsj3 

-sficd; 

sfis8i 

cA 

(230) 


(231) 


(232) 
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By  inspection  of  the  transformation  matrix,  the  partial  derivative  of  the 
first  column  is  equal  to  the  second  column. 


—(bZig')=bZ2c 

v  1  /  2 


(233) 


Similarly,  this  is  captured  in  the  dynamic  orientation  matrix  [6,  pp.  464- 
467],  where  A  is  introduced  as  a  steering  logic  matrix. 


hB  =  A8h?  where  A  ■ 


(234) 


Upon  close  inspection,  A  matrix  is  equivalent  to  the  second  column  of  the 


Z  matrix  utilized  within  this  thesis. 


A  =  Z: 


(235) 


An  alternative  proof  follows  using  dyadics. 

~\B 

Vi-, 


(236) 


V&, 


Using  g,  =  8g2  from  Equation  (183),  Equation  (236)  becomes 


%%Pi  =  b2-gh  urn 

k  •  L 


(237) 


In  summary,  all  variations  are  equivalent 

=  XJ,$S  =  AJ.OfS  =  ^(’Z?)j,Qcg 

“z:;  =  =  .4<i'  =  —  ( ”z.°  is 

1  2  d8[  1  ’ 


(238) 
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The  total  wheel  torque  from  Equation  (229)  now  becomes 


Tb  =  BZGJ  8G  +  BZGJ  Qg8 

W  3  **  W  Q  2  **  w  g 


(239) 


Assuming  like-sized  wheels  such  that  hw  =  jWnp.x  =  A21  02  =  jWni  Qn 
becomes  a  scalar  term,  the  equation  can  be  simplified  into 

<  =  jJZG8G+jwaGBZG8  (240) 


Equation  (240)  shows  the  final  torque  equation  for  the  wheels  of  a  CMG 
array.  The  equation  is  dominated  by  the  second  term  [29].  The  first  term  is  dependent 
upon  gimbal  accelerations,  which  are  generally  slow  [27].  The  second  term  is  dependent 
upon  a  slow  gimbal  angular  velocity,  but  it  is  amplified  significantly  by  the  very  high 
wheel  rate.  This  is  the  “torque  amplification”  that  the  literature  suggests  [27,  p.  1] 

[28,  p.  160]. 

In  dyadic  fonn,  a  single  CMG’s  wheel  torque  can  be  found  using  the  time 
derivative  of  Equation  (223). 


dt 


(K)=  Luai(8n'b  k)k  +iw,nQ/(£n  *k)k  +Lhfli{girk)k 
+  A,,  A  (in  *4)4+  A;,  A'  (in  *  k  )k  +  A;nA  (in  *  k  )k 

+  A/,  A-  (in  * k)k  +  A,, A-  (in  * k ) k  +  A/,  A/  (in  * k )k 

+ a  iji  ( I  /  3 -k)k  +  Lji  ( I  /  3  •  k  )k + Jwi3A  ( 1/3  •  k  )k 

+  jmji  ( I  /  3  •k)k  +  a, A  (l,-3  •  k  )k  +  jmJi  (l/3  •  k  )k 
+  JwhA  ( I  /  3  *  k  )k  +  JmJi  ( 1/3  *  k  )k  +  Jwj,  ( I  /  3  •  k )  k 


(241) 


The  inertia  of  the  wheel  is  not  changing  with  respect  to  the  gimbal  frame. 
Also,  the  orientation  of  the  3  axis  of  gimbal  frame  is  not  moving  with  respect  to  the  body 
frame  as  determined  in  Equation  (183).  The  wheel  speed  is  assumed  constant. 


A,,  =  0  »  An,  =  0  ,  In  =  <5,1,2 ,  1,3  =  0  ,  A  =  0 


(242) 


Equation  (241)  now  becomes 
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(243) 


Bd 


dt  (K )  =  UP  A  (gi2  •  4)4 + UP  A  (gi2  •  4)4 + j^pA  (&2  •  4)4 
+  jwiJi  (i,3  *4)4+  Lji  (gi3  *  4  )4  +  jwiJ,  (<4  •  4  )4 
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The  first  term  again  validates  the  claim  made  in  Equation  (237). 

3.  Total  Momentum  and  Torque 

The  combined  CMG  gimbal  and  wheel  momentum  and  torque  are  now 


kmc=k+k 


=  ( “z,  J,Sr; ) + ( "z°  +  "zf  ) 

=(A,+^)"z3»5»+i/z^ 


—  5  i  _  B 

^ CMG  ~  ^  g 
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V  3  Jg33  g  )  V  3  Jw?,ig  2  Jwu  g” ) 

=  BZG(j  +j  )8G  +  BZG(j  Qc)«5 

3  \  ^  833  Jw33  /  8  2  \*'wll  8  / 

=  (/  +/  )BZG8G  +  (j  Q G)BZG8 

\Jg} 3  Jw33  J  3  u{  g  ) 

By  Euler’s  transport  theorem,  the  inertial  derivative  is 

"dK°) 


-+NcoBxhB 


dt 


dt 

i'B  ,  N 


CMC 


"'CMG  ~  " CMG  ^  n cmG 


(244) 


(245) 


(246) 


(247) 


Equation  (247)  can  be  rewritten  by  substituting  the  angular  momentum  terms 
hN=  (j  +j  )BZG8G+(j  QG)BZG8 

8  #33  *^33  /  3  g  y-'H’n  g  J  2 


+  NcoG  x 


((j  +j  ) BZG8G  +  j  BZG8G)' 

\\''  g33  Jw33  )  3  g  Jwn  3  g  I 
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(248) 


Like  reaction  wheels,  CMGs  are  subject  to  the  same  pseudo-inverse  law.  Much 

like  the  pseudo-inverse  of  the  reaction  wheel  orientation  matrix,  BZG  (or  A)  can  be 

Z+ 

inverted  using  the  same  method.  Here  will  be  used  to  map  the  three  body  axes  back  to 
n  CMGs. 


-xN  =  { j  +j  )GzG8G+(j  aG)BzG8 

req  \J gyi  ■'W33  )  3  vg  \Jwt,  g)  2” 


+  N(0G  x 


((jg,3  +  Jw33)GZ3^g  +  j. 


B 


Z  8g 


(249) 


Although  Equation  (249)  shows  the  full  dynamics,  the  gyroscopic  coupling  of  the 
body  motion  with  the  gimbal  rate  can  be  assumed  as  very  small  [29].  Additionally  the 
gimbal  acceleration  can  be  assumed  small  [29].  Therefore,  Equation  (249)  can  be  reduced 
to  Equation  (250)  and  rearranged  per  Equation  (25 1)  using  the  pseudo-inverse. 

-rN  =  {j  Qg)bZg8 

''req  ^tvn  g  )  2  v 

=  BZGhG8  (250) 

B  rjG  G 

—  Z0  T 

2  vv 

=(ft°r'['z?T(-o  <25i> 

=(*.")"  K) 


This  is  the  same  steering  logic  proposed  in  [6,  Algorithm  7.147].  Therefore,  the 
assumptions  above  were  taken  as  well. 


D.  MODELING  CMGS 

In  order  to  apply  the  physical  actuator  limitations  of  the  CMG  as  part  of  the 
simulation  model,  the  torque  required  of  the  assembly  must  be  broken  down,  individual 
gimbal  and  wheel  limit  saturations  applied,  and  then  the  torque  of  the  CMG  assembly 
rebuilt.  Within  Figure  63,  (251)  is  used  to  break  down  the  torque  required  and  allocate  it 
across  the  four  wheels.  The  orientation  matrix  Z  block  consists  of  Figures  60-62. 
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CMG  Steering  Law 

(saturation  occurs  in  red) 


This  CMG  model  can  be  further  expanded  to  include  null  motion  in  the  presence 
of  near-singular  conditions.  Additionally,  the  components  of  Equation  (249)  can  be 
modeled,  but  Equation  (25 1)  is  modeled  here  for  simplicity. 

E.  CMG  MOMENTUM  SPACE 

The  boundary  created  by  the  maximum  momentum  in  any  given  direction  for  a 
particular  CMG  configuration  is  much  more  difficult  to  define  than  that  of  a  reaction 
wheel  configuration.  The  momentum  space  is  no  longer  a  three-dimensional  geometric 
envelope  with  all  interior  points  possible.  The  physical  limitation  is  still  the  vector  sum  of 
the  maximum  capacitates  of  the  individual  CMGs,  where  momentum  saturation  occurs. 
Unlike  reaction  wheels,  the  surface  is  not  a  polyhedron.  Instead,  the  CMG  momentum 
envelope  is  contains  singular  cavities  that  pass  through  the  center  from  one  edge  of  the 
momentum  space  to  the  other. 

Like  the  reaction  wheel  case,  maximum  momentum  space  for  an  evenly  spaced 
(Act  =  90)  four-CMG  array  is  when  /3  =  54.73°  [30,  p.  866].  The  resulting  surface  is 
shown  in  Figure  64.  Since  the  gimbal  axis  g3  and  the  wheel  angular  momentum  g,  are 

always  perpendicular,  a  reduced  amount  of  wheel  momentum  is  available  on  the  surface 
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of  the  momentum  envelope,  which  is  set  within  the  body  frame 

(&,  x  kmgA  =  1  or  8i,  * k,rzA  =  Si2  * h^A  =  O) .  These  regions  are  identified  by  white 

cavities  in  Figure  64,  in  which  2"  cavities  exist  [28,  p.  161].  The  identification  of 
singular  surfaces  will  be  discussed  in  Section  E. 


Figure  64.  Momentum  Space  of  a  4  CMG  System  Given  (3=54.73° 


The  shape  within  Figure  64  is  unique  to  this  four-CMG  system.  The  shape  alone 
reveals  configurations  of  maximum  momentum.  The  momentum  peaks  help  identify  the 
best  possible  axes  to  rotate  about.  These  maximums  occur  at  =  3.266 hmax ,  revealing 
that  when  using  CMGs  in  comparison  with  similarly  sized  reaction  wheels,  40%  more 
momentum  is  achievable  about  certain  axes. 
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The  concept  of  placing  a  momentum  sphere  directly  inside  the  minimum  radius  of 
this  shape  does  not  make  sense  because  the  shape  is  concave.  This  leads  the  discussion  to 
CMG  singularities. 


F.  CMG  SINGULARITIES 

Within  a  singular  gimbal  angle  configuration,  it  is  not  possible  to  produce  a 
torque  in  certain  directions.  Therefore,  the  required  torque  may  or  may  not  be  produced. 
Singular  states  occur  when  the  rank  of  BZ2  is  not  full  (less  than  n  CMGs).  This  condition 
can  occur  because  BZ2  is  state  dependent. 

b  V  1G- 

-saic8i  -  ca  c/fsfi 

BZ2<  (<5, )  =  caic8i  -  saicfis8I  (253) 

sj5s8i 

If  a  configuration  of  delta  leads  two  CMGs  to  yield  the  same  columns  within  their 
respective  Z  matrix,  then  a  singular  state  occurs.  The  entire  array  for  the  [0,90,180,270] 
arrangement  is  shown  in  Equation  (254). 

b  r  ~\G 


-c(5s8t 

~c82 

c(5s83 

-c84 

c8l 

-cj5s8n 

~c8. 

c[3s84 

(254) 

spsS, 

s(3sS2 

s/5s83 

s[3s84 

For  example,  when  5,  =  0°  and  S3  =  180° ,  two  columns  are  the  same.  Therefore, 
rank(  BZf)  =  2,  which  is  not  full.  Therefore,  a  singular  condition  results  during  the 
pseudo-inverse. 

B  - 

0  ~cS2  0 

bZ2g(0°,<52,180°,<S4)=  l  -cps82  l 

0  s(3s82  0 
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This  is  just  one  of  many  different  configurations  can  reduce  the  rank  and  cause  a 
singular  condition. 


1.  Solving  Singularity  Conditions 

A  detailed  procedure  and  set  of  equations  for  determining  the  singularity 
conditions  and  associated  singular  momentum  vectors  is  given  in  [6,  pp.  676-678],  CMG 
singularity  occurs  when  the  magnitude  of  the  dot  product  of  the  wheel  momentum  vector 
and  the  spacecraft  rotation  axis  is  one:  g.  •hfmgA  =  hf  •hBmgA  =  ±1 .  In  this  thesis,  only  the 

results  of  the  four  CMG  case  will  be  analyzed.  The  singularity  condition  for  each  CMG 
can  be  written  as 

k  x  hq„gA )  x  k  _  +  (g,-3  x  hBcmgA )  x  l3 
h?  ‘KngA  gh  x  h"mgA 


(256) 


e- 


hB  =  ±- 


The  vectors  et  contain  the  cross  product  of  the  individual  CMG  gimbal  vectors 

i'  B 

and  the  total  body  momentum  unit  vector  hcmgA  .  Both  positive  and  negative  variations  are 
accounted  for  within  Equation  (257)  by  the  sign  function  of  Equation  (258). 


G  =  ^ 1  “  ( ' SP  +  CP  ) 

e2  =  ±J\-(sf3  hBcmgAy  +  cfi  hBmgA^ ) 


(257) 


e3  =  (~SP  hcmgA,  +  CP  hBcmgAz  ) 


e,  =  +J1-  - 


-sPh*mgAy  +  cph*mgA) 


ei  =  sign(ei)  =  ± 1  (258) 

The  total  CMG  array  momentum  in  the  body  frame  is  created  using  Equation 


(259). 


Ca=  5X  =  Z 


_^(g;3xfcgA)xg/3 


(259) 


n—l 


n= 1 
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Individually  accounting  for  each  body  axis  separately,  Equation  (259)  can  be 
written  as  Equation  (260). 


_cp(cph^mgA^sf5hBmgA)  ^  h»mgK  ^  cp(cphBmgA-  sph*mgA  )  ^  h 


K  = 


hB  ,  c 

hB  =  ^nS^  + 


B 

cmgAx 


P[sphBcmgArcPKngA)  |  hBcmgAy  |  cp(sph*mgA  +cphBmgAy) 


hB  = 


P{*P  hBmgAz  ~  cp  hBcmgAx )  |  hBcmgA  -  c/3  hBmgAy ) 


(260) 


+  ■ 


P ( sfi hBcmgAz  +  cp hBcmgAs )  ^  .v/j(.v/j hBcmgA;  +  c/3 hBmgAy ) 


Given  four  CMGs,  there  are  16  (2")  combinations  of  e(  that  produce  different 

singular  conditions.  When  e(.  has  the  same  sign  for  all  CMGs  |e  =  ±[  1  1  1  1  ]  j , 

the  outer  singularity  surface  shown  in  Figure  64  is  produced.  Other  combinations  result  in 
interior  singularities:  states  that  are  contained  within  the  volume  of  the  CMG  momentum 
envelope. 

Interior  singularities  mean  that  at  some  non-maximum  momentum  (within  the  set 
defined  by  the  geometry  of  Figure  64),  the  required  torque  magnitude  or  direction  cannot 
be  created  exactly,  creating  an  off-eigenaxis  rotation.  Two  types  of  internal  singularities 
exist:  hyperbolic  and  elliptic  singularities.  The  names  derive  from  the  shape  of  the  gimbal 
angle  motion  in  the  neighborhood  of  a  singular  point.  These  are  discussed  in  detail  within 
[3,  pp.  23-30]. 


2.  Hyperbolic  Singularities 

Hyperbolic  singularities  are  the  milder  of  the  two  internal  singularities, 
considered  passible  in  [3,  pp.  23-30].  This  is  the  case  when  the  CMG  array  cannot 
exactly  produce  the  required  torque.  A  specific  gimbal  signal  can  be  added  to  the  steering 
logic  whenever  a  singular  state  is  detected.  This  null  motion  produces  no  net  torque  on 
the  system,  but  assists  in  passing  through  hyperbolic  singularities.  They  occur  when 
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exactly  two  of  the  et  have  mismatched  signs.  The  six  combinations  are  shown  in  the 
following  three  cases. 


e,-=± 

1 

1 

-1 

»  ei  =  ± 

1 

-1 

1 

-  ei=± 

1 

-1 

-1 

-1 

-1 

1 

(261) 


The  hyperbolic  singularities  make  up  the  very  center  of  the  shape,  extending  from 
the  origin  (zero  momentum)  to  a  maximum  momentum  magnitude  of  1.897  at  its  largest. 

A  graphical  depiction  of  case  2  |e  =  ±[  1  -1  1  -1  ]  j  is  shown  in  Figure  65.  It  is 

noteworthy  that  it  is  symmetric  about  all  three  body  axes. 


-05 


Hyperbolic 
Singularity 
(case  2) 


e.  =  ± 


Figure  65.  Hyperbolic  Singularity  Surface  for  a,=±[  I  -1  1  -1]T 
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Figure  66  demonstrates  all  three  cases  plotted  together.  Each  of  the  cases 
extended  from  the  origin  to  1 .897  hmax . 


Figure  66.  Flyperbolic  Singularities 


3.  Elliptic  Singularities 

Elliptic  Singularities  occur  when  the  torques  allowable  by  the  physical  layout  of 
the  CMGs  are  perpendicular  to  the  required  torque  [11,  p.  179].  They  are  impassible  [3]. 
Elliptic  singularities  occur  when  only  one  e(  sign  is  mismatched.  The  eight  possible 
combinations  are  shown. 
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e,=± 


(262) 


-1 

1 

1 

1 

1 

-1 

1 

1 

c  =  + 

£.  =  + 

c  =  + 

1 

?  s  — 

1 

?  s  — 

-1 

1 

1 

1 

1 

-1 

The  elliptic  singularities  extend  across  the  width  of  the  shape,  extending  from  a 
minimum  of  1.0069  to  a  maximum  of  2.998  hmm  at  its  largest.  A  graphical  depiction  of 

case  1  |e  =  ±|^  — 1  1  1  1  ]  j  is  shown  in  Figure  67.  It  is  noteworthy  that  the  singular 

surface  is  oriented  at  the  skew  angle  (/3  =  54.73°)  due  to  the  alignment  of  the  CMG. 


-2  -t  0  1  2 


Figure  67.  Elliptic  Singularity  Surface  (case  1) 


Given  all  cases  at  once,  independent  plots  extend  across  the  momentum  space, 
shown  in  Figure  68. 
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b2 


b3 


Figure  68.  Elliptic  Singularities 


Most  interestingly,  the  eight  interior  horns  shown  in  Figure  68  align  with  the  eight 
hollow  circles  shown  in  the  perimeter  volume  in  Figure  64.  This  union  is  shown  in  Figure 
69. 
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Elliptic  and  Surface  Singularities  Combined 


Figure  69.  Elliptic  and  Surface  Singularities 


G.  CMG  LARGEST  INSCRIBED  SPHERE 

This  thesis  will  not  focus  on  complicated  CMG  singularity  avoidance  algorithms, 
an  ongoing  and  developing  research  area  in  its  own  regard.  In  CMG  systems,  it  may  be 
easiest  to  treat  these  singularities  simply  as  impassible  in  order  to  avoid  exotic 
controllers.  To  create  a  volume  that  is  guaranteed  reachable,  the  minimum  momentum 
magnitude  of  the  elliptic  singularity  region  must  be  completely  avoided.  Since  this 
minimum  occurs  at  the  calculated  values  of  1.0069 /imax  (1 ,0/imax  in  [3]),  this  also  creates 
the  largest  sphere  of  guaranteed  momentum  capability.  Superimposing  a  sphere  into  the 
elliptic  region  creates  the  space  of  guaranteed  mobility,  shown  in  Figure  72. 
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Largest  Guaranteed  Passible  Sphere 


Figure  70.  CMG  Largest  Sphere 


This  reveals  that  the  maximum  guaranteed  momentum  capacity  is  1 .0069  /zmax , 
69%  smaller  than  originally  estimated  in  Equation  (252). 


h 


1.0069  h„ 
3.266  h 


=  0.3083 


(263) 


This  means  the  available  momentum  in  the  CMG  case  is  actually  less  than  that  of 
the  reaction  wheel  cases:  recall  1.633  /i  for  the  full  reaction  wheel  and  1.333  \i  for 

IIldA  1 1 1UX 

the  pseudo-inverse  limited  reaction  wheel. 


H.  SUMMARY 

This  chapter  was  used  to  develop  the  CMG  equations  of  motion  from  the  first 

principle  of  conservation  of  angular  momentum.  The  additional  moving  frame  causes 

CMG  dynamics  to  be  much  more  complex  than  reaction  wheel  dynamics.  The  orientation 

matrix  was  shown  to  be  state-dependent,  which  leads  to  singularities.  Certain  CMG 

arrangements  lead  to  different  types  of  singularities  that  can  cause  controller  failure. 

Elliptic  singularities  are  non-passable  and  hyperbolic  singularities  are  passable. 
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The  singularities  were  individually  mapped  in  three-dimensional  momentum 
space  to  demonstrate  where  they  occur.  The  momentum  envelope  was  combined  with 
both  of  the  interior  hyperbolic  and  elliptic  singularities  to  demonstrate  their  intersection. 
Using  the  simplification  that  elliptic  singularities  must  be  avoided  and  hyperbolic 
singularities  are  passable,  a  sphere  of  accessible  momentum  was  identified.  The  radius  of 
this  sphere  guarantees  the  momentum  of  1.0069  hmm  about  any  axis  of  rotation. 
Therefore,  this  momentum  value  can  be  used  to  maneuver  the  CMG  about  any  eigenaxis. 

Despite  the  lost  momentum  capacity  and  required  complex  ACS  architecture, 
torque  amplification  remains  the  primary  driver  for  CMG  desirability  over  reaction 
wheels.  For  CMGs,  a  small  gimbal  input  results  in  a  large  torque  output.  This  relationship 
is  quantified  in  Chapter  VII. 


133 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


134 


VII.  REACTION  WHEEL  AND  CMG  COMPARISON 


The  previous  chapters  were  used  to  develop  the  necessary  background  and 
equations  of  motion  to  support  a  detailed  discussion  within  this  chapter.  This  chapter  is 
used  to  help  bridge  the  performance  gap  between  reaction  wheels  and  CMGs.  Instead  of 
insisting  that  CMGs  must  be  used  for  highly  agile  spacecraft,  this  chapter  will  introduce 
some  concepts  for  sizing  reaction  wheels  to  attain  the  required  slew  performance.  In 
some  cases  CMGs  will  be  more  desirable,  but  this  chapter  will  illustrate  that  for  some 
satellite  systems,  reaction  wheels  can  potentially  outperform  CMGs. 


A.  MOMENTUM  SPACE  AS  A  VISUALIZATION  TOOL 


In  Chapter  IV,  the  reaction  wheel  momentum  space  was  discussed.  Due  to  the  use 
of  the  pseudo-inverse  for  control  allocation,  the  entire  momentum  and  torque  space  was 
not  available  for  control.  Also,  there  is  no  real  reason  to  stay  on  the  eigenaxis,  despite  the 
fact  the  eigenaxis  maneuvers  have  been  historically  implemented  on  spacecraft.  The  input 
can  be  shaped  in  quaternion  space  to  move  around  the  momentum  space,  searching  for 
the  axes  of  maximum  momentum  and  torque  (vertices  from  Chapter  IV).  This  allows  the 
slew  time  to  be  reduced  and/or  optimized,  regardless  of  the  slew  distance  and  path 
traveled.  Recall  from  Chapter  V  that  for  a  successful  maneuver,  it  is  the  boundary 
conditions  that  must  be  satisfied,  not  the  path  itself. 


Consider  a  120°  eigenaxis  slew  about  e  = 


V2V3  0  2/V3 


T 

.  The  same 


spacecraft/reaction  wheel  configuration  from  Chapter  V,  Sections  F-G  will  be  used. 
A  trapezoidal  input  is  fed  into  the  system,  accounting  for  the  limitations  of  the  pseudo¬ 
inverse.  The  dynamic  response  of  the  reaction  wheel  array  is  used  to  build  Figures  71-72. 
The  vectors  in  the  body  frame  are  nonnalized  against  the  individual  wheel  maximum. 


h“=h“/h. 


iB 

norm 

mB 


r  It 


(264) 


For  a  reaction  wheel  system,  the  shapes  of  the  momentum  and  torque  envelopes 


are  exactly  the  same.  Thus,  normalization  via  (264)  allows  both  momentum  and  torque 
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signals  to  be  compared  on  the  same  scale.  Figure  73  shows  the  scaled  system  response  to 
the  shaped  input.  Only  the  relative  magnitude  is  shown  with  no  regard  for  direction.  The 
system  follows  the  signal  exactly:  bang-coast-bang.  The  only  deviation  is  a  discontinuity 
at  Point  5.  This  occurred  because  the  shaped  trajectory  implemented  cut  the  velocity  off 
when  the  angle  reached  0.999  of  the  final  angle,  causing  a  split  second  of  unpredictable 
acceleration,  and  is  an  artifact  of  the  simulated  model  only. 


Momentum 

Torque 


Time  (sec) 


Figure  7 1 .  Normalized  Momentum  and  Torque  for  a  Trapezoidal  Response 


The  same  maneuver  is  next  mapped  onto  the  momentum  and  torque  space  in  Figure 
72.  Since  the  maneuver  is  an  eigenaxis  maneuver,  the  momentum  travels  from  zero  to 
maximum  momentum  within  the  envelope  defined  by  the  pseudo-inverse  control 
allocation,  yet  the  axis  of  rotation  remains  fixed.  The  numbered  points  reflect  those  from 
Figure  71.  On  the  right  is  the  torque  mapping.  The  locations  of  Points  1,2  are  maximum 
torque  in  the  direction  of  the  eigenaxis.  Point  5  demonstrates  a  strong  torque  reversal  to 
decelerate  the  system. 
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regions  mathematically  unavailable 
by  use  of  pseudo-inverse  ^ 


Figure  72.  Momentum/Torque  Space  During  Eigenaxis  Slew 


The  results  shown  in  Figure  72  are  uninteresting,  and  it  is  not  surprising  that  this 
envelope  mapping  is  not  used  to  help  visualize  conventional  reorientations.  The  same 
maneuver  will  now  be  illustrated  using  the  optimal  control  approach  described  in  Chapter 
V,  Figures  51-53. 
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Figure  73.  Normalized  Momentum  and  Torque  for  an  Optimal  Maneuver 


Now  the  temporal  element  of  the  identified  points  becomes  important.  The  following 

is  a  description  of  the  instantaneous  rotation  axis  as  it  pertains  to  Points  1-5  given  in 

Figures  73  and  74. 

Point  1 .  The  array  pulled  far  off  eigenaxis  at  the  point  of  maximum  torque,  a  maximum 
vertex  at  2.309  Tmax .  This  allowed  the  system  to  gain  momentum  as  quickly  as 
possible. 

Point  2.  The  array  quickly  moved  along  the  edge  of  the  torque  space  to  a  minor  vertex  at 
2-000  Tmax. 

Point  3.  The  array  rapidly  changed  directions  again,  this  time  to  another  maximum  vertex 
at  2.309  Tmax .  Next,  the  torque  decelerated  the  system  as  it  pulled  the  rotation 
angle  close  to  the  eigenaxis. 

Point  4.  The  rotation  angle  wobbled  around  the  eigenaxis,  forming  a  lasso-like  loop  near 
maximum  momentum.  The  torque  remained  close  to  zero,  making  a  few 
adjustments  necessary  to  keep  the  orientation  direction  on  track. 

Point  5.  Point  1  is  reversed  to  decelerate  the  array. 
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Figure  74.  Momentum/Torque  Space  During  Optimal  Maneuver 


The  type  of  action  shown  in  Figure  74  is  common  for  most  optimal  maneuvers.  The 
activity  appears  unintuitive,  but  trends  develop.  Torque  seeks  high  ground  to  raise  the 
momentum  quickly.  Torque  and  momentum  follow  the  edges  and  vertices  as  they  move 
around.  Rarely  does  the  maneuver  stay  very  close  to  the  actual  eigenaxis.  Figures  75  and 
76  represent  two  more  maneuvers  to  further  illustrate  some  of  the  interesting  behavior, 
obtained  as  a  result  of  optimizing  the  slews.  The  same  momentum  limits  are  used,  but  the 
torque  is  doubled  to  0.2  Nm  to  illustrate  different  maneuvers. 

The  optimal  use  of  the  full  momentum  and  torque  envelopes  is  extremely  interesting. 
The  momentum  and  torque  move  like  water  encased  in  a  polyhedron,  flowing  to  high 
surfaces.  The  system  always  seems  to  tend  towards  maximum  torque  in  a  particular 
direction  to  maximize  the  momentum  in  order  to  rapidly  maneuver  the  spacecraft. 
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2.5 


Figure  75.  Optimal  Slew  Example  2,  q0=[0.36  0  0.79  0.5],  qf=[0  0  0  1] 


Figure  76.  Optimal  Slew  Example  3,  qo=[0.23  0.17  0.82  0.5],  qt=[0  0  0  1] 
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B.  MOMENTUM  AND  TORQUE  RELATIONSHIP 


Many  variables  go  into  detennining  the  torque  output  of  reaction  wheel  and  CMG 
systems.  This  section  uses  market  data  to  develop  a  “rule  of  thumb”  relationship  that  can 
be  utilized  in  size  comparisons  between  the  two  very  different  systems. 


1.  Simplified  Dynamics 

In  order  to  size  a  MED  array,  a  reasonable  approach  would  be  to  use  the  largest 
inertia.  In  the  same  way  that  the  minimum  momentum  limits  angular  rate,  the  largest 
inertia  limits  rate.  A  starting  point  for  sizing  wheels  might  begin  without  the 
complications  of  gyroscopic  coupling  per  Equation  (265). 


h 


N 

medA 


T 


N 

medA 


(265) 


In  Chapter  III,  it  was  shown  that  the  spacecraft  can  be  modeled  as  a  simple  double 
integrator.  This  creates  a  linear  relationship  between  spacecraft  rate/acceleration  and 
momentum/torque,  as  both  use  the  same  proportionality  factor  J.  This  allows  a  very  basic 
relationship  between  the  dynamics  of  the  spacecraft  and  the  momentum  exchange  device 
to  be  considered.  As  the  engineer  narrows  down  selection  within  a  few  reaction  wheels, 
subsequent  iterations  should  include  gyroscopic  coupling  in  the  sizing  analysis. 


2.  Market  Research  and  Simplified  Momentum  and  Torque  Equations 

Market  data  may  be  used  to  develop  a  relationship  between  momentum  and 
torque  for  reaction  wheel  and  CMG  systems.  Appendix  B  consists  of  a  table  of 
performance  data  collected  from  company  specifications  sheets  and  published  papers 
across  the  Internet.  Reaction  wheels  are  more  common  and  over  fifty  data  points  were 
readily  available.  CMGs  are  far  less  prevalent  as  less  than  twenty  systems  could  be 
identified.  Figure  77  offers  a  glimpse  at  this  momentum  and  torque  comparison  for 
reaction  wheels  and  CMGs.  Reference  [4,  Figure  4]  offers  a  similar  chart  and  the  six 
CMG  data  points  presented  in  this  paper  are  overlaid  on  the  chart  for  comparison  with  the 
information  presented  in  this  thesis.  The  market  data  from  Appendix  B  are  marked  with 
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“X’s.”  The  CMG  data  from  [4]  are  marked  with  circles.  Common  points  contain  both 
marks,  but  are  counted  only  once  to  facilitate  curve-fitting  of  the  data. 
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Figure  77.  Market  Research  Torque  vs.  Angular  Momentum  for  Reaction  Wheel  and 

CMG  Systems 

The  data  in  Figure  77  is  curve  fit  to  obtain  the  trend  for  each  type  of  system.  A 
few  things  become  clear.  First,  there  is  a  near  proportional  balance  between  CMG 
momentum  and  torque,  Equation  (266).  This  is  because  the  individual  CMG  torque  is 
determined  by  the  product  of  gimbal  rate  (commonly  1  rad/s  [27,  p.  8])  and  momentum. 
On  the  other  hand,  reaction  wheels  seem  to  have  an  approximately  square  root 
relationship  between  momentum  and  torque  with  a  proportionality  constant  of  less  than 
three  percent,  Equation  (267).  The  market  data  confirms  that  reaction  wheels  have  far 
less  torque  capability  than  CMGs. 

CMG:  llng  -  hllg  since  8  ~lrad/s  (266) 
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(267) 


Reaction  Wheel:  rZ  ~  0.027 

where  torque  units  are  N-m  and  momentum  units  are  N-m-s. 

Figure  77  is  created  on  a  per-wheel/CMG  basis.  This  means  that  MED  arrays  are 
not  accounted  for,  so  array  calculations  using  orientation  matrices  are  still  necessary. 
Another  notable  figure  is  the  large  torque  difference  at  high  momentum  levels.  This 
reveals  why  so  many  CMGs  are  being  developed  in  the  100+  Nms  range.  On  the  other 
hand,  micro-CMGs  are  rare,  seemingly  for  good  reason.  With  the  shape  of  the  curves, 
there  actually  may  be  a  torque  advantage  for  reaction  wheels  around  the  0.001  Nms 
momentum  level.  Regardless,  these  curves  offer  a  basis  for  comparison.  Therefore,  given 
momentum  storage  of  a  reaction  wheel  or  CMG,  the  torque  available  will  be  described  as 


CMG: 

„  !  n  \  1 .0019 

tg  =  1.0967  hG 

cmg  \  cmg ) 

(268) 

Reaction  Wheel: 

w  /  w  \0.5351 

%l  =  0.0270  (/£) 

(269) 

C.  COMPARISON  DEVELOPMENT 


For  comparison  purposes,  a  propagation  of  effective  eigenaxis  trajectories  for  the 
optimal  models  will  be  completed,  similar  to  Figures  51  to  53  from  Chapter  V.  The 
optimal  maneuvers  will  be  propagated  in  a  bang-bang  or  bang-coast-bang  trajectory 
about  its  effective  eigenaxis  (Chapter  V,  Section  H).  Already  derived  previously  in  this 
thesis,  the  following  spherically-available  momenta  will  be  used. 
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Ka  =1-333  h 


max 

w 

rvvm„Y 


K-pSeoPtA  =  1-577/1. 


hLoptA  =1-805  hlm 


KmgA  =  1 .007  hn 


(270) 


The  torque  equations  within  Equations  (268)  and  (269)  will  be  used  to  compute 
reaction  wheel  and  CMG  output  torques.  Transfonnation  of  torque  available  to  the  body 
frame,  given  a  4  reaction  wheel  or  4  CMG  array,  is  computed  as  follows. 
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=  1.333^ 

'  ,lmax 

.B 
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=  1 .007  xG 

'  ^max 

(271) 


For  1  Nms  wheels,  the  data  in  Table  3  can  be  determined  using  Equations  (270) 
and  (271). 


Wheel/Gimbal  Frame 
hmax  (Nms)  Tmax  (Nm) 

Body  Frame 
hmax  (Nms)  tmax  (Nm) 

RW 

1 

0.027 

1.333 

0.036 

RW  pseudo-opt 

1 

0.027 

1.577 

0.036 

RWopt 

1 

0.027 

1.807 

0.036 

CMG 

1 

1.097 

1.007 

1.104 

Table  3.  Body  Momentum/Torque  Comparison  Table  -  1  Nms  Momentum 


Now  a  spacecraft  size  class  must  be  established.  Since  the  MED  momentum 
storage  is  the  starting  point,  either  the  spacecraft  size  or  desired  slew  rate  must  be  chosen. 
Spacecraft  rate  will  be  chosen  next,  because  we  are  focusing  on  agile  spacecraft.  In  the 
space  industry,  a  2-3  °/s  slew  is  considered  agile,  even  for  CMGs.  For  example, 
Worldview  1  slews  at  2.3°/s  [2].  Therefore,  a  rate  of  NcoB  =  3 °  /  s  will  be  the  benchmark 
with  which  the  spacecraft  sizing  will  be  detennined.  The  base  reaction  wheel  will  be  used 
for  sizing,  and  all  other  systems  will  be  required  to  slew  the  same  spacecraft. 


/  B  ,  B  ■  N 

™ rwA  —  'T-req  ~  J sic  ® 
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N,.B 

CO 


1 333  Nms 
3  °/s 


=  25.5  kg  ■  in 


(272) 


Using  these  properties,  the  individual  spacecraft  rates  for  the  various  systems  can 
be  solved  for. 
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(273) 


(274) 


The  data  from  equations  (268)  through  (274)  are  compiled  in  Table  4.  Notice  the 
higher  momentum  capacity  for  reaction  wheels  corresponds  directly  to  the  larger 
maximum  spacecraft  rotation  rates.  Note  that  this  rotation  rate  only  refers  to  the  effective 
rotation  rate  about  the  eigenaxis.  The  actual  rotation  rate  of  the  optimal  case  about  the 
instantaneous  axis  of  rotation  can  potentially  be  much  faster  so  the  values  in  Table  4  are 
conservative  estimates.  The  body  acceleration  of  the  reaction  wheels  is  more  than  an 
order  of  magnitude  lower  than  the  acceleration  possible  using  the  CMGs. 


Wheel/Gimbal  Frame 
hmax  (Nms)  imax  (mNm) 

S/C 

js/c  (kg-m2) 

Body  Frame 

hmax  (Nms)  imax  (Nm)  w  (%)  a  (°/s2) 

RW 

1  0.027 

25.458 

1.333  0.036  3.000  0.081 

RW 

’  ~  pseudo-opt 

1  0.027 

25.458 

1.577  0.036  3.549  *  0.081  * 

RW„P, 

1  0.027 

25.458 

1.807  0.036  4.067  *  0.081  * 

CMG 

1  1.097 

25.458 

1.007  1.104  2.266  2.485 

*effective  rotation  and  acceleration  rates  given  about  an  effective  eigenaxis.  Refer  back  to  Section  H  of  Chapter  6  for 
description  of  effective  eigenaxis.  The  actual  body  rates  and  acceleration  can  be  much  higher  throughout  the  maneuver. 

Table  4.  Body  Rate/Acceleration  Comparison  for  a  25.5  kg-m2  Class  Spacecraft 
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Given  the  effective  body  rate  and  acceleration  data  from  Table  4,  the  data  can 
now  be  propagated  through  the  “effective”  or  equivalent  eigenaxis  slews  using 
trapezoidal  input  trajectories.  This  is  the  same  as  what  was  done  for  the  optimal  reaction 
wheel  case  from  Chapter  V.  Recall  the  curve  from  Chapter  V  was  built  using  a  specific 
reaction  wheel  of  0. 1  Nm  torque,  much  higher  than  the  one  used  here  per  the  relationship 
in  Equation  (269).  Therefore,  the  curves  are  built  using  different  acceleration  data,  but  the 
concept  is  the  same. 

The  plot  of  the  eigenaxis  slew  angle  versus  slew  time  for  all  four  cases  is  shown 
in  Figure  78.  This  figure  demonstrates  the  shape  of  the  curves.  Since  the  reaction  wheel 
arrays  have  higher  momentum  storage,  the  time/angle  slope  of  the  momentum-limited 
region  is  much  smaller  than  the  one  for  the  CMG  case.  This  causes  the  reaction  wheel 
systems  to  eventually  catch  up  with  the  CMG  system.  The  time-optimal  case  would  not 
catch  up  with  the  CMG  until  the  maneuver  exceeds  180°,  which  is  an  unrealistic  slew. 
For  a  25.5  kg-m2  class  spacecraft,  the  control  system  comparison  shows  CMGs  to  be 
more  favorable  for  agile  maneuvering. 
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Now,  suppose  the  wheel  momentum  is  decreased  by  an  order  of  magnitude 
( hrmx  =  0.1  Nms  ).  Using  the  same  N(0B  =  3 °  / s  requirement,  the  corresponding  spacecraft 
is  now  2.55  kg-m2.  The  new  comparison  data  is  shown  in  Table  5. 


Wheel/Gimbal  Frame 
hmax  (Nms)  imax  (Nm) 

S/C 

js/c  (kg-m2) 

Body  Frame 

hmax  (Nms)  imax  (Nm)  w  (°/s) 

a  (°/s2) 

RW 

0.1 

0.008 

2.546 

0.133  0.010 

3.000 

0.236 

RW 

T  ’  pseudo-opt 

0.1 

0.008 

2.546 

0.158  0.010 

3.549  * 

0.236  * 

RWopt 

0.1 

0.008 

2.546 

0.181  0.010 

4.067  * 

0.236  * 

CMG 

0.1 

0.109 

2.546 

0.101  0.110 

2.266 

2.474 

*effective  rotation  and  acceleration  rates  given  about  an  effective  eigenaxis.  Refer  back  to  Section  H  of  Chapter  6  for 
description  of  effective  eigenaxis.  The  actual  body  rates  and  acceleration  can  be  much  higher  throughout  the  maneuver. 

Table  5.  Body  Rate/Acceleration  Comparison  for  a  2.55  kg-m2  Class  Spacecraft 


Notice  the  spacecraft  rates  did  not  change  across  the  two  models.  Instead,  the 
spacecraft  acceleration  is  the  important  change.  Due  to  the  roughly  1:1  ratio  CMG 
momentum  to  torque  ratio  (  r ~  h™  )  from  Equation  (268),  the  CMG  trajectory  did  not 

change  very  much.  On  the  other  hand,  the  reaction  wheel  has  a  square  root  relationship, 
( T^,  -0.027 -JO™,  )  from  (269).  This  causes  the  reaction  wheel  torque  to  momentum  ratio 

to  improve  as  the  wheel  size  is  decreased.  The  reaction  wheel  acceleration  rates  are  now 
about  three  times  higher  than  for  the  larger  wheel  (see  Table  4).  Generating  the  effective 
eigenaxis  trajectories  for  each  system  gives  Figure  79.  Using  these  curves,  it  becomes 
evident  that  any  slews  greater  than  77°  will  favor  the  time-optimal  reaction  wheel  control 
system  over  the  CMG. 
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Notice  the  significant  change  in  the  reaction  wheel  curves.  The  spacecraft  rates 
did  not  change  across  the  two  models,  but  the  improved/larger  spacecraft  acceleration 
obtained  from  the  reaction  wheel  system  increased  by  about  a  factor  of  three.  Because  the 
angular  rates  are  the  same  as  before,  the  time/angle  slope  of  the  momentum-limited 
region  did  not  change  at  all.  Due  to  the  higher  available  torque,  the  momentum-limited 
region  starts  at  a  lower  angle  for  the  reaction  wheels.  This  causes  an  earlier  intersection 
with  the  CMG  curve.  To  better  illustrate  these  transition  points,  the  reaction  wheel  curves 
are  looked  at  as  a  function  of  angle,  after  which  the  performance  of  the  reaction  wheel 
assembly  exceeds  that  of  the  CMG  assembly.  This  threshold  will  be  referred  to  for  the 
remainder  of  the  thesis  as  the  reaction  wheel  advantage  angle,  shown  in  Figure  80. 


148 


Figure  80.  Maneuver  Angles  where  Reaction  Wheel  Systems  Outperform  CMGs  for  a 

2.55  kg-m2  Class  Spacecraft 

The  trend  of  rapidly  increasing  reaction  wheel  performance  is  continued  when 
wheel  momentum  is  decreased  by  another  order  of  magnitude  (hmm=0.0\Nms). 

Using  the  same  N(Ob  =  3°/s  requirement,  the  corresponding  spacecraft  size  is  now 
0.255  kg-m2.  The  comparison  is  shown  in  Table  6  and  Figures  83-84,  illustrating  that 
the  reaction  wheel  advantage  angles  decrease  as  body  accelerations  provided  by  the 
wheels  increase. 


Wheel/Gimbal  Frame 
hmax  (Nms)  imax  (mNm) 

S/C 

js/c  (kg-m2) 

Body  Frame 

hmax  (Nms)  Tmax  (Nm)  w  (%) 

Ct  (7s2) 

RW 

0.01 

0.002 

0.255 

0.013  0.003 

3.000 

0.689 

RW 

*  ’  pseudo-opt 

0.01 

0.002 

0.255 

0.016  0.003 

3.549  * 

0.689  * 

RWopt 

0.01 

0.002 

0.255 

0.018  0.003 

4.067  * 

0.689  * 

CMC 

0.01 

0.011 

0.255 

0.010  0.011 

2.266 

2.464 

*effective  rotation  and  acceleration  rates  given  about  an  effective  eigenaxis.  Refer  back  to  Section  H  of  Chapter  6  for 
description  of  effective  eigenaxis.  The  actual  body  rates  and  acceleration  can  be  much  higher  throughout  the  maneuver. 

Table  6.  Body  Rate/Acceleration  Comparison  for  a  0.255  kg-m2  Class  Spacecraft 
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CMG  RW 


Figure  81.  Slew  Angle  vs.  Time  Comparison  for  a  0.255  kg-m2  Class  Spacecraft 
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Figure  82.  Maneuver  Angles  where  Reaction  Wheel  Systems  Outperform  CMGs  for  a 

0.255  kg-m2  Class  Spacecraft 


Referring  to  Figures  78,  79,  and  81,  it  is  apparent  that  when  using  like-sized 

reaction  wheels  and  CMGs,  the  larger  momentum  capacity  of  reaction  wheels  is  much 
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more  advantageous  for  small  spacecraft  than  for  larger  spacecraft.  In  the  case  shown  in 
Figure  83,  the  optimal  reaction  wheel  control  system  becomes  more  agile  than  the  CMG 
control  system  for  all  maneuvers  greater  than  24.1°.  Standard  operation  of  the  reaction 
wheel  array  does  not  begin  to  outperfonn  CMGs  until  the  maneuver  size  exceeds  30.3°. 
If  only  the  reaction  wheel  advantage  angle  is  plotted  against  the  wheel  momenta,  the 
curves  in  Figure  83  are  obtained.  The  curves  reveal  that  reaction  wheels  overcome  CMGs 
at  relatively  small  angles  for  small  rotor  momenta.  As  the  spacecraft  momentum 
increases,  the  reaction  wheel  advantage  angle  increases  rapidly  and  reaction  wheel 
control  systems  promptly  lose  their  edge. 


Figure  83.  Reaction  Wheel  Advantage  Angle  vs.  Wheel/Rotor  Momentum:  Where  RWs 

Overcome  CMGs  (co=3°/s) 


Likewise,  if  the  spacecraft  inertia  versus  slew  angle  is  plotted  (assuming  co=3°/s), 
the  trend  is  shown  again  (Figure  84).  For  small  spacecraft,  reaction  wheels  seem  to 
outperform  CMGs,  especially  for  large  angle  maneuvers. 
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Figure  84.  Reaction  Wheel  Advantage  Angle  vs.  Spacecraft  Inertia:  Where  RWs 

Overcome  CMGs  (co=3°/s) 

Thus  far,  the  results  of  this  chapter  show  that  reaction  wheels  outperform  CMGs 
for  small  spacecraft.  In  addition,  the  two  optimal  control  approaches  described 
significantly  outperformed  the  heritage  eigenaxis  maneuver.  For  small  spacecraft,  the 
reaction  wheel  advantage  angle  was  reduced  by  only  a  few  degrees  by  using  optimal 
reaction  wheel  methods  (Figure  85).  At  larger  slew  angles,  40-60°  improvement  can  be 
achieved  by  optimal  systems  over  the  heritage  eigenaxis  maneuver  (Figure  85).  This 
means  for  large  angle  slews,  the  optimal  reaction  wheel  system  outperforms  the  CMG 
system  significantly  more  quickly  than  if  heritage  control  logic  were  used. 
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Figure  85.  Zoom  in:  Reaction  Wheel  Advantage  Angle  vs.  Spacecraft  Inertia:  Where 

RWs  Overcome  CMGs  (co=3°/s) 


Figure  86  shows  these  reductions  in  reaction  wheel  advantage  angles  as  a 
percentage  improvement.  The  pseudo-inverse  optimal  method  reduces  the  reaction  wheel 
advantage  angle  by  roughly  15%  over  the  heritage  eigenaxis  shaped  trajectory,  while  the 
full  time-optimal  solution  reduces  the  reaction  wheel  advantage  angle  by  18%  across  the 
entire  spacecraft  inertia  regime. 
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Figure  86.  Optimal  Improvement  over  normal  Reaction  Wheel  Systems 


D.  MOMEMTUM  EXCHANGE  DEVICE  PERFORMANCE  AS  A  FUNCTION 

OF  MASS 

In  the  analysis  of  Section  C,  equally  sized  reaction  wheels  and  CMGs  were  used, 
i.e.,  both  had  the  same  momentum  storage  capacity.  This  does  not  have  to  be  the  case. 
Suppose  the  systems  engineer  gives  a  mass  restriction  for  the  entire  MED  array  to  the 
ADCS  engineer.  The  ADCS  engineer  can  quickly  determine  whether  reaction  wheels  or 
CMGs  would  be  better  given  this  limitation.  In  this  section,  it  will  be  shown  how  much 
more  momentum  can  come  of  a  reaction  wheel  array  in  comparison  with  CMGs  given  a 
mass  constraint.  Figure  87  shows  the  market  research  for  individual  MED  mass  and  the 
output  momentum  the  actuator  can  achieve.  The  data  was  curve  fit  to  develop  a  general 
trend. 
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Figure  87.  Market  Research  Angular  Momentum  vs.  Mass 


The  mechanical  complexity  of  CMGs  requires  additional  mass  and  volume  within 
the  spacecraft.  Using  this  knowledge  and  the  relationships  from  Section  C,  the  full  array 
mass  and  body  momentum  can  be  created  using  the  trend  lines  (Figure  88).  This  analysis 
shows  that  for  the  same  array  mass,  an  order  of  magnitude  more  momentum  storage  is 
typically  available  from  reaction  wheel  arrays  than  from  CMG  arrays.  This  implies  that  if 
mass  is  a  critical  design  constraint  for  the  spacecraft,  then  engineers  should  take  an  even 
harder  look  at  reaction  wheels. 
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Array  Mass  (kg) 


Figure  88.  Array  Momentum  vs.  Mass 


Torque  output,  on  the  other  hand,  is  better  for  large  mass  CMGs.  Figure  89 
demonstrates  the  transition  where  a  CMG  array  will  begin  to  create  more  torque  than  a 
reaction  wheel  array.  This  is  an  important  point  because  it  helps  illustrate  why  few  CMGs 
have  been  developed  on  the  small  scale.  From  the  market  research  shown  in  Figure  87, 
only  one  CMG  was  actually  smaller  than  0.898  kg. 
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Array  Mass  (kg) 


Figure  89.  Array  Torque  vs.  Mass 


E.  SUMMARY 

In  this  chapter,  it  was  shown  that  reaction  wheels  seem  to  perfonn  better  than 
CMGs  for  small  spacecraft.  As  the  spacecraft  size  decreases,  the  reaction  wheel 
advantage  increases.  This  increase  can  be  further  improved  by  the  use  of  optimal  control 
algorithms. 

The  CMG  torque  advantage  is  evident  for  large  moments  of  inertia.  On  the  other 
hand,  the  torque  gap  lessens  as  the  spacecraft  size  (therefore  the  MED  array  momentum 
capacity)  decreases.  Although  the  gap  does  not  completely  close,  the  ratio  of  CMG 
torque  to  reaction  wheel  torque  decreases  significantly  with  spacecraft  size. 

Given  arrangements  of  reaction  wheels  and  CMGs  with  similar  individual  MED 
momentum,  the  reaction  wheel  array  will  have  a  larger  momentum  storage  capacity.  This 
causes  reaction  wheel  systems  to  catch  up  with  CMG  systems  over  large  angle  slews.  The 
reaction  wheel  advantage  angle  is  the  slew  size  in  which  a  reaction  wheels  become  a 
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better  alternative  than  CMGs.  For  large  spacecraft,  the  reaction  wheel  advantage  angle 
angle  is  very  high,  possibly  over  180  degrees. 

Time-optimal  shaped  maneuvers  of  reaction  wheel  systems  can  be  used  to  lower 
this  reaction  wheel  advantage  angle  by  13  to  21%,  depending  upon  algorithms  used  and 
spacecraft  size.  By  shaping  only  the  quaternion  and  body  rate  trajectories,  the 
mathematical  limitations  of  the  pseudo-inverse  control  allocation  must  still  be  ascribed 
to.  This  allows  a  reduction  of  the  time-optimal  reaction  wheel  advantage  angle  by 
approximately  15%  over  the  advantage  angle  of  the  shaped  eigenaxis  maneuver.  The 
addition  of  shaped  time-optimal  wheel  torques  reduces  the  advantage  angle  by 
approximately  18%. 

Considering  the  mass  of  the  reaction  wheels  and  CMGs,  the  reaction  wheels 
produce  a  much  higher  momentum  limit  than  CMGs  for  a  given  mass.  Interestingly,  the 
additional  order  of  magnitude  of  momentum  equates  to  additional  torque  in  the  similar- 
weight  MEDs.  For  MEDs  smaller  than  0.898  kg,  the  reaction  wheel  array  will  actually 
produce  more  torque  than  a  CMG  array.  This  means  that  given  a  small  spacecraft  with 
weight  restrictions,  reaction  wheels  are  the  better  alternative  for  both  torque  and 
momentum. 
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VIII.  CONCLUSIONS  AND  FUTURE  WORK 


A.  CONCLUSIONS 

This  thesis  is  an  attempt  to  bridge  the  gap  between  reaction  wheels  and  CMGs, 
particularly  for  small  satellites.  Using  the  software  upgrades  and  relative  sizing  concepts 
proposed,  reaction  wheel  arrays  have  the  ability  to  outperform  CMGs  for  some  classes  of 
spacecraft.  The  relatively  low  cost  of  reaction  wheels  can  be  attractive  once  again. 

Reaction  wheels  are  typically  used  inefficiently  within  the  spacecraft  industry. 
The  standard  maneuver,  a  shaped  bang-coast-bang  maneuver  about  the  eigenaxis,  should 
be  reconsidered  in  the  presence  of  new  control  algorithms.  The  results  of  this  thesis  show 
that  the  speed  of  this  eigenaxis-restricted  maneuver  can  be  improved  18%  in  each  of 
momentum  and  torque  by  applying  the  pseudo-inverse  limited  time-optimal  shaped 
quaternion  and  rate  trajectories.  Furthermore,  the  speed  of  the  standard  eigenaxis 
maneuver  can  be  improved  by  36%  by  delivering  the  optimal  torque  inputs  directly  to  the 
reaction  wheels  themselves.  These  two  software  solutions  reduce  the  reaction  wheel 
advantage  angle  by  15%  and  18%,  respectively.  The  CMG  performance  gap  is  closed  and 
surpassed  in  many  cases. 

Due  to  their  high  torque  capacity,  CMGs  dominate  reaction  wheels  for  small 
angle  maneuvers.  This  is  because  a  spacecraft  using  CMGs  can  accelerate  to  maximum 
angular  velocity  very  quickly.  Although  CMGs  have  a  lot  more  torque  potential,  they  are 
exceedingly  momentum-limited.  A  spacecraft  using  reaction  wheels,  on  the  other  hand,  is 
slower  to  accelerate,  but  has  a  much  greater  maximum  velocity.  Considering  like-sized 
wheels,  i.e.,  the  same  momentum  capacity,  this  thesis  shows  that  reaction  wheels  catch  up 
and  surpass  CMGs  due  to  their  higher  rate  limits,  provided  the  maneuver  is  large  enough. 
The  ramifications  are  particularly  large  for  smallsats  and  cubesats  as  they  require  very 
small  MEDs.  For  small  satellites,  this  thesis  has  shown  that  reaction  wheels  have  the 
ability  to  outperform  like-sized  CMGs  at  relatively  small  angle  maneuvers.  As  wheel  and 
spacecraft  sizes  decrease  proportionally,  reaction  wheels  catch  up  to  CMGs  even  more 
quickly.  Optimal  control  further  improves  this  performance. 
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B.  FUTURE  WORK 

This  thesis  did  not  include  an  exhaustive  approach  at  validating  the  models  of  the 
analytical  effective  eigenaxis  of  each  of  the  optimal  control  models.  Only  100  random 
data  points  were  used  for  a  single  combination  of  spacecraft/wheel  sizing.  The  results  of 
this  work  should  be  validated  across  the  entire  space  of  reaction  wheel  and  spacecraft 
sizes.  Gimbal  momentum/torque  and  gyroscopic  coupling  relationships  should  be 
accounted  for  into  this  “rule  of  thumb.”  For  instance,  as  body  rate  increases,  the 
gyroscopic  coupling  increases.  Most  likely,  the  optimal  maneuvers  will  use  this  as  an 
advantage  to  widen  the  perfonnance  gap  over  standard  maneuvers. 

The  reaction  wheel  and  CMG  data  should  be  analyzed  further  to  account  for 
variations  in  momentum  size  and  torque  output.  To  further  understand  the  CMG 
momentum  and  torque  relationship,  research  must  be  conducted  to  determine  why  a 
1 : 1  relationship  between  momentum  and  output  torque  is  so  commonly  used. 

The  different  control  approaches  should  also  be  explored  from  non-rest  initial  and 
final  conditions.  Given  some  calculable  motion  within  a  system,  the  optimal  methods 
should  be  able  to  capitalize  on  this  motion  for  further  gains  over  the  heritage  maneuvers. 

Optimal  control  of  the  CMG  should  be  explored.  The  optimal  maneuver,  by  its 
very  nature,  should  be  able  to  avoid  the  singularities  altogether.  This  would  allow  the 
unused  capacity  of  the  CMG  momentum  space  to  be  accessed  without  complicated 
singularity-avoidance  algorithms.  One  consideration  during  this  expansion  would  be  the 
confidence  level  in  the  spacecraft  inertia  and  states.  If  knowledge  of  these  uncertainties  is 
poor,  then  the  optimal  control  approach  would  still  have  a  pseudo-inverse  limited  and 
singularity-ridden  feedback  system  to  deal  with.  Therefore,  methods  such  as  those 
described  in  [31]  could  be  tried. 

Further  analysis  should  be  done  on  the  momentum  and  torque  envelope 
visualization  method  for  optimal  control  maneuvers.  Looking  at  these  envelopes  closely, 
it  becomes  obvious  that  the  torque  and  momentum  vectors  tend  toward  the  vertices  and 
edges  of  the  envelopes.  If  this  space  is  studied  in  detail,  it  should  be  possible  to  predict 
the  optimal  slew  times  much  more  accurately.  Just  as  the  optimal  eigenaxis-restricted 
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maneuver  is  a  bang-coast-bang  optimal  path,  certain  characteristics  of  optimal  maneuvers 
could  be  derived  for  maneuver  size  and  eigenaxis  orientation.  This  could  allow  quick  on¬ 
board  processing  of  semi-optimal  maneuvers,  without  actually  computing  the  exact 
computer-intensive  optimal  path.  A  controller  much  like  the  eigenaxis  quaternion 
feedback  could  potentially  be  developed  to  drive  the  system  toward  the  appropriate 
vertices  and  edges  of  the  momentum  and  torque  envelopes. 

A  comparison  should  be  done  with  various  reaction  wheel  and  CMG 
arrangements.  For  example,  if  it  is  known  that  the  pseudo-inverse  limited  time-optimal 
control  logic  is  going  to  be  used,  then  it  may  be  possible  to  adjust  the  skew  angle  to 
maximize  its  envelope  volume,  thereby  creating  a  larger  momentum/torque  radius  (per 
Chapter  V,  Section  G).  Additionally,  sizing  wheels  with  different  arrangements  of  wheel 
sizes  could  be  considered.  For  example,  a  skew  arrangement  of  five,  six,  or  possibly  ten 
smaller  wheels  may  be  more  capable  than  four  larger  wheels,  not  to  mention  more 
redundancy.  The  shape  would  have  more  flat  facets,  but  less  loss.  Like  a  cube  and  a 
regular  octahedron  of  similar  volumes,  the  octahedron  allows  a  larger  inscribed  sphere. 
Mixed  configurations  of  larger  and  smaller  wheels  could  also  be  considered.  Finally, 
hybrid  reaction  wheel  and  CMG  configurations  could  be  considered  to  maximize  torque, 
yet  maintain  a  high  momentum  envelope. 

Explore  the  impact  of  external  disturbances  on  optimal  maneuvers.  Integrate 
gravity  gradient,  solar  pressure,  and  aerodynamic  torques  into  the  optimal  control 
problem.  Compare  with  the  eigenaxis  maneuvers.  Most  likely,  the  optimal  maneuvers 
will  be  able  to  exploit  the  external  disturbances  to  widen  the  performance  gap  between 
the  shaped  eigenaxis  trajectory  and  the  time-optimal  trajectories. 
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APPENDIX  A.  OPTIMAL  CONTROL  PROBLEM  IN  DESIGNER 
UNITS 


Minimize:  3[jc(-), «(•)>*/]  =  tf  where  3  =  3T 

=>  3[m,mJfy  =  tfr  (275) 
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APPENDIX  B.  MARKET  RESEARCH 


A.  REACTION  WHEELS 


Manufacturer 

Model 

Mass 

(kg) 

**max 

(N-m-s) 

T 

*  max 

(N-m) 

Reference 

Bradford  Engineering 

W18 

4.95 

18 

0.248 

[  32  ] 

Bradford  Engineering 

W18E 

5.5 

22 

0.248 

[  32  ] 

Bradford  Engineering 

W18ES 

5.85 

25 

0.248 

[  32  ] 

Bradford  Engineering 

W45 

6.45 

40 

0.248 

[  32  ] 

Bradford  Engineering 

W45HT 

6.5 

22 

0.403 

[  32  ] 

Bradford  Engineering 

W45E 

6.65 

55 

0.248 

[  32  ] 

Bradford  Engineering 

W45ES 

6.86 

70 

0.248 

[  32  ] 

Dynacon  MSCI 

Micro  Wheel  4000 

5.2 

4 

0.05 

[  33  ] 

Dynacon  MSCI 

Micro  Wheel  1000 

1.4 

1.1 

0.03 

[  34  ] 

Dynacon  MSCI 

Micro  Wheel  200 

1 

0.18 

0.03 

[  35  ] 

Honeywell 

HR0610 

5 

12 

0.055 

[  36  ] 

Honeywell 

HR0610 

3.6 

4 

0.055 

[  36  ] 

Honeywell 

HR12-12 

6 

12 

0.1 

[  37  ] 

Honeywell 

HR  12-25 

7 

25 

0.15 

[  37  ] 

Honeywell 

HR  12-50 

9.5 

50 

0.2 

[  37  ] 

Honeywell 

HR  14-25 

7.5 

25 

0.1 

[  37  ] 

Honeywell 

HR  14-50 

8.5 

50 

0.15 

[  37  ] 

Honeywell 

HR  14-75 

10.6 

75 

0.2 

[  37  ] 

Honeywell 

HR  16-50 

6 

50 

0.1 

[  37  ] 

Honeywell 

HR  16-75 

7 

75 

0.15 

[  37  ] 

Honeywell 

HR  16- 100 

9.5 

100 

0.2 

[  37  ] 

L3 

RWA-15 

14 

20 

0.68 

[  38  ] 

L3 

MWA-50 

10.5 

67.8 

0.07 

[  39  ] 

Maryland  Aerospace 

MAI- 101 

0.64 

0.0011 

0.000635 

[  40  ] 

Maryland  Aerospace 

MAI-201 

0.73 

0.0108 

0.000625 

[  41  ] 

Maryland  Aerospace 

MAI-300 

0.317 

0.0076 

0.000625 

[  42  ] 

NASA  GSFC 

SMEX  RW 

3.3 

4 

0.14 

[  43  ] 

Table  7.  Market  Research  Reaction  Wheels 
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Manufacturer 

Model 

Mass 

(kg) 

frmax 

(N-m-s) 

T 

''max 

(N-m) 

Reference 

Rockwell  Collins 

RSI  01-5/15 

0.6 

0.04 

0.005 

[  44  ] 

Rockwell  Collins 

RSI  01-5/28i 

0.7 

0.12 

0.005 

[  44  ] 

Rockwell  Collins 

RSI  12-75/60 

4.85 

12 

0.075 

[  44  ] 

Rockwell  Collins 

RSI  4-75/60 

3.7 

4 

0.075 

[  44  ] 

Rockwell  Collins 

RSI  15-215/20 

7.7 

15 

0.215 

[  44  ] 

Rockwell  Collins 

RSI  45-75/60 

7.7 

45 

0.075 

[  44  ] 

Rockwell  Collins 

RSI  12-220/45 

6 

12 

0.22 

[  44  ] 

Rockwell  Collins 

RSI  18-220/45 

6.45 

18 

0.22 

[  44  ] 

Rockwell  Collins 

RSI  25-220/45 

7.15 

25 

0.22 

[  44  ] 

Rockwell  Collins 

RSI  50-220/45 

9.2 

50 

0.22 

[  44  ] 

Rockwell  Collins 

RSI  30-280/30 

9.2 

30 

0.28 

[  44  ] 

Rockwell  Collins 

RSI  68-170/60 

9.5 

68 

0.17 

[  44  ] 

Rockwell  Collins 

RSI  68-170/60 

9.5 

68 

0.17 

[  44  ] 

Rockwell  Collins 

RDR  23-0 

5.5 

23 

0.09 

[  44  ] 

Rockwell  Collins 

RDR  57-0 

7.6 

57 

0.09 

[  44  ] 

Rockwell  Collins 

RDR  68-3 

7.6 

68 

0.075 

[  44  ] 

Rockwell  Collins 

MWI  100-100/100 

16.5 

100 

0.1 

[  44  ] 

Rockwell  Collins 

MWI  30-400/37 

15.3 

30 

0.4 

[  44  ] 

Sinclair 

Picosatellite  RW 

0.09 

0.007 

0.001 

[  45  ] 

Sinclair 

Picosatellite  RW 

0.12 

0.01 

0.001 

[  45  ] 

Sinclair 

Nanosatellite  RW 

0.185 

0.03 

0.002 

[  45  ] 

Sinclair 

Microsatellite  RW 

0.225 

0.06 

0.005 

[  45  ] 

Sinclair 

Microsatellite  RW 

0.225 

0.06 

0.025 

[  45  ] 

Sinclair 

Microsatellite  RW 

0.97 

1 

0.1 

[  45  ] 

SunSpace 

SunStar 

1.986 

0.65 

0.05 

[  46  ] 

Surrey  Space  SSTL 

Micro  Wheels  10SP-M 

0.96 

0.42 

0.011 

[  47  ] 

Surrey  Space  SSTL 

Micro  Wheels  100SP-0 

2.6 

1.5 

0.11 

[  47  ] 

Surrey  Space  SSTL 

Smallsat  Wheels  200SP-M/O 

5.2 

12 

0.24 

[  48  ] 

Surrey  Space  SSTL 

Smallsat  Wheels  200SP-0 

5.2 

12 

0.24 

[  48  ] 

Teldix 

20  cm 

3.4 

6.5 

0.2 

[  5,  p.  398] 

Teldix 

26  cm 

6 

20 

0.2 

[  5,  p.  398] 

Teldix 

35  cm 

8 

80 

0.2 

[  5,  p.  398] 

Teldix 

50  cm 

12 

300 

0.3 

[  5,  p.  398] 

Teldix 

60  cm 

37 

1000 

0.6 

[  5,  p.  398] 

Table  8.  Market  Research  Reaction  Wheels  (cont.) 
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B. 


CMGS 


Manufacturer 

Model 

Mass 

(kg) 

1_ 

**max 

(N-m-s) 

r 

•'max 

(N-m) 

• 

Q 

"max 

(rad/s) 

Reference 

Astrium 

CMG  4-6S 

13 

4 

6 

1.5 

L  49  ] 

Astrium 

CMG  15-45S 

18.4 

15 

45 

3 

[  50  ] 

Astrium 

CMG  10-30 

18.4 

10 

30 

3 

[  51  ] 

Goodrich/Ithaco 

Violet  CMG 

1.28 

0.31 

0.3100 

1 

[  52  ] 

Honeywell 

M50 

33.1 

75 

75 

1 

L  53  ] 

Honeywell 

M95 

38.6 

129 

129 

1 

[  53  ] 

Honeywell 

M160 

44 

217 

217 

1 

[  53  ] 

Honeywell 

M225 

54 

305 

305 

1 

L  53  ] 

Honeywell 

M325 

61.2 

441 

441 

1 

L  53  ] 

Honeywell 

M325D 

61.2 

441 

441 

1 

[  53  ] 

Honeywell 

M715 

89.8 

969 

969 

1 

[  53  ] 

Honeywell 

M600 

81.6 

813 

813 

1 

[  53  ] 

Honeywell 

Ml  300 

125 

2000 

2000 

1 

[  53  ] 

Honeywell 

Ml  400 

132 

2000 

2000 

1 

[  53  ] 

Surrey  Space  SSTL 

BILSAT-1  CMG 

2.2 

0.28 

0.0559 

0.20 

L  54  ] 

Swampsat 

Swampsat 

0.5 

8.00E-04 

8.00E-04 

1 

[  55  ] 

Table  9.  Market  Research  Control  Moment  Gyroscopes 
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